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The quantum-mechanical treatment previously applied 
to benzene, naphthalene, and the hydrocarbon free radicals 
is used in the calculation of extra resonance energy of 
conjugation in systems of double bonds, the dihydro- 
naphthalenes and dihydroanthracenes, phenylethylene, 
stilbene, isostilbene, triphenylethylene, tetraphenylethyl- 


ene, biphenyl, 0, m, and p-diphenylbenzene, and 1,3,5- 
triphenylbenzene. The calculated values, which are in 
approximate agreement with empirical values from thermo- 
chemical data, are used in the discussion of chemical 
properties and in the formulation of rules regarding 
conjugation in overlapping systems. 


ITH the aid of the recently developed 

method! of calculating matrix elements 
for various electronic structures of molecules, 
based on Slater’s treatment,? we have carried out 
a quantum-mechanical discussion of a number of 
conjugated systems involving double bonds and 
aromatic nuclei. This discussion, which is pre- 
sented in the following pages, is closely similar to 
that previously given benzene and anthracene 
and the hydrocarbon free radicals? and has 
similarly led to an illuminating qualitative and 
semi-quantitative explanation of the properties of 
the systems. 


ADDITION Propucts oF NAPHTHALENE 


The measurements of heats of combustion 
made by Roth and Auwers‘ show that liquid 


* Contribution No. 370. 

‘Linus Pauling, J. Chem. Phys. 1, 280 (1933). The same 
calculations could be made by the use of the method of 
Eyring and Kimball, J. Chem. Phys. 1, 239 (1933). 

*J. C. Slater, Phys. Rev. 38, 1109 (1931). 

*Linus Pauling and G. W. Wheland, J. Chem. Phys. 1, 
362 (1933), 


Pt A. Roth and K. v, Auwers, Liebig’s Ann. 407, 145 


1,2-dihydronaphthalene, in which the benzene 
ring is conjugated with a double bond, is 4.9 
kcal. more stable than liquid 1,4-dihydronaph- 
thalene, in which this conjugation does not occur. 
It is probable that these isomers differ by about 
0.1 kcal. in heat of vaporization, the boiling 
points being 206.5 and 212.0°C, respectively, 
leading to a difference in energy of the gaseous 
molecules of 5.0 kcal. This information permits us 
to discuss the effect of change in carbon-carbon 
distance on the exchange integral between ad- 
jacent atoms. A canonical set of fourteen inde- 
pendent structures for 1,2-dihydronaphthalene 
(in all of which single bonds between adjacent 
atoms are left unchanged) is shown in Fig. 1. Of 
these the two unexcited structures 1 and 2 will 
contribute most to the normal state of the 
molecule, so that the carbon-carbon distance for 
the 3-4 bond and for the six bonds of the 
unbroken benzene ring will be the double-bond 
distance of 1.38A, and the corresponding ex- 
change integral will be that previously® desig- 
nated by the symbol a, and found to have a value 
of about —1.5 v.e. The carbon-carbon distance 
connecting the 3-4 bond with the benzene ring 
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will be somewhat larger, lying somewhere be- 
tween 1.38 and 1.54A, inasmuch as the structures 
which make this bond a double bond are first- 
excited structures, and hence play a subordinate 
réle in the normal state of the molecule. The 
corresponding exchange integral, which we may 
designate as 8, will thus have a somewhat 
different value from a (presumably smaller in 
absolute value). The remaining three carbon- 
carbon distances will have the single-bond value 
1.54A, for neither the unexcited nor the first- 
excited structures make them double bonds. The 
corresponding exchange integral 6’ will hence 
differ from a by an amount somewhat larger than 
for 

If we assume that structures 1 and 2 have the 
same coefficients in the wave function for the 
normal state, and similarly 3, 4, 5; 6, 7, 8; 
9, 10, 11; and 12, 13, 14, we obtain the following 
secular equation, in which x represents the 
quantity (Q—W’)/a, and b the ratio B/a. The 
energy W,,2 of the molecule is equal to W’ — 3/28, 
the latter term being omitted in the secular 
equation for convenience. 


5x/2+17/2—5b/4 3x-+12—3b/2 
3x+12—3b/2 


7x/4+25/4—b/8 9x/4+33/4 
x+4+b 


Because of the small contribution of the second 
and third-excited states, we may neglect them, 
which reduces the equation to a cubic, of the 
form 


b*(1982— 63) +5(2342?+ 1322—390) 
+ (682*+ 442? —8842—582)=0, (1) 


in which z= —x—b. 

If as an approximation we assume )=1(6= 8’ 
=a), we find as the pertinent root of this 
equation z=2.3062, corresponding to 
+1.8062a. A similar treatment applied to 1,4- 
dihydrobenzene (which involves no interaction 
between the benzene ring and the double bond) 
leads to Wi4=Q+3.6055a— 28’, or, in this case, 
Q+1.6055a. (In this paper as in the previous one 
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7x/4+25/44+7b/4 7x/4425/4--b/8 x«+4+0 
9x/2+27/2—9b/4 9x/4+33/4+9b/4 9x/4+33/4 
7x/4+25/4+7b/4 9x/4+33/4+9b/4 11%/24+7/44+11b/2 5x/2+53/84+5b/2 5x/2+7/44+5b/2 | =0. 
5x/2+53/8+5b/2 11x%/2+5/2—11b/4 11x/8+4+11b/8 
3x/2+9/2+3b/2 5x/2+7/44+5b/2 11%/8+4+11b/8  5x/2—11/4+5b/2 


Fic. 1. Complete set of canonical structures for 1,2-dihydro- 
naphthalene. 


3x/2+9/2+3b/2 


interactions involving the single bonds between 
adjacent atoms are neglected, on the reasonable 
assumption that they are effectively the same for 
all the structures considered.) Hence the 1,2 
compound is calculated to be more stable than 
the 1,4 compound by —0.2007a or, for a= —1.5 
v.e., 0.301 v.e. This is in rough agreement with 
the experimental value 5.0 kcal. or 0.22 v.e. given 
above. The comparison indicates that 6 and §’ 
differ from a by a few percent. If we assume, as is 
reasonable because of the relative interatomic 
distances, that B lies midway between 8’ and a, 
and that a safe lower limit for the extra stability 
of the 1,2 compound is 4.6 kcal. or —0.14a, we 
write as a second equation connecting } and 2 


2=4.2455—2b. (2) 
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Fic. 2. 


The simultaneous solution of Eqs. (1) and (2) is 
shown in Fig. 2, in which the two curves represent 
the two relations between the variables. Their 
intersection corresponds to B=0.935a and £’ 
=0.87a, which are thus minimum values for 
these exchange integrals. 

It is probable that the errors introduced by the 
neglect of all except unexcited and first-excited 
structures and by equating 8 and #’ to a@ cancel 
each other to a considerable extent, especially for 
more complicated molecules than the dihydro- 


x+5/2 x/2+7/2 
x/24+7/2 17x/8+103/16 
35x/8+361/16 19x/2+47 
19x/16-+47 


The pertinent root of this is x = —4.599, corre- 
sponding to the energy W1.2=Q—xa—3a/2=Q 
+3.099a, of which Q+2.976a results from reso- 
nance within the naphthalene nucleus (calculated 
by neglecting all structures other than the three 
unexcited and the sixteen naphthalene first- 
excited structures), and the difference, 0.123a, is 
the energy of conjugation of the naphthalene 
nucleus and the double bond. The total energy of 
1,4-dihydroanthracene, in which this conjugation 
does not occur, is Wi,4=Q+2.976a, and that of 
9,10-dihydroanthracene, involving two unconju- 
gated benzene rings, is Wo,10= Q+3.190a (Table 
II). Table II could be easily extended to include 


35x/8+361/16 
19x/2+47 
557x/8+4627/16 427x/32+2599/32 
35x/16+403/32 427x/32+2599/32 119x/8+595/16 


naphthalenes. These approximations are made 
throughout the rest of this paper. 

Calculated resonance energies for 1,2 —CioH 10, 
1,4—Cy Ho, 1,5—CyoHi0, and 9,10—CiH 0, the 
last two involving a chain of four and two chains 
of two conjugated double bonds (discussed in a 
following section), are given in Table I. No data 
are available regarding the two least stable 
dihydronaphthalenes. 


TABLE I. Calculated resonance energies 
for dihydronaphthalenes. 


W Ww- Wi,2 Ww- Wi,2 
1,2—CioHi0 0 0 


1,4—CyoHi0 +1.6055a —0.2007a 0.301 v.e. 
1,5—CioHi0 +1.225a —0.58la 0.872 
9,10—CioHi0 +0.964a —0.842a 1.263 


ADDITION PRopUCTS OF ANTHRACENE 


The secular equation for 1,2-dihydroanthracene 
resulting from the consideration of (I) the 
symmetrical unexcited structure, (II) the two 
unsymmetrical unexcited structures, (III) the 
sixteen first-excited structures not involving 
conjugation of the double bond and the naphtha- 
lene nucleus, and (IV) the six first-excited 
structures involving this conjugation, is 


19x/16-+7 
35x/16-+403/32 


additional less stable and as yet unprepared 
dihydroanthracenes and higher addition prod- 
ucts. No thermal data are available for com- 
parison with the calculated energy values of 
Table II. 


TABLE II. Calculated resonance energies 
for dihydroanthracenes. 


W— Wo,10 W—Wa,10 


9,10—CisHis 0 0 
1,2—CuHie —0.091a 0.137 v.e. 
1,4—CuHw —0.214a 0.321 v.e. 


The resonance energies are essentially independent of 
the nature of the atoms attached to the nucleus; hydrogen 
is mentioned for convenience. 
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It is, however, reasonable to assume that, in 
the case of closely similar reactions, that reaction 
leading to the most stable product will occur 
most easily. This hypothesis leads to a complete 
explanation of the observed courses of hydro- 
genation of benzene, naphthalene, and anthra- 
cene, which are represented by the following 
schemes, showing only the reactions for which 
data are available.® 


CsHo+H21,2 —CeHs+He 
21,2,3,4-— C,H 1 +H. 


CipHs+H2-1,2 —CyoH 19 + H2-1,2,3,4 —CroH iz. 


CuHi0+H2 9,10 —CuHie+He 


—1,2,3,4,5,6,7,8 CyH 


With this hypothesis, the calculated stability 
of the 9,10 addition compound of anthracene 
provides an explanation of the ease of attack of 
the 9,10 positions for addition. A similar calcu- 
lation for phenanthrene shows that for this 
molecule too the 9,10 positions should be most 
reactive. 


CONJUGATED DouBLE-BoND SYSTEMS 


In Table III are given calculated values of the 
resonance energy for straight conjugated chains 
containing two, three, and four double bonds and 
a branched chain with three double bonds (the 
list could be easily extended). In each case the 
secular equation was solved as a quadratic, all 
first-excited structures being given the same 


TABLE III. Resonance energies in systems of double bonds. 


Number of 
first-excited 
System structures Resonance energy 
1 —0.232a=0.348 v.e. 
ses 3 —0.482a=0.723 v.e. 
=—-=—-=—= 6 —0.725a=1.09 v.e. 
» 2 —0.438a =0.657 v.e. 


5 W. Hiickel, Theoretische Grundlagen der orgonischen 
Chemie, 1, p. 359, Akad. Verlag. M. B. H., Leipzig, 1931. 


coefficient. No reliable thermochemical data are 
available for comparison with these values. It 
may be mentioned that our hypothesis provides 
no information regarding addition in various 
positions (1,2 vs. 1,4 in dienes), which must 
instead be discussed on the basis of detailed 
mechanisms. 


PHENYLETHYLENE, STILBENE, ETC. 


The secular equation for phenylethylene, in 
which a phenyl group is conjugated with a 
double bond, is the same as for 1,2-dihydro- 
naphthalene, and leads to the same value for the 
resonance energy. 


O---0 0---0 0---0 0---0 


Q---O @---O O---@ O---@ 
O-=-O O-=-O O-=-Y O-=-O 


Q=3-O O=2-O 
Q=2-O O=5-O 
O-s20 
O-<=0 


Fic. 3. Unexcited and first-excited canonical structures of 
sym.-diphenylethylene. 


For stilbene, sym.-diphenylethylene, we con- 
sider the thirty-seven structures shown in Fig. 3, 
grouping together the four unexcited structures, 
the twelve first-excited structures not involving 
conjugation, the twelve in which there is conju- 
gation of the double bond and a phenyl group, 
and the nine in which there is conjugation of the 
two phenyl groups by way of the double bond. 
This leads to the biquadratic secular equation 
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25x/4+30 15x+81 35x/4+453/8 49x/16+399/16 
15x+81 81x/2+207 87x/4+1185/8  63x/8+519/8 
35x/4+453/8 87x/4+1185/8 269x/8+1107/8 77x/4+753/8 
49x/16+399/16 63x/8+519/8  77x/4+753/8 121x/4+759/8 
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which on solution as a biquadratic leads to 
W=(Q+5.5889a. The energy the molecule would 
have without resonance, found by equating the 
upper left-hand quadratic part of the secular 
determinant to zero, is W’=Q+5.1898a, corre- 
sponding to an extra resonance energy of con- 


The secular equation for isostilbene, in which 
there is no conjugation of the two phenyl groups 
with each other, is just the cubic equation 
obtained from the stilbene determinant by 
omitting the fourth row and column. 

The secular equation for triphenylethylene, 
corresponding to 8, 36, 36, and 36 structures of 


scaendaniaabiiaias the four kinds described for stilbene, is =: 

125x/8+1675/16  225x/4+3285/8 525x/16+4395/16 245x/16-+4921/32 
225x/4+3285/8  1755x/8+25,569/16  1935x/16+17,019/16 231x/4+2379/4 
525x/16+4395/16 1935x/16+17,019/16 2385x/16+31,653/32 3423x/32+12,465/16 
245x/164+4921/32 231x/44+2379/4 3423x/32+12,465/16 2959x/16+3985/4 


That for tetraphenylethylene, with 16, 96, 196, and 144 structures, respectively, of these four 


kinds, is 

625x/16+5375/16  375x/2+1725 875x/8+17,975/16 1225x/16+7315/8 
375x/24+1725 3825x/4+35,685/4  4275x/8+91,815/16 3045x/8+18,699/4 
875x/8+17,975/16 4275x/8+91,815/16 9175x/16+41,245/8 12,845x/32+273,911/64 
1225x/16+7315/8  3045x/8+18,699/4  12,845x/32+273,911/64 74,761x/64+557,465/64 


= = 


Ky 


The energy values obtained on solution of these 
equations are given in Table IV. 


CONJUGATED BENZENE RINGS 

For biphenyl we consider four unexcited 
structures, twelve first-excited structures not 
involving conjugation, and nine involving con- 
Jugation, as shown in Fig. 4. These lead to the 
secular equation 
25x/4+215/8 15x+147/2  49x/8+301/8 
1Sx+147/2  81x/24-747/4 63x/4+393/4 |=0, 
49x/8-+301/8 63x/4+393/4 121x/4+319/4 


whose solution is given in Table V. 

For ortho, meta, and para diphenylbenzene we 
consider eight unexcited structures of the type 
indicated by the letter A in Fig. 5, thirty-six 
first-excited structures of type B, and thirty-six 
of type C, in which there is conjugation with the 
central ring. This includes all first-excited struc- 
tures for m-diphenylbenzene. For ortho and para 
diphenylbenzene, however, there are in addition 
nine first-excited structures of type D, in which 
the two end rings are conjugated together. The 
secular equation for the ortho compound is 
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The secular equation for the meta compound is obtained from this by omitting the fourth row 


and column. That for the para compound is: 


125x/8+775/8  225x/4+765/2 245 
225x/4+765/2  1755x/8+11,907/8 231 


For 1,3,5-triphenylbenzene, with 16 unexcited, 96 unconjugated first-excited, and 108 conjugated 


first-excited structures, the secular equation is 


625x/16+10,125/32 375x/2+6525/4 


375x/2+6525/4 3825x/4+67,545/8 


125x/8+775/8 225x/4+765/2 245x/8+3941/16 
225x/44765/2 1755x/8+11,907/8 231x/24+1917/2 
245x/84+3941/16  231x/2+1917/2  5771x/32+7813/8 
245x/64+1169/32 231x/164+4527/32 539x/16+3367/16 


«/8+3941/16 49x/16+497/16 
399x/32+981/8 
245x/8+3941/16 231x/24+1917/2 2959x/16+16,151/16 77x/24+1891/8 
49x/164+497/16  399x/324+981/8  77x/2+1891/8 


x/24+1917/2 


3675x/32+9135/8 


9135x/16+23,481/4 


2452/64 + 1169/32 
231x/164+4527/32 
539x/16+3367/16 
121x/4+110 


121x/4+110 


3675x/32+9135/8  9135x/16+23,481/4 50,265x/64-+799,227/128 


é 7? 
06 00 O 
& 
a “a 
OO DO 
TABLE IV. Calculated resonance energies for phenyl- 
7 48 
Molecule w w’ WwW’ -W 
Ox) =) Phenylethylene 3.3062a Q+3.105Sa —0.2007a =0.301 v# 
20 2 a2 23 Sym. diphenyl- 
ethylene 5.5892 Q+5.190a —0.389a 
Unsym. diphenyl- 
Ox) ethylene O+ 5.545@ Q45.190a -0.355a =0.533 | 
Triphenylethylene Q+ 7.804@ Q+7.285a —0.519a =0.779¥% | 
Tetraphenyl- 
ethylene Q+10.080a Q0+9.362a —0.719a =1.08 Vo 


Fic. 4. Unexcited and first-excited structures of biphenyl. 
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The energy values found on solution of these 
equations and of the non-conjugated parts of 
them, together with their differences, the energy 
of conjugation, are given in Table V. 


TABLE V. Calculated resonance energies for conjugated 
benzene rings. 


Molecule Ww Ww’ WwW’ -W 
Biphenyl 0+4.867a Q0+4.690a —0.177a =0.266 v.e. 
m-Diphenylbenzene Q+7.086a QO+6.785a —0.301la =0.452 v.e 
p-Diphenylbenzene Q4+7.094a 0+6.785a —0.309a =0.464 v.e 
o-Diphenylbenzene Q+7.103a Q+6.785a —0.318a =0.477 v.e 
1,3,5-Triphenylbenzene Q+9.273a Q+8.862a  -—0.41la =0.617 v.e. 


COMPARISON WITH EXPERIMENT 


The meager existant results of experiments on 
heats of combustion, given in the previous paper 
of this series,® are in rough agreement with the 
theoretical calculations, the discrepancies in most 
cases lying within the experimental error. The 
calculated energy of conjugation of a benzene 
ring and a double bond, 0.30 v.e., is in excellent 
agreement with the extra resonance energy 
0.29 v.e. for phenylethylene. For stilbene and 
isostilbene’? the empirical values 0.67 and 0.45 
v.e. agree less well with the calculated values 
0.58 and 0.53 v.e., respectively. For biphenyl the 
agreement is good, 0.35 v.e. empirical, 0.27 v.e. 
calculated; for 1,3,5-triphenylbenzene the differ- 
ence of 0.5 v.e. between the empirical value 1.1 
v.e. and the calculated value 0.62 v.e., while 
large, may be due to error in the reported heat of 
combustion or in part in the estimated heat of 
sublimation. 


CONCLUSIONS REGARDING CONJUGATION 


The calculations provide us with sufficient 
information to permit the formulation with con- 
siderable confidence of a set of rules regarding 
conjugation. 

A phenyl group is 20 or 30 percent less effective 
In conjugation than a double bond, and a 
naphthyl group is less effective than a phenyl 
group. Thus the resonance energy of bipheny] is 
only 77 percent of that for two conjugated 


*L. Pauling and J. Sherman, J. Chem. Phys. 1, 606 
(1933), 

"The value for this compound, omitted from the 
Previous paper, is calculated from the heat of combustion 
by the methods described there. 
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double bonds, and the ratios of energy of con- 
jugation of a double bond with another double 
bond, a phenyl group, and a naphthyl group are 
1 : 0.86 : 0.53. 

If a conjugated system can be resolved into 
two overlapping parts with a common double 
bond or benzene ring, and with no conjugation 
from one part to the other by way of the double 
bond or benzene ring, the total conjugation 
energy is 10 or 15 percent less than the sum of the 
conjugation energies for the two parts. Examples 
are isostilbene, with conjucation energy 1.77 
times that for phenylethylene, metadiphenyl- 
benzene, 1.70 times the biphenyl energy, and 
branched triene, 1.89 times the diene energy. 
The conjugation energy of 1,3,5-triphenyl- 
benzene, in which the central benzene ring forms 
a part of three biphenyl systems, is 23 percent 
less than three times the biphenyl conjugation 
energy. 

In an overlapping system in which there is 
conjugation from one part to another by way of a 
double bond this extra conjugation may counter- 
act the effect of overlapping, so that the total 
conjugation energy is approximately equal to the 
sum for the parts. Thus the conjugation energy of 
stilbene is 1.94 times that for phenylethylene, 
that for unbranched triene 2.08 times that for 
diene, and for unbranched tetraene 3.12 times 
that for diene. A phenyl group is not so effective 
as a link in a conjugated chain, however; no 
doubt for the reason that, because of resonance 
between the two Kekulé structures, the average 
path between ortho or para positions is three 
bonds long. In consequence the conjugation 
energy for p-diphenylbenzene is only three 
percent greater than for the meta compound, and 
that for the ortho compound (in which the 
bridge is more effective) only six percent greater. 

Conjugation of phenyl groups and double 
bends is of quite appreciable effect in stabilizing 
molecules, as is, of course, evident from the early 
recognition of the phenomenon by chemists. The 
extra resonance energy of tetraphenylethylene, 
amounting to 1.08 v.e., is more than enough to 
overcome the unsaturation of the double bond, 
amounting to 0.74 v.e. (the difference in energy 
of two carbon-carbon single bonds and a double 
bond), in consequence of which the properties of 
the substance differ radically from those of an 
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unsaturated compound such as ethylene.® Using 
the bond-energy values given in earlier papers,® 
we calculate for AH for the reactions Co(CsHs). 
+X, the values —0.30 v.e., —0.42 
v.e., and +0.24 v.e. for X2=He, Cle, and Bro, 
respectively. Assuming that AF and AH do not 
differ greatly, these values suggest that tetra- 
phenylethylene will add hydrogen and chlorine 
but not bromine.’ This is in agreement with 
experiment; tetraphenylethylene has been hydro- 


8 We are indebted to Professor G. E. K. Branch and to 
Dr. G. W. Wheland for pointing this out to us. 

9L. Pauling and J. Sherman, J. Chem. Phys. 1, 606 
(1933); L. Pauling, J. Am. Chem. Soc. 54, 3570 (1932). 

10 The values of AH for addition to a non-conjugated 
double bond are 1.08 v.e. less than these, so that such 
addition reactions are all exothermic. 


genated catalytically to tetraphenylethane," and 
has been shown to add chlorine directly,” but it 
does not add bromine, and attempts to prepare 
tetraphenyldibromethane have led to tetra- 
phenylethylene itself." 


1 W. H. Zartman and H. Adkins, J. Am. Chem. Soc. 54, 
1668 (1932). 

2H, Finkelstein, Ber. Dtsch. Chem, Ges. 43, 1533 
(1910); J. F. Norris, R. Thomas and B. M. Brown, ibid. 
43, 2940 (1910). See also J. Schmidlin and R. v. Escher, 
ibid. 43, 1153 (1910). The resultant dichlortetraphenyl- 
ethane decomposes easily, mainly with the liberation of 
HCI and the formation of tetraphenylethylene with Cl 
substituted in the para position in one of the phenyl 
groups. 

13H, Biltz, Ann. d. Chemie 296, 231 (1897). H. Finkel- 
stein, reference 12. 
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The Virial and Molecular Structure 
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The virial theorem can be applied to a molecule, assum- 
ing that external forces are applied to keep the nuclei fixed. 
The mean electronic kinetic energy is then the negative of 
the mean total internal energy, plus a term arising from 
the virial of the external force. Since the external force is 
derivable from the knowledge of the total internal energy 
as function of nuclear position, this theorem gives a means 
of finding kinetic and potential energy separately for all 
configurations of the nuclei, as soon as the total energy is 


known, from experiment or theory. By using simple forms 
of curve for total energy, for the diatomic molecule, the 
kinetic and potential energy are derived and discussed. 
The curves give direct indication of the formation of 
the covalent bond, in cases where this occurs, by describing 
the removal of charge from the atoms to the region between 
the atoms, with resulting decrease of kinetic energy, 
increase of potential energy. 


HE study of molecular structure has been 
too much based on particular models, 
rather than on general principles. One such 
general principle which has not been used and yet 
which seems capable of shedding considerable 
light on molecular problems, is the virial theorem. 
It is known that this theorem holds in quantum 
as in classical mechanics (see the Appendix for a 
simple proof) and its application to the theorem 
that the kinetic energy of an atom is the negative 
of the total energy is familiar. We shall use it, 
however, in a rather different way, applying it in 
the case where there are impressed forces. 


THE VIRIAL THEOREM WITH IMPRESSED FORCES 


The familiar statement of the virial theorem is 


where the term on the right is a summation over 
all coordinates of the system, and F, represents 
the force acting on the x coordinate, etc. We 
apply the theorem to a molecule consisting of 
electrons and nuclei attracting and repelling 
according to the inverse square law. Then, by a 
familiar argument, the virial (xF,+---) of the 
inverse square forces equals the potential energy 
of those forces, so that if the system is left to 
itself, the mean kinetic energy is minus half of 
the mean potential energy. But we shall not leave 
the system to itself. We shall rather assume that 
external forces of constraint act on the nuclei, 


just sufficient to keep the nuclei at rest, instead of 
allowing them to vibrate as they would if they 
were free. In applying the virial theorem, we 
shall now have to compute merely the electronic 
kinetic energy, since there is no nuclear kinetic 
energy of vibration. The part of the virial 
coming from the inverse square forces will be, as 
usual, equal to the potential energy of these 
forces. But, in addition, there will be a contri- 
bution to the virial coming from the external 
forces of constraint. Thus the kinetic energy will, 
on the average, be different from minus a half of 
the mean potential energy. It will be equal to 
this quantity only in two cases: at infinite 
separation of all atoms from each other and at 
the equilibrium configuration, in which the atoms 
could exist without vibration at the absolute 
zero. For, in both these cases, the atoms would 
remain in equilibrium without external forces of 
constraint, these forces would have no virial and 
only the virial of the inverse square forces need be 
considered. This is true only to the extent to 
which we can neglect zero point vibration; it 
seems that the arguments of the present paper 
are all accurate only to that approximation, 
since it is not exactly legitimate to replace the 
coupling between electronic and nuclear motion 
by an external force acting on the nuclei. 

To fix our attention on a definite problem, we 
shall consider a diatomic molecule, though the 
results we shall obtain can easily be generalized 
to any molecule or solid. Let the impressed force 
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on nucleus 1 be F,, on nucleus 2 F2, where Fs 
=-—F,, to keep the molecule as a whole in 
equilibrium. Then 


The vector r2—17; is simply the vector from the 
first nucleus to the second. Since F? is parallel to 
this vector, with a diatomic molecule, the result 
is rF, where r is the internuclear distance, F the 
force to hold the nuclei in equilibrium. If the 
nuclei are attracting each other, this force is 
outward and the product is positive, while if 
they are repelling it is negative. Thus for a 
diatomic molecule in an attractive region, the 
mean kinetic energy is less than minus half the 
mean potential energy of the electrostatic forces, 
while in a repulsive region it is greater. 


THE KINETIC ENERGY OF A DIATOMIC MOLECULE 


J. C. SLATER 


lation of E will also permit the calculation of 
T and V directly from the model, without use of 
the virial theorem. These virial equations, then, 
will provide a valuable and independent check of 
the correctness of the wave functions of the 
model. This check in general is much more 
sensitive than that obtained from an approxi- 
mately correct energy. For it is well known that 
the energy is stationary with respect to a slight 
variation of the wave function, so that almost 
the correct energy will be found from a poor wave 
function, whereas the relation between kinetic 
and potential energy changes rapidly for even 
small changes in wave function. 


KINETIC AND POTENTIAL ENERGY IN TYPICAL 
CASES 


First we consider an attractive state of a 
molecule. At distances greater than the equi- 


7 ; librium distance, T is less than —£, going 
Again, confining ourselves to the simple case of through a minimum in general, rising again, and [| 
the diatomic molecule, we can easily see that the equalling —E at equilibrium. It is rising rapidly | 


results of the last section allow an unambiguous with decreasing distance at equilibrium. It does | 
determination of the mean electronic kinetic pot continue to rise indefinitely, however. For as si 
energy of the molecule as a function of inter- the nuclei of the two atoms are pushed into | 4 
nuclear distance, as soon as the curve of total coincidence, the atoms coalesce into a single [| ‘i 


atom whose atomic number is the sum of the 


electronic energy as function of distance (the | 
atomic numbers of the two atoms and, of ot 


i familiar ‘‘potential energy curve” of the mole- 


cule) is known. Let this electronic energy, the course, this combined atom has a finite kinetic J 
sum of the electronic kinetic and potential energy. This is in agreement with the result of i Si 
energy (T+V) be called E. Then the force the virial theorem. For, at very small distances, j = 
exerted by the nucleus 2 is —dE/dr and the force the total energy E will be the sum of two terms: | of 
F which equilibrates this is dE/dr. Hence we the repulsion Z,Z2e?/r of the two nuclei of in 
f" : have, dropping the average signs for convenience, atomic numbers Z; and Zs, and the energy E’ of [ 2 
tee the electrons of the combined atom. Thus for 
very small r’s, we have 
as 
with the equation T= 
—r(d/dr)(E' +Z,Z2e?/r) = 
T=—-—E-r(dE/dr), V=2E+r(dE/dr). sta 
V=2(E' +-1(d/dr)(E’ exc 
These important equations determine the mean =2E'+Z,Z002/'. ay 
kinetic and potential energies as functions of 7, ene 
one might almost say, experimentally, directly That is, the kinetic energy of the electrons is the but 
from the curves for E as function of r which can negative of the total energy of the combined gOir 
be found from band spectra. The theory is so atom, a finite amount, while the potential energy 
simple and direct that one can accept the results is the sum of twice the total energy of the 
without question, remembering only the limi- combined atom and the nuclear repulsion, al! 
tation of accuracy mentioned above. Any model just as we should expect. For a repulsive state the 
of molecular structure which permits the calcu- situation is similar except that the preliminary 
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decrease of kinetic energy at large r is not 
present. 

These results become more definite if we take 
some of the well-known approximate formulas 
for interatomic energy and compute kinetic and 
potential energy from them. We should state at 
the outset that no one of these formulas has the 
correct behavior at r=0. Hence we can expect a 
reasonable behavior only in the neighborhood of 
the equilibrium distance and for larger distances. 
First we take the power series formula used by 
Born and others for crystal structure. This gives 


E=E,—a/r"+b/r", 
T = —E,,—(n—1)a/r"+(m—1)b/r", 


where E,, is the energy at infinite separation. 
Assuming that m>n, and that both m and n are 
greater than unity, we see that the kinetic 
energy first decreases with decreasing 7, as we 
have stated, then goes through a minimum 
similar to that of the total energy, rising again 
and becoming infinite as r=0. The minimum of 
kinetic energy is easily shown to be at a larger r 
than the equilibrium distance, by a factor 
{(m—1)/(n—1)}“"-™), which, for example, is 
about 1.2 in the typical case m=9, n=5. 
Similarly, V rises for large r, goes through a 
maximum at a slightly larger 7 than the minimum 
of the kinetic energy, then decreases without 
limit. This, of course, neglects the fact which we 
have pointed out, that at still smaller distances 
the potential energy goes through a minimum 
and again increases, becoming positively infinite 
as r approaches zero. In Fig. 1: we show curves for 
E, T and V, for the case m=9, n=5. The case 
n=1, where the attractive forces are electro- 
static, as in an ionic crystal or molecule, is 
exceptional. For in that case the second term in 
T, involving the factor n—1, is zero. The kinetic 
energy then does not go through its initial dip 
but begins to increase immediately as 7 decreases, 
going up with the same power of 1/r as the 
repulsive potential. The term —a/r, appearing in 
E, comes, in this case, entirely from the potential 
energy. 

_ A purely repulsive state can be approximated, 
in terms of inverse powers, by a single term 


Fic. 1. Kinetic, potential and total energy of diatomic 
molecule as function of 7. Energy is assumed to be of 
form —a/r°+b/r°, 


b/r™, obtained in the formulas above by setting 
a=0. Here we see that the kinetic energy 
increases, and the potential energy decreases, 
with decreasing 7, according to the same inverse 
power as the total energy. In this case again, the 
behavior for very small 7 is not correctly given. 
One interesting feature, which is met also in the 
attractive case, is that if m is large, the kinetic 
and potential energies separately are much larger 
than the total energy, which appears as the small 
difference of two large quantities. 

The features of which we have spoken are not 
peculiar to the particular energy formula which 
we have used. If the energy is approximated by a 
sum of two exponential terms, as in the Morse 
curve, or in a repulsive case by a single ex- 
ponential, the behavior in the neighborhood of 
the minimum and for large r is essentially 
similar to what we have described. The formulas 
are easily derived and it is hardly necessary to 
discuss them in detail. 


INTERPRETATION OF KINETIC AND POTENTIAL 
ENERGY CURVES 


The curves for kinetic and potential energy 
which we have discussed in the preceding section 
and plotted in Fig. 1, seem at first sight to be 
unexpected. We consider first the attractive case, 
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taking homopolar valence, for which m would be 
about 5 (or for which we would use a Morse 
curve, which gives results like Fig. 1). For we see 
that the decrease of energy which is responsible 
for binding seems to come from a decrease of 
kinetic energy, rather than of potential energy. 
The theory of homopolar binding suggests, on 
the other hand, that the mechanism of this 
binding is one in which charge concentrates 
between the nuclei and is attracted electro- 
statically by both nuclei, this attraction pro- 
ducing the binding. Thus we should expect that 
the binding would result from a decrease in 
potential energy, rather than kinetic. The answer 
to this paradox is not far to seek, however. We 
have just said that the first step in homopolar 
binding is the concentration of charge in the 
region between the nuclei. Now it requires a 
relatively small amount of energy (as we can see 
from the minimum principle for energy) to pull 
a certain amount of charge from the inner part of 
the atom, where it is ordinarily found, out to the 
periphery, as for example the region between the 
atoms. It can very well be that the decrease of 
potential energy made possible by this change, 
coming from the attraction of the electron of one 
atom by the nucleus of the other, can decrease 
the total energy enough to make this the stable 
configuration and produce binding. But when the 
charge is pulled out to the edge of the atom, its 
potential energy increases greatly, as it is 
separated from the nucleus, and its kinetic energy 
correspondingly decreases. Thus the rise of 
potential energy and the fall of kinetic energy, 
noticed in Fig. 1, are produced. As the atoms are 
pushed further together, however, a point is 
reached where the charge to be shared is fully 
removed to the region between the atoms. 
Further decrease of r then produces no further 
removal of charge to the periphery, no further 
decrease of kinetic and increase of potential 
energy. Instead, the other, more fundamental 
influences come into play. The kinetic energy 
begins to rise as the electrons are forced into 
smaller space. This may be likened to the 
adiabatic rise of energy of an electron gas 
compressed in a box. And the potential energy 
falls, as the negative charge concentrated be- 
tween the atoms and the positive charge left 
over on the atoms come closer together. The 
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attractive force of the electrical attractions 
balances the repulsion of the electronic motion at 
equilibrium, the positive slope of the potential 
energy curve indicating large attraction, and the 
negative slope of the kinetic energy large 
repulsion. The curves show, then, that in the 
neighborhood of the equilibrium position, and 
even more at smaller distances, the behavior of 
the valence bond is rather simple and quite as we 
should expect. The peculiar behavior of potential 
and kinetic energies at larger r can be definitely 
associated with the formation of the shared 
electron pair which is responsible for the bond. 
This interpretation can be made rather unam- 
biguously without use of models; for it is hard to 
see any meaning which can be given to an 
increase in electronic potential energy, decrease 
in kinetic, with much smaller change in total 
energy, other than a shift of electrons toward the 
outside of the atom. 

Our identification of the minima and maxima, 
respectively, of the kinetic and potential energy 
for large r with the formation of a shared electron 
pair becomes much more reasonable when we 
consider the other cases of energy curves. We 
have already seen that these minima and maxima 
do not occur in the case where the attractive 
force is electrostatic, as in ionic molecules (7 = 1), 
or in repulsive molecular states. In both these 
cases, the usual theory is that no electron pair is 
formed. The only action is the electronic re- 
pulsion when the atoms come closer together and, 
in the ionic case, the electrostatic attraction. 
This fits in perfectly with the behavior of the 
kinetic energy which rises rapidly, as_ the 
repulsive part of the energy does, at small 
distances. 

The considerations which we have suggested 
indicate the importance of the virial in discussing 
molecular structure. There has been no attempt 
in this paper to deal exhaustively with its 
importance but only to point it out as a method 
which should be used in future work. And it 
should be recalled that it is as useful in treating 
more complicated molecules and solids, as with 
diatomic molecules. With a solid, for instance, if 
we subject it to a hydrostatic pressure p, the 
theorem becomes T= —E+3u, if the internal 
forces are electrostatic, as is well known. Thus, 
here, the kinetic energy is greater than —E for 
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positive pressures, less for negative pressuresand, atomic forces can be drawn here as for molecules. 
in general, shows just the sort of behavior as Such conclusions were drawn before the wave 
function of volume that we have discussed for mechanics by Schottky,! but his paper seems to 
diatomic molecules as a function of r. The same have had but small influence in the theory of 
conclusions regarding the nature of the inter- solids. 


APPENDIX 
Let Schrédinger’s equation be 
h? 
—+(V-Ey=0. 


é 822m; Ox ;? 


We differentiate with respect to x;, and multiply by x;y, obtaining 
h? dy 


xy + 
Ox 


oy 
Ox; 0. 


Xj 


Now by Schrédinger’s equation for y, the last term is equal to 


We substitute this term in and sum over x;, obtaining 


i 822m; i Ox Ox; 


We now integrate this equation over the coordinate space. In doing this, we integrate the first term 
by parts. We do this by means of the identity 


Ox Ox ;? Ox; Ox; Ox, 


which can be proved, after a little manipulation, by performing the differentiations indicated on the 
right. In integrating the right side of this over the coordinates, the derivative integrates to zero if we 
assume the system to be a closed one, for then we take the limits of integration at infinity, where 
¥=0. Hence we are left with only 


h? _ oy —dV\_ 
( Ww dr. 


8 3m; OX; 


The left side is recognized as the mean kinetic energy, the right side as —} times the mean of 
Lix;(—8V/dx,), which is just the virial, when we remember that F;= —0V/dx;. 


‘Schottky, Phys. Zeits. 21, 232 (1920). 
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The Free Energy, Entropy and Heat Capacity of Bromine and of Hydrogen Bromide 
from Spectroscopic Data 


The free energy, entropy and heat capacity of gaseous 
bromine and of hydrogen bromide have been computed for 
the temperature range 200—-1600°K. The equilibrium con- 
stants for the reaction Br.=2 Br, obtained from these 
numbers and the spectroscopic heat of dissociation, differ 
from Bodenstein’s values by more than his apparent 
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experimental error. Calculations for the dissociation of 
hydrogen bromide give results which agree with the e.m.f. 
measurements of Bodenstein and Geiger, but are in definite 
disagreement with the experiments of von Falckenstein 
when the latter are corrected for the dissociation of bro- 
mine. 


N A recent paper,! Brown, who had effected 

the rotational analysis of the bromine spec- 
trum,” reported that a calculation of the equi- 
librium constant for the dissociation of diatomic 
into monatomic bromine ‘‘by the Gibson and 
Heitler equation”’ leads to values of log K which 
were uniformly about 0.15 less than the experi- 
mental results of Bodenstein.’ This discrepancy 
was so surprising in view of the excellent agree- 
ment found by Gibson and Heitler‘ in the case of 
the iodine. dissociation that it seemed worth 
while to recompute the equilibrium. Brown gave 
no details of his calculation, but presumably 
treated the rotational-vibrational levels of the 
bromine molecule as those of a rigid rotator 
combined with a harmonic oscillator; it did not 
seem likely that the error introduced by this 
simplification could account for the serious 
discrepancy (0.7 cal./deg. in R In K) which he 
reported, but we have nevertheless computed the 
free energy, entropy and heat capacity for the 
temperature range 200—1600°K, taking account 
of the interaction and stretching terms in the 
formula for the energy levels. Since the dis- 
sociation of hydrogen bromide into hydrogen 
and bromine has been measured by von Falcken- 
stein,» we have also carried out similar calcu- 
lations for hydrogen bromide. 


1W. G. Brown, J. Am. Chem. Soc. 54, 2394 (1932). 
2 W. G. Brown, Phys. Rev. 39, 777 (1932). 

3 Bodenstein, Zeits. f. Elektrochemie 22, 327 (1916). 
4 Gibson and Heitler, Zeits. f. Physik 49, 465 (1928). 


5yvon Falckenstein, Zeits. f. physik. Chemie 68, 270 


(1909). 
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Brown? gives for Br??—Br®!, €,=323.2(v+}) 
—1.07(v+3)*, B,=0.08077—0.000275v, 
= —2.03X10-°. Owing to the existence of the 
two isotopes of bromine in almost equal pro- 
portions (0.505 : 0.495), the three diatomic va- 
rieties will exist at high temperatures in the 
proportions 0.255 of 79-79, 0.500 of 79-81, 0.245 
of 81-81. As far as the calculation of thermo- 
dynamic reaction quantities is concerned, it has 
been shown‘ that it is permissible to calculate the 
free energy and entropy as if the substance were 
non-isotopic, provided that this procedure is 
adopted for each of the reactants and products, 
and that the spectroscopic constants used in the 
computation are averages, the result of properly 
weighting the constants for each of the isotopic 
varieties. Calculation shows that such effective 
constants for bromine are (within the limits of 
experimental error) just those for 79-81 and the 
calculation can thus be carried out with these 
constants as if for a non-isotopic bromine of 
molecular weight 159.83; it is necessary, of 
course, to add the usual symmetry correction for 
a homonuclear diatomic molecule, —R In 2, in 
spite of the fact that all rotational levels are 
present in 79-81 with full statistical weight. 
The free energy, entropy and heat capacity 
entered in Table I’ were obtained from Gordon 


6 Giauque and Overstreet, J. Am. Chem. Soc. 54, 1731 
(1932); see also Gordon and Barnes, J. Phys. Chem. 36, 
2292 (1932) and reference 1. 

7In the calculations, the values of the universal con- 
stants are those given in Int. Crit. Tab., Vol. I. 
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THERMODYNAMIC PROPERTIES OF BRz AND HBR 


and Barnes,’ Eqs. (6), (9) and (11). For tempera- 
tures less than 700°K, the quantities (In Q, 
+eohc/kT), S,/R, C,/R, 5, v*, 51, $2, C1, C2 were 
obtained by interpolation in their Tables I—XI; 
for temperatures $700°K, w,hc/kT lies beyond 
the range of the tables and therefore the calcu- 
lation was carried out as shown in their Example 
4, the integrations starting from v= 11, as in that 
example. There is an additional small correction 
that contributes to the heat capacity at the 
highest temperatures but not to the entropy or 
free energy; this is due to the excited electronic 
> tes? 14,734 and 15,891 cm! above the normal 
'Y s ate of the molecule. The first of these excited 
states is unclassified and no correction can be 
applied as yet but the second is presumably *II 
with e,’=169.7(v' +3) —1.913(v’'+3)? and B’ 
=(0.0592. To a good enough approximation, the 
resulting contribution to the heat capacity is 
given by 


kT 


By De-evhelk T 


where A =15,891. The numerical value of this 
term is only 0.004 for 1600° and 0.002 for 1500° 
and is included in the numbers entered in Table I. 

In the case of hydrogen bromide, the con- 
stants!” are e, = 2647(v+3) —44(v+3)*, B, =8.342 
—0.225v, D=—0.000,325, and the quantities 
necessary for the calculation can be obtained for 
the whole temperature range from the tables. 
The effect of the H? isotope of hydrogen on the 
effective spectroscopic constants has been ig- 
nored, since this is vanishingly small!! except for 
artificially enriched samples. 

None of the values of —(F°—E,°)/T and S° in 
Table I includes the nuclear spin contribution. 
For bromine (since j,= 3/2) this would amount to R 
In 16 and must be added to the entries, as must 
the isotope contribution, if ‘‘absolute” values are 
desired; for hydrogen bromide, the spin con- 
tribution would be R In 8. 


*Gordon and Barnes, J. Chem. Phys. 1, 297 (1933). 

* W. G. Brown, Phys. Rev. 38, 1179 (1931), also reference 
2; Weizel, Bandenspektren, p. 385, Leipzig (1931). 

. Czerny, Zeits. f. Physik 44, 235 (1927); Kratzer, 
ne f. Physik 3, 289 (1920); Weizel, Bandenspektren, 


" Urey and Rittenberg, J. Chem. Phys. 1, 137 (1933). 
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TABLE I, 
Bromine Hydrogen Bromide 
petit 
T T 

200 47.822 55.291 8.272] 37.786 44.703 6.959 
250 49.509 57.162 8.486} 39.330 46.256 6.961 
298.1 50.868 58.667 8.618 40.550 47.481 6.964 
300 50.918 58.722 8.622} 40.594 47.526 6.964 
350 52.130 60.058 8.713| 41.663 48.599 6.970 
400 53.196 61.226 8.777| 42.589 49.530 6.983 
450 54.146 62.263 8.824| 43.407 50.354 7.006 
55.006 63.194 8.859} 44.139 51.094 7.039 

550 55.789 64.040 8.888| 44.802 51.767 7.085 
600 56.509 64.814 8.911| 45.409 52.386 7.141 
650 57.176 65.528 8.931 45.968 52.960 7.204 
700 57.797 66.191 8.948| 46.487 53.496 7.274 
750 58.377 66.809 8.963} 46.971 54.001 7.350 
800 58.922 67.388 8.977| 47.426 54.477 7.426 
850 59.437 67.932 8.989| 47.854 54.930 7.503 
900 59.923 68.446 9.001} 48.259 55.361 7.580 
950 60.384 68.933 9.011} 48.644 55.773 7.656 
1000 60.823 69.396 9.022} 49.010 56.168 7.730 
1050 61.242 69.836 9.031 49.360 56.547 7.801 
1100 61.643 70.257 9.041} 49.695 56.911 7.869 
1150 62.026 70.659 9.050} 50.016 57.263 7.934 
1200 62.394 71.044 9.059} 50.325 57.602 7.996 
1250 62.747 71.414 9.068} 50.623 57.929 8.055 
1300 63.087 71.770 9.077| 50.911 58.247 8.111 
1400 63.732 72.443 9.094} 51.456 58.851 8.214 
1500 64.334 73.071 9.113} 51.969 59.421 8.307 
1600 64.898 73.660 9.132 52.451 59.960 8.390 


It is, perhaps, of interest to note what error 
would have been introduced in Table I if the 
second order terms in the energy level formula 
had been ignored. If, for 1600°K, bromine had 
been treated as a harmonic oscillator of frequency 
w(1—2x)=321.1 plus a rigid rotator (B 
=(.080,77), the entries! in Table I would have 
been 64.826, 73.494 and 8.927. In the case of 
—(F°—E,°)/T, neglect of the anharmonic term 
in x causes an error of 0.036 in the final result, 
neglect of the interaction terms which take 
account of the change in B with », an error of 
0.022 and neglect of the stretching terms 
involving D, an error of 0.014. In the case of the 
heat capacity at the same temperature, the 
corresponding errors are 0.112, 0.060 and 0.030, 
respectively, which (with the contribution 0.004 
from the excited *II state) just account for the 
difference between the value 8.927 and the entry 
in the table. In the case of hydrogen bromide at 


12Using the notation of reference 8 for bromine at 
1600°K, we find that (In Q,+«0hc/kT)=1.4055, S,/R 
=2.2941, C,/R=1.0493, v=3.16, #®=23.9, s,=6.98, 
S.=78.9, =8.20, c = 186; =0.003,405, 6. =1.1610-, 
d =d,(1+0.0068%) = 0.00709. 
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1600°, the assumption that its spectrum could be 
treated as that of a harmonic oscillator of 
frequency 2559 plus a rigid rotator with B 
= 8.342, would lead to the values 52.418, 59.884 
and 8.260 instead of the numbers entered in the 
table. 


THE THERMAL DISSOCIATION OF BROMINE 


From Brown’s value of the heat of dissoci- 
ation,® 1.962 volts, AE,° for the reaction Bre 
=2 Br is 45,230 cal. The values of —(F° — E°)/T 
for monatomic bromine can be readily computed 
since the ‘State Sum”’ for monatomic bromine is 
simply (4+2-exp [—3685hc/kT ]) where 3685 
cm is the separation” of the ?P; and ?P; states 
of the bromine atom. The resulting —R In K, 
=R In (Pp;,)/(Psr)? are entered in Table II and 
correspond to the continuous curve of Fig. 1; the 
figure also shows the experimental results of 
Bodenstein® and of Perman and Atkinson" as 
recalculated by Lewis and Randall. It is 
evident that there is still a difference of about 


TABLE II. Ki=(Pp,)*/(Ppr,) and 


Br2(g) =2Br 3H2+3Br2(g) =HBr(g) 
AF° —AF 


T°K —RinK;, kcal. Rin kcal. 
298.1 128,92 38.43 43.12 12.85 
300 127.95 38.39 42.86 12.86 
350 106.09 37.13 37.09 12.98 
400 89.68 35.87 32.75 13.10 
450 76.89 34.60 29.36 13.21 
500 66.65 33.33 26.63 13.32 
550 58.26 32.04 24.40 13.42 
600 51.26 30.76 22.53 13.52 
650 45.33 29.47 20.94 13.61 
700 40.25 28.17 19.58 13.70 
750 35.83 26.87 18.39 13.79 
800 31.96 25.57 17.35 13.88 
850 28.55 24.27 16.43 13.97 
900 25.51 22.96 15.61 14.05 
950 22.78 21.64 14.88 14.13 
1000 20.33 20.33 14.21 14.21 
1050 18.11 19.02 13.61 14.29 
1100 16.09 17.70 13.06 14.37 
1150 14.24 16.38 12.57 14.45 
1200 12.54 15.05 12.11 14.53 
1250 10.97 13.72 11.68 14.61 
1300 9.53 12.39 11.30 14.68 
1400 6.94 9.72 10.60 14.84 
1500 4.69 7.04 9.99 14.99 
1600 2.71 4.34 9.46 15.14 


13 Turner, Phys. Rev. 27, 397 (1926). 

4 Perman and Atkinson, Zeits. f. physik. Chemie 33, 
215, 577 (1900). 

8 Lewis and Randall, Thermodynamics, p. 514 (1923). 
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Fic. 1. The equilibrium constant for the thermal dis- 
sociation of bromine. K:i=(Ppr)?/(Ppr;). 


half a calorie per degree between Bodenstein’s 
results and the entries in Table II, and that this 
discrepancy is considerably larger than his 
apparent experimental error. This result is rather 
surprising, since there is no reason to distrust the 
spectroscopic data and Bodenstein’s figures give 
every internal evidence of self-consistency. 


THE DISSOCIATION OF HYDROGEN BROMIDE 


From the entries of Table I, Giauque’s data 
for hydrogen'® and the calorimetrically deter- 
mined heat of formation of gaseous hydrogen 
bromide, viz., 12,100 cal., it follows that AF°s9.,; 
for thereaction }H2+43Br2(g) = HBr(g) is— 12,850 
cal. and that is — 11,990 cal.!” The resulting 
values of R In Ke=R In (Pusr)/(Pu,)*: (Par,)! 
and of AF° for the reaction are given in Table II; 
these numbers agree with those given by Nernst’s 
free energy equation!® (which was based on the 
e.m.f. measurements of Bodenstein and Geiger)" 
within about a quarter of a calorie per degree in 
R In Ke up to 1500°K. 

Direct measurements on this equilibrium have 
been carried out by von Falckenstein,® the 
partial pressure of the hydrogen in the equilib- 


16 Giauque, J. Am. Chem. Soc. 52, 4816 (1930). 

1 Datta, Zeits. f. Physik 77, 404 (1932), gives 87,700 as 
the heat of dissociation of hydrogen bromide from measure- 
ments on continuous absorption. Combining this value 
with the heat of dissociation of hydrogen (102,900 cal., 
according to Richardson, Proc. Roy. Soc. A123, 466 (1929)) 
and of bromine, AEo° is — 13,600 cal. Datta’s value, how- 
ever, is based on a considerable extrapolation. 

18 Nernst, Zeits. f. Elektrochemie 15, 687 (1909). 

19 Bodenstein and Geiger, Zeits. f. physik. Chemie 49, 
70 (1904). 
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rium mixture being determined by means of a 
semi-permeable platinum diaphragm. He found, 
as mean values for the forward and reverse 
reaction, Py, to be 1.91, 2.78 and 4.17 mm at 
1297°, 1381° and 1495°, respectively,”° with total 
pressures of 761, 764 and 777 mm. For these 
three temperatures, 1/K, is 126, 41.5 and 11.2 
from Table II, so that Px,, is 0.807, 0.830 and 
0.707 mm and Pye, is 756, 756 and 765 mm. 
The resulting R In K» for the three temperatures 
are 12.74, 12.34 and 12.11, respectively, in 
definite disagreement with the values found by 
interpolation from Table IJ, viz., 11.32, 10.72 and 
10.02. Thus the supposed agreement between 


* The dissociation of diatomic into monatomic hydrogen 
is negligible at these temperatures; see reference 16. 


von Falckenstein’s results and those of Boden- 
stein and Geiger comes solely from the fact that 
the calculations based on von Falckenstein’s 
measurements ignored the dissociation of dia- 
tomic into monatomic bromine. It should be 
noted in passing that Giauque and Overstreet® 
found a similar disagreement in the case of 
Léwenstein’s measurements”! on the dissociation 
of hydrogen chloride. The entries in Table II are 
of course uncertain to the extent of the error in 
the calorimetrically determined heat of formation 
of hydrogen bromide but von Falckenstein’s R In 
Ko, after correction for the bromine dissociation, 
cannot be made self-consistent by any adjust- 
ment in the heat of formation. 


*1 Lowenstein, Zeits. f. physik. Chemie 54, 715 (1906). 
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The Reactions of the Hydroxyl Radical in the Electrodeless Discharge in Water Vapor 


W. H. RopesusH AND M. H. Want, University of Illinois 


In a study of the reactions occurring in the electrodeless 
discharge in water vapor consistently high yields of hydro- 
gen peroxide were obtained. A number of reactions have 
been studied and by a process of elimination a mechanism 
for the formation of hydrogen peroxide from hydroxyl has 
been suggested. In the course of some spectrographic 


(Received June 19, 1933) 


some new bands, with heads at 3564 and 3328A, were 
observed which are probably due to an ionized hydroxyl. 
These bands were not present in the uncondensed electrode 
discharge but appeared in the condensed electrode dis- 
charge at low pressures. 


HEN a discharge is passed through water 

vapor the so-called ‘‘water’’ bands are 
excited. These bands are known to be due to 
the hydroxyl radical and it is supposed that an 
electron, having an energy of about nine volts, 
dissociates a water molecule into a hydrogen and 
an excited hydroxyl, which then emits the char- 
acteristic radiation, 


HOH =H+OH. (1) 


Even with the weakest discharge the hydroxyl 
bands are emitted with a surprising intensity. 
This would indicate the possibility of obtaining 
a high concentration of hydroxyl and numerous 
attempts have been made to study the reactions 
of this molecule. In this laboratory an attempt 
was made to obtain a molecular beam of hy- 
droxyl but the attempt failed because no chem- 
ical reaction characteristic of hydroxyl was dis- 
covered and it was therefore uncertain as: to 
whether hydroxy] itself were present in the beam. 
Abroad, Steiner! and Bonhoeffer? made similar 
attempts to study the chemical behavior of 
hydroxyl, likewise with negative results. Finally 
Bonhoefier and Pearson? published an account of 
an unsuccessful attempt to obtain the hydroxyl 
bands in absorption. They reached the conclusion 
that the hydroxyl concentration was extremely 
low because of the rapid disappearance of hy- 


1Steiner and Bay, Zeits. f. physik. Chemie 3B, 149 
(1929). 

? Bonhoeffer and Pearson, Zeits. f. physik. Chemie 14B, 
2 (1931). 
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droxy! by the reaction, 
OH+0OH =H,0+0. (2) 


It occurred to us that if this reaction takes place, 
a high concentration of atomic oxygen should be 
present in the vapor which has passed through 
the discharge and it should therefore be possible 
to obtain a molecular beam of atomic oxygen 
which could be identified by well-known meth- 
ods.* In carrying out this result however we 
discovered that there was very little oxygen, 
either atomic or molecular, present in the dis- 
charge and that large yields of hydrogen peroxide 
were obtained, whereas Bonhoeffer and Pearson 
had reported only traces of hydrogen peroxide 
to be formed. This made it appear doubtful that 
reaction (2) was occurring and made necessary a 
careful study of the reactions actually taking 
place in the water vapor discharge. 


EXPERIMENTAL 


In the study of the chemical reactions the 
electrodeless discharge was used. The apparatus 
by which the discharge was produced was essen- 
tially the same as that described in a previous 
article.t A flow system was used in all the experi- 
ments. The water reservoir was made from 5 mm 
Pyrex tubing which had previously been cali- 
brated so that the amount of water used during 
a run could be measured. The height of the water 
column was followed with a cathetometer. The 

3 Kurt and Phipps, Phys. Rev. 34, 1357 (1929). 


‘Rodebush and Klingelhoefer, J. Am. Chem. Soc. 55, 
130 (1933). 
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flow of water vapor was regulated by a capillary 
leak and by the temperature of the bath sur- 
rounding the reservoir. Conductivity water was 
always used in the reservoir. 

The discharge was produced in a 200-cc Pyrex 
flask. Optimum conditions for high yields of 
hydrogen peroxide were at pressures of 0.1 to 
0.2 mm of mercury. Under best conditions the 
discharge had a vivid red appearance, but if the 
discharge was too intense blue rings appeared, 
which followed the turns of the coil, indicating 
the presence of oxygen atoms. This was con- 
firmed by spectrographic evidence. A_ better 
yield of hydrogen peroxide was obtained when 
the discharge bulb was cooled by a blast of air. 

It was found expedient to use a U-tube trap 
for collecting the sample. In the early course of 
the work ground glass joint traps, which necessi- 
tate the use of stopcock grease, were employed 
and it was found impossible to obtain con- 
sistently large yields of hydrogen peroxide. A 
second U-tube trap was placed in series with the 
first to determine if all the sample was caught in 
the first trap. Two runs indicated that no con- 
densable products escaped through the first 
trap. 

After a run had proceeded the desired length 
of time and the sample frozen out in the liquid 
air-cooled trap, the apparatus was shut off from 
the pumps by a large stopcock. The vacuum was 
broken, drawing the air slowly through a liquid 
air trap so no moisture would be carried into the 
trap containing the sample. When the pressure 
within the apparatus was atmospheric the sample 
was permitted to melt. The trap was then cut off, 
weighed, and the sample titrated with 0.1 N. 
potassium permanganate. After a thorough rins- 
ing the trap was dried and reweighed. The 
sampling process was carried out as hurriedly as 
possible in order to minimize decomposition of 


peroxide. A set of representative data is given in 
Table I. 


TABLE I, 

Millimols Mol % water con- Mol % water con- 
water verted to hydrogen verted to hydrogen 
used peroxide and hydrogen and oxygen 

1 18.67 50.5 4.9 

2 18.50 50.8 4.6 

3 18.77 48.8 2.0 

4 20.30 58.2 2.6 
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Effect of surface 


The reaction tube was filled with Pyrex tubing 
so that the surface-volume ratio was increased 
more than fivefold. There was no increase in 
yield of peroxide, showing that the peroxide 
forming reaction is not influenced by the surface 


of the confining vessel. Results are given in 
Table II. 


TABLE II. 
Surface 
volume Time Wt.sample Wt.H,O. in 
(cm) (min.) (g) (g) sample 
2.09 75 0.60 0.318 53.0 
11.80 60 0.45 0.221 49.0 


From time to time during the course of this 
work it was noticed that the peroxide yields were 
lowered whenever the glass became very warm. 
The reaction tube was placed in an ice bath and 
several runs made to determine if decomposition 
occurred when the wall was at room tempera- 
ture. There was no variation in the yield of 
peroxide. 


Effect of distance 


A study was made of the yield of hydrogen 
peroxide as a function of the distance between 
the discharge bulb and the trap. The results are 
shown in Table III. 


TABLE III. 

Distance Wt. % Distance Wt. % 
(cm) in sample (cm) in sample 
160 50.2 17 50.0 
100 50.7 3 52.0 

40 51.6 a 47.0 


* This run was made with a very high rate of flow. 


Action of catalysts 


By the action of a silver screen catalyst, an 
attempt was made to destroy the active species 
at varying distances down the reaction tube and 
determine what effect it would have on the 
peroxide yields. As soon as the discharge was 
turned on the catalyst became hot and some 
metal evaporated to the walls of the tube, causing 
it also to be catalytically active. Analysis of the 
sample showed only a very small fraction of a 
percent of peroxide present. It seemed plausible 
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that, even if the peroxide were formed before it 
reached the catalyst, it would be decomposed by 
the hot metal. A concentrated solution of hy- 
drogen peroxide (85 percent) was prepared by 
the method given by Rice, Reiff and Kilpatrick.* 
The vapor from this solution was passed through 
the apparatus, which was entirely free of any 
catalyst, and a sample collected and analyzed. 
This procedure was repeated with the catalyst in 
the reaction tube. The sample in the first case 
tested about 65 percent by weight of peroxide 
and in the latter only 3 percent, showing that 
practically a quantitative decomposition of per- 
oxide was brought about by the catalyst. A 
chemically deposited silver mirror was found to 
be equally effective for the decomposition of per- 
oxide. 

Taylor and Lavin® found that dehydrogenation 
catalysts, such as zinc chromite, were selective 
for atomic hydrogen recombination but did not 
destroy the hydroxyl molecules. The discharge 
bulb and reaction tube were coated with a 
similar dehydrogenation catalyst, copper chro- 
mite, and a discharge run on water vapor. The 
collected sample showed no peroxide present. 

It would be very desirable to know whether 
the hydrogen peroxide is formed in the vapor 
immediately after the latter leaves the discharge 
or only by condensation in the liquid air trap. 
A possible method for detecting hydrogen per- 
oxide vapor is the photography of the ultra- 
violet absorption spectrum. The absorption spec- 
trum of hydrogen peroxide vapor was found by 
Urey, Dawsey and Rice’ to be continuous in the 
region from 2150 to 2750A. An attempt was 
made to obtain the absorption spectrum of 
peroxide as formed in the water vapor discharge. 
The discharge bulb was placed near a long ab- 
sorption tube but arranged so the discharge 
would not be photographed. Quartz windows 
were sealed on each end of the absorption tube. 
A small hydrogen arc was used as a source of 
continuum and was found to give consistently 
reproducible intensities. Several plates were 
taken with the discharge in operation but no 


5 Rice, Reiff and Kilpatrick, J. Am. Chem. Soc. 48, 3019 
(1926). 

6 Taylor and Lavin, J. Am. Chem. Soc. 52, 1910 (1930). 

7Urey, Dawsey and Rice, J. Am. Chem. Soc. 51, 1371 
(1929). 
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noticeable absorption was obtained. A blank run 
was made passing the vapor from a concentrated 
solution of peroxide through the absorption tube 
and several pictures were taken. A very small 
amount of absorption was indicated. The partial 
pressure of peroxide in these blank runs was 
about 10 times greater than that occurring in 
the water vapor discharge. It is evident that the 
partial pressure of peroxide as found in the water 
vapor discharge is too small to give any notice- 
able absorption. 

Since it was not possible to determine either 
exactly where the hydrogen peroxide is formed or 
what the mechanism is by which it is formed 
there remained only the alternative of settling 
the mechanism by a process of elimination. 
Accordingly a number of chemical reactions 
which might produce hydrogen peroxide were 
tried out. 


Condensation of oxygen atoms with water in the 
trap 

Water vapor was rapidly drawn through the 
apparatus and frozen out in the trap, taking 
care to obtain a reasonably uniform layer of 
water on the walls of the trap. A discharge was 
then run with oxygen, at about 0.2 mm pressure, 
and a sample collected over a period of 60 
minutes. The sample gave no reaction with 
potassium permanganate indicating the absence 
of hydrogen peroxide. It did have an ozone odor 
and colored a solution of potassium iodide. 
Three different trials indicated the same results. 

The apparatus was changed so that the dis- 
tance between the discharge bulb and the trap 
was about 1.6 meters. In traversing this distance 
the oxygen atoms would have completely dis- 
appeared, so that there would be no possibility 
of a reaction involving oxygen atoms in the trap. 
A run of some hours’ duration was made with 
the water vapor discharge and a sample collected. 
Titration with permanganate showed the sample 
to be 50.2 percent by weight peroxide. The 
possibility of a condensation reaction in the trap 
involving oxygen atoms appears to be excluded. 


Reaction of hydrogen atoms with oxygen 


The oxygen was introduced into the active 
species from the hydrogen discharge a short F7 
distance from the discharge bulb, taking pre [7 
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cautions to prevent back diffusion of the oxygen 
into the bulb. The spectrum of the discharge did 
not change on admittance of the oxygen into the 
reaction tube so it was assumed that back 
diffusion was negligible. The flow of gases was 
measured by gas burettes and regulated by 
capillary leaks. Results are shown in Table IV. 


TABLE IV. 
Flow of gas 
Time (cc./min. at Wt. Wt. % H.02 
atmos. conditions) sample H.O, in 
No. (min.) O» (g) (g) sample 


1 98 22.0 1.30 0.70 0.022 3.1 
Z 80 15.5 3.70 0.06 0.0002 0.33 
3 110 6.3 1.15 0.06 0.039 64.5 
4 60 14.0 1.10 0.05 0.020 40.0 


At the higher rates of flow the yields were very 
noticeably decreased, but, the reaction tube in 
the region of the entrance of the oxygen became 
very warm. It was consistently observed in the 
course of this work that whenever the glass 
became heated the yield of peroxide was always 
decreased. 


Reaction of oxygen atoms with hydrogen 


Several trials were made, employing different 
rates of flow, but in no case was any appreciable 
amount of hydrogen peroxide present. Usually 
one to three drops of standard permanganate 
would give a pink coloration. Good yields of 
water were obtained and a small amount of 
ozone was noted. The heating effect in the re- 
action tube was quite small. 


Reaction of oxygen atoms with water 

Runs were made varying the pressure of 
oxygen in the discharge bulb and the rate of 
flow of water vapor but no trace of peroxide was 
found in any of the samples. Because of the small 
sample obtained runs of several hours duration 


were made. Small amounts of ozone were de- 
tected. 


Reaction of hydrogen atoms with water 


No evidence of any sort of a reaction was 
noticed. No peroxide was found in the samples 
and there was no heating effect at the entrance 
of the water vapor into the active species from 
the hydrogen discharge. 
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vapor discharge and oxygen 


A discharge was run on the water vapor and 
the oxygen added to the active species in the 
reaction tube. The results are given in Table V. 
The weight of the sample is larger than the 
amount of water used in each measurable case. 
Beside the reactions occurring in the water 
vapor discharge the reaction between hydrogen 
atoms and oxygen molecules occurs, and it was 
found to yield both peroxide and water. 


TABLE V. 

Wt. H.02 Wt. Wt. 

Time used Flow sample % 

No. (min.) (g) (cc. /min.) (g) in sample 
1 90 0.83 0.160 53.1 
2 90 2.20 0.160 47.2 
3 45 0.372 1.13 0.410 55.3 
4 60 0.532 1.13 0.588 55.0 


Effect of premixing hydrogen and water vapor 


The hydrogen was added to the water vapor 
before passing into the discharge bulb. Two 
trials were made and in each case the yield was 
about the same as in the water vapor discharge. 
(See Table VI.) This indicates that there is very 
little oxygen present with which the hydrogen 
may react. 


TABLE VI. 
Flow of Wt. H.O Wt. Wt. % 
Time He used sample H,0, in 
(min.) (cc./min.) (g) (g) (g) sample 
60 1.50 0.150 0.070 46.6 


60 1.86 0.436 0.421 0.206 48.8 


Effect of premixing oxygen and water vapor 


Several runs were made and in each case the 
sample, which was larger than the amount of 
water used, tested about 50 percent peroxide. 
It appears that neither oxygen atoms nor mole- 


TABLE VII. 
Flow of Wt.H.O Wt. Wt. 
Time used sample H.O, in 
No. (min.) (cc./min.) (g) (g) (g) sample 
1 40 2.30 0.12 0.059 49.7 


2 45 1.34 0.336 0.362 0.164 45.3 
3. «640 0.715 0.337 0.350 0.187 51.6 
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Reaction of the active species from the water 
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700 W. H. RODEBUSH 
cules participate in the peroxide forming re- 
action, for if they did the yield of peroxide should 
have been noticeably increased. (See Table VII.) 


Spectrographic studies 

The spectrum of the electrodeless discharge in 
water vapor was photographed with a Hilger 
quartz spectrograph (E38) with Eastman 40 
plates. An exposure of a few seconds gave the 
well-known water bands at 3064 and 2811A and 
in addition some hitherto unreported bands. The 
strongest of these had a head at 3564A and was 
degraded toward the red. A fainter band had a 
head at 3328A. In order to ascertain if possible 
why these bands had not been previously ob- 
served the emission spectrum of water vapor 
was photographed in a Wood tube with alu- 
minum electrodes. With a 60 cycle discharge no 
trace of the bands was obtained but with a 
condensed discharge the new bands came out 
very strongly at pressures in the neighborhood 
of 0.1 mm. At higher pressures they disappeared. 

Several plates were made of the bands on a 
large spectrograph with a dispersion of 3A per 
mm. Preliminary assignments were made but the 
analysis proved to be very difficult and was not 
carried further. The bands have a hydride 
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Fic. 1. Electrodeless discharge. 0.1 mm pressure. 


structure and since they appear only in the con- 
densed discharge are to be attributed to an 
ionized radical. As water of high purity was used 
to supply the vapor it seems probable that the 
carrier is OH*. 

These bands always appear with a condensed 
discharge at pressures low enough for the dis- 
charge to have a red color. At higher pressures 
the discharge is bluish white and these bands are 
not obtained. 

The spectrum of the electrodeless discharge at 
low pressures, taken on a small Hilger spectro- 
graph is shown in Fig. 1. In Fig. 2 the same 
spectrum is shown with large dispersion. The 


Fic. 2. Electrodeless discharge. 0.1 mm pressure. 


new bands at 3564 and 3328A are easily seen as 
well as the well-known hydroxyl] bands. Figs. 3, 
4 and 5 are photographs of the discharge in 
water vapor in a Wood tube, with aluminum 
electrodes. The new bands appear only at low 
pressures and when a condensed discharge is 
used. 


8 Professor F. W. Loomis of the Department of Physics 
has undertaken the analysis of these bands. A preliminary 
note, reporting their discovery was published; J. Am. 
Chem. Soc. 55, 1742 (1933). 


Stewart and Lavin® have reported that they 
obtained the OH bands in emission apparently 
as an afterglow from the dark part of the tube 
30 cm from the discharge. This point was in- 
vestigated and it was found that under ordinary 
conditions a faint discharge travelled throughout 
the apparatus which gave the hydroxyl bands 
with considerable intensity. When elaborate pre- 
cautions were taken in the way of shielding to 


* Stewart and Lavin, Proc. Nat. Acad. Sci. 15, 829 (1929). 
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Fic. 3. High frequency electrode discharge. 
0.1 mm pressure. 


Fic. 4. High frequency electrode discharge. 
Several mm pressure. 


Fic. 5. 60 cycle electrode discharge. 0.1 mm pressure. 


prevent the discharge, no trace of the hydroxyl 
bands appeared after fifteen minutes exposure on 
the small spectrograph. 


DISCUSSION 


The reaction mechanism proposed by Bon- 
hoeffer and Pearson to account for the dis- 
appearance of hydroxyl cannot be correct since 
large amounts of hydrogen peroxide were ob- 
tained in these experiments. The difference in the 
experimental results is undoubtedly due to the 
fact that the electrodeless discharge was used 
and the powerful catalytic action of traces of 
sputtered metal avoided. 


The failure of Bonhoeffer and Pearson to 
obtain the absorption spectrum should not be 
difficult to explain. On the other hand we are 
inclined to agree with their conclusion that 
hydroxyl disappears by some one fast reaction 
mechanism so that its concentration is always 
low and the molecule is not available for other 
reactions. A low concentration of hydroxyl is 
not inconsistent with a high intensity of emission 
of the hydroxyl bands since it is the water mole- 
cule which is excited to give the hydroxyl bands. 
We suspect that the reaction by which hydroxyl 
disappears is the rapid homogeneous bimolecular 
condensation to hydrogen peroxide, 


OH+0OH =H,0.. (3) 


The energy liberated here is about 41 large 
calories. This amount of energy would corre- 
spond to the excitation of the hydroxyl to the 
fifth vibrational state and the disposal of this 
energy does not seem to be a serious obstacle to 
the reaction. 

Heterogeneous reactions do not occur rapidly 
unless adsorption takes place upon the walls of 
the vessel and experience has shown that mini- 
mum of catalytic wall activity is obtained when 
considerable amounts of water vapor are present. 
Since large concentrations of hydrogen atoms 
were shown to reach the liquid air trap the wall 
activity was certainly slight. The only alternative 
to bimolecular reactions are those taking place 
in triple collisions and these must of necessity 
be slow reactions. 

The above reaction mechanism should be very 
rapid, causing the disappearance of hydroxyl 
immediately after leaving the discharge. About 
50 percent of the water is converted to hydrogen 
peroxide. It is not likely that much more than 
50 percent of the water molecules are dissociated 
in the first place and there are reactions which 
could account for the reconversion of consider- 
able amounts of hydrogen peroxide to water. 
With a more intense discharge a considerable 
amount of atomic oxygen is produced probably 
by the dissociation of free hydroxyl. But at 
moderate intensities the amount of oxygen pro- 
duced is negligible. The results of Bonhoeffer and 
Pearson who obtained only traces of hydrogen 
peroxide are accounted for by the presence of 
particles of metal sputtered from the electrodes. 
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Our experiments show that hydrogen peroxide 
vapor is decomposed quantitatively by small 
amounts of exposed metal surface into water 
and oxygen. 

The evidence for the above reaction is negative 
and falls into two classes. First, one must 
assume that if hydrogen peroxide were not 
formed by simple, rapid, mechanism that com- 
peting reactions would take place. The most im- 
portant of these is the reaction with atomic 
hydrogen. 

H+0OH =H,0—115 cal. (4) 


The energy to be disposed of here is probably too 
large to permit a bimolecular mechanism. 

Finally it is necessary to consider other possible 
reactions by which hydrogen peroxide may be 
formed. It is necessary to consider only those 
mechanisms which are sufficiently simple to be 
inherently probable. There are two mechanisms 
which form hydrogen peroxide directly. 


O+H,.0 =H,0,—33 cal. (5) 
OH+H,0=H,0.+H+72 cal. (6) 


The first reaction is inherently improbable as 
it involves a change from a triplet to a singlet 
state and it was shown in this work not to take 
place either in the main part of the apparatus or 
in the liquid air trap. The second reaction cannot 


take place because of the large amount of heat 


absorbed. 

It was shown in this work and has been shown 
previously by Bates!® and others that atomic 
hydrogen reacts with molecular oxygen to give 
hydrogen peroxide by some mechanism which is 
presumably complex and hence not very rapid. 
There are two mechanisms by which molecular 
oxygen might be formed. 


OH+0OH =H:+0, (7) 
OH+OH =H,0+0. (8) 
0+0 =O, 


The first reaction has a high heat of activation 


10 Bates and Salley, J. Am. Chem. Soc. 55, 110 (1933). 
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since the reverse reaction has a high heat of 
activation and the heat of reaction is small. The 
second mechanism is controverted by the experi- 
mental evidence as to the concentration of atomic 
oxygen. Furthermore the formation of molecular 
oxygen by these two steps must be exceedingly 
slow. Also to produce a 50 percent conversion to 
hydrogen peroxide 100 percent dissociation of 
water molecules would be necessary. There is 
reason to believe as has been mentioned that the 
conversion to hydrogen peroxide is much greater 
than observed since the reaction 


H+H,0, =OH+H,.0 


must take place to some extent particularly on 
the walls of the trap. 

The experiment with the discharge only three 
cm from the trap with a flow rate of 9 meters 
per second shows that the reaction probability 
per collision from hydroxyl by (3) must be of 
the order 10~*. It is unlikely that any competing 
reaction can be found which is as fast. 

There is of course no direct proof that the 
condensation of hydroxyl does not take place on 
the walls of the trap and undoubtedly this re- 
action does take place to some extent. The fact, 
however, that no variation in yield of peroxide 
with distance of trap or rate of flow was ob- 
served, is most easily explained by assuming the 
disappearance of hydroxyl by a gas phase re- 
action which is faster than any competing re- 
action. 

If we assume that the new bands observed are 
due to OH* then the fact that they appear only 
at low pressures tends to confirm the hypothesis 
of the reaction mechanism. While the absolute 
amount of hydroxyl present will decrease with 
decreasing pressure the relative amount will in- 
crease and this will result in the formation of 
an increased amount of OH+. This conclusion 
can only be true if the disappearance reaction for 
hydroxyl is of second or higher order. 

The authors wish to express their gratitude to 
Professor F. W. Loomis for the loan of spectro- 
graphic equipment and for advice in the interpre- 
tation of spectra. 
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Complex Formation due to Polarization. I. The System Krypton and 
Hydrogen Chloride 


GEORGE GLOCKLER, CHARLES P. RoE! AND DONALD L. FULLER, University of Minnesota 
(Received July 10, 1933) 


Rare gases form compounds of a loose complex nature 
such as hydrates, solvates or mercury complexes. Such 
molecular aggregates contain a molecule of permanent 
dipole moment or an excited atom or molecule of dipole 
character as the other reaction partner besides the rare 
gas. The forces holding the complex molecule are the 
attraction between the permanent dipole and the dipole 
induced in the rare gas molecule due to its polarizability. 
Five lines of evidence are cited which indicate the forma- 
sion of such complex molecules: (1) Band spectra of complex 


HgA and HgKr are known. (2) Rare gas hydrates have 
been described. (3) Isoelectronic systems similar to the 
rare gases have proton affinity. (4) Rare gas hydride ions 
have been found in the mass-spectrograph. (5) The P-V-T 
relation of gaseous mixtures show that interaction between 
the unlike molecules exists. The second Virial coefficient 
is determined experimentally for Kr—HCl mixtures as a 
function of composition and it is shown that interaction 
exists between the rare gas atom Kr and the permanent 
dipole HCl. 


INTRODUCTION 


HE rare gases are usually considered to be 
chemically nonreactive although period- 
ically the preparation of some chemical com- 
pounds involving rare gases is announced.? In 
this paper we wish to discuss this question and 
make clear the distinction between “primary 
valency compounds” and ‘‘complex molecules” 
involving the rare gases. Whenever the possibility 
exists for the formation of an electron pair,* then 
the usual valency forces are acting. However, the 
further possibility exists that atoms which cannot 
react in the normal state with other atoms or 
molecules may still form complex structures on 
account of the van der Waals force field4 which 
surrounds them. It will be shown that hydrate 
formation and the behavior of rare gas molecules 
toward dipole structures in general give evidence 
of such loosely bound molecular complexes held 
together by polarization forces. 


‘This article is based upon a thesis presented to the 
faculty of the Graduate School of the University of 
Minnesota by Charles P. Roe in partial fulfillment of the 
requirements for the degree of Doctor of Philosophy. 

*Gmelin-Kraut, Handbuch Inorg. Chem. 8th Edition. 
Verlag Chemie, Leipzig-Berlin (1926). Mellor, Inorg. Chem. 
Vol. VII, Longmans, Green & Co., London (1927). 

*G. N. Lewis, Valence and the Structure of Atoms and 
Molecules. Chemical Catalog Co., N. Y. (1923). Heitler 
and London, Zeits. f. Physik 44, 455 (1927). 

*F, London, Zeits. f. physik Chemie 11B, 222 (1931). 


One of the latest attempts to produce chemical 
compounds involving a rare gas has been de- 
scribed by Antropoff’ who has studied the 
reaction between chlorine and krypton in elec- 
trical discharge. He believes that a compound 
“krypton chloride’ is possible. That such a 
compound cannot be held together in the 
ordinary way by primary valency forces, may be 
shown quite definitely by the use of thermo- 
dynamic arguments involving a Born-Haber 
cycle. Grimm and Herzfeld® show that hypo- 
thetical KrCl is endothermic to the extent of 
95,000 calories. If it be supposed that the 
compound be salt-like with chlorine as the 
positive ion in the crystal-lattice, then it can be 
shown that the hypothetical “‘chlorine kryptide” 
would be endothermic to the extent of 165,000 
calories, if the electron affinity of Kr is zero and 
the grating energy is that of RbCl. Either of the 
two structures considered seems then energetic- 
ally impossible mainly because of the high 
ionization potentials. If the two elements krypton 
and chlorine really combine, then the resulting 
compound must have a different structure from 
the salt-like configurations assumed so far. The 
bond of the hypothetical molecule KrCl may be 
made between an excited state of krypton and 


5 A. von Antropoff, Naturwiss. 20, 688 (1932). 
6 Grimm and Herzfeld, Zeits. f. Physik 19, 141 (1923). 
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the normal state of chlorine, and the resulting 
energy rich molecule may at least have transient 
existence. Compounds of the rare gases of the 
ordinary valence type are then not to be expected 
and we shall next consider the evidence for more 
loosely bound structures or complex molecules of 
a van der Waals type. 


EVIDENCE FROM BANDSPECTRA 


Complex molecules which have their origin 
in the combination of one excited and one 
unexcited atom are well known in band spectro- 
scopy. 

Oldenberg’ has studied band spectra excited by 
fluorescence with \2537A in mixtures of rare 
gases and mercury vapor, and he explains these 
spectra by means of the assumption that mole- 
cule formation is possible not only between 
normal rare gas atoms and excited mercury 
atoms 


A('S) +Hg(?P)-A: Hg*, 


but that normal atoms of the two species con- 
sidered can also combine into loosely bound 
structures of the type HgKr and HgA: 


A(‘S)+Hg('!S) ~HgA; AH =0.025 e.v. 
Kr('!S)+Hg('!S) >HgKr; AH = 0.035 e.v. 


The small value of the heat of dissociation is 
readily understood, for the forces acting between 
these atoms are due to their mutual polarization. 
Oldenberg favors the view that the normal state 
of these molecules has a definite system of 
vibrational levels and a definite heat of dis- 
sociation. Another interpretation due to Condon® 
is to suppose that the normal state of the 
complex molecule A-Hg is essentially repulsive 
and that maxima of intensity in the true con- 
tinuum are really not bands in the ordinary sense 
but are due to the oscillatory behavior of the 
wave functions describing the system and its 
transitions. The possibility of such interaction 
exists because the excited mercury atom pos- 
sesses a high dipole moment, and it can therefore 
polarize the rare gas atom. The same type of 
interaction might then be expected between 
normal rare gas atoms and other normal atoms or 


7 O. Oldenberg, Zeits. f. Physik 55, 1 (1929). 
8 E. U. Condon, Phys. Rev. 32, 858 (1928). 
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molecules possessing permanent dipole moments 
such as hydrogen chloride and water resulting in 
the formation of a structure of greater com- 
plexity. 
RARE GAs HyYDRATES 

The formation of hydrates of gases was studied 
very early® and the heavier rare gases form such 
hydrates. It is interesting to note that the 
stability of the rare gas hydrates increases with 
increasing atomic weight and polarizability of the 
rare gas. 


ISOELECTRONIC SYSTEMS 


However, an entirely different point of view 
gives rise to the’same ideas. This point of view is 
based on octet theory of atomic structure and has 
been stated by Grimm.'® An atom and a mole- 
cule are said to be isoelectronic systems if they 
have the same outer electronic configuration. 
They are also called pseudo-atoms. For example 
hydroxyl ion (OH~) is a pseudo atom of fluorine. 
Table I gives such related structures. The 


TABLE I, Jsoelectronic systems. 


Na* 
Ne NeH* 
F HF H.F(?) 
O OH H,O H,0* 
N NH NH. NH; NH,* 
CH CH, CH; CH, CH,* 


formation of NH,* and H;0* follows very simply 
on the basis of these considerations. A study of 
the table immediately suggests that perhaps neon 
and methane may possess a slight proton affinity 
like water and ammonia. This entirely different 
view then leads also to an understanding of the 
existence of hydrates of the rare gases and it 
suggests possible complex structures between 
rare gas molecules and hydrogen chloride. This 
point of view would, however, lead one to con- 
sider these structures as “primary valency” 
compounds. 


APPLICATION OF QUANTUM MECHANICS 


Recently quantum mechanics has been applied 
to this problem with considerable success. The 


9M. de Frocrand, Comptes Rendus 135, 959 (1902); 
M. P. Villard, Thesis, Paris, 1896. 
10H, Grimm, Zeits. f. Electrochemie 31, 474 (1925). 
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simplest case of proton affinity has been studied 
by Unséld" who has given a theoretical solution 
of the problem: H+H+—-H,*. The problem of 
weak or van der Waals molecules has been 
considered by Ejisenschiitz and London” who 
give a perturbation method for the calculation of 
the mutual energy of two molecules at large 
separation. In this calculation the following 
forces are considered: (1) the mutual interaction 
from the electronic motions of the molecule; (2) 
the polarization of the electron shell of one atom 
by a permanent field of the other when the latter 
possesses electric asymmetry; and (3) the ordi- 
nary electrostatic forces between two asymmetric 
molecules. London‘ was thus able to calculate the 
order of magnitude of van der Waals a and to find 
the heats of such reactions as 


Hg+Hg—-Hg2, AH = 2000 cal. 
Hg+A—>HgA, AH = 1450 cal. 
Hg+Kr—-HgkKr, AH = 1720 cal. 


In all cases the results are in good agreement 
with the band spectra evidence mentioned above.’ 
Margenau"™ has extended the problem by calcu- 
lating interactions due to quadrupoles. 


MaAss-SPECTROSCOPY 


Consideration of Table I shows that the rare 
gas atoms may have a tendency to combine with 
charged ions to form more complex gaseous ions 
of the type HeH+, NeH?, etc. It is interesting to 
note that workers in the field of mass-spectro- 
scopy have found evidence of such structures. 
The following ions have been observed 
(HeH')*(5), (HeH?)*(6), (HeH'H')* (6), 
(NeH)*+(23) and (AH)*+(41). The experimental 
evidence favors the view that the proton is 
outside the nucleus because in the case of helium 
for example the mass found is 5.01 and not 5.00 


" A. Unsild, Zeits. f. Physik 43, 563 (1927). 

” Eisenschiitz and London, Zeits. f. Physik 60, 491 
(1930), 

“H. Margenau, Phys. Rev. 38, 747 (1931). 
ome R. Hogness and E. C. Lunn, Phys. Rev. 26, 44 


was Lukanow and W. Schiitze, Zeits. f. Physik 82, 610 
*K. T. Bainbridge, Phys. Rev. 43, 103; 44, 57 (1933). 


"H. Kallmann and W. Lasareff, Naturwiss. 20, 206 
(1932), 


which latter value would be expected if the ion 
considered were an isotope of helium. 


INTERMOLECULAR FORCES 


It was expected that some evidence of molecule 
formation could be obtained at high pressures by 
studying the pressure-volume-temperature re- 
lations of gaseous mixtures of krypton and 
hydrogen chloride. By such experimental studies 
we hope to contribute evidence for the existence 
of loose complex molecules involving the rare 
gases. A good summary of the theory of imperfect 
gases has been given by Fowler.'* According to 
this treatment the theory of molecular inter- 
action in a binary gas mixture leads to the 
following equation 


PV=RT+BP (for one mole of gas) (1) 


where 


(1 —x)?+2b,2(1 —x)x+Doox? (2) 
and 


f (3) 
Vap 


and P=pressure, 7=abs. temperature, R=gas 
constant, B=second virial coefficient, x= mole 
fraction of gas No. 1. 

E.s=interaction energy between molecule a 
and 6, where a and 6 may refer to the same or to 
different species. From these equations it is seen 
that bag will be zero if Eug is zero. In other words 
the existence of interaction between molecules of 
type a and type B can be determined by studying the 
second virial coefficient B and noting whether or not 
the interaction term by. is necessary for the de- 
scription of B as a function of composition. This 
study is made by investigating the pressure- 
volume-temperature relations of a gaseous mix- 
ture such as krypton and hydrogen chloride. 
Krypton was chosen because it has a high 
polarizability and no chemical reactivity in the 
ordinary sense; hydrogen chloride was chosen 
because it has a large permanent dipole moment. 


APPARATUS 


An apparatus has been designed for studying 
the pressure-volume-temperature relations of 


18R. H. Fowler, Statistical Mechanics, Ch. VIII, page 
171, Eq. (487). Cambridge University Press, 1929. 
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Te Gauge ToPomsey Burrito 


Fic. 1. General plan of experimental system. 


systems between —40°C and +250°C and at 
pressures from 1 to 200 atmospheres. A general 
outline is given in Fig. 1. An improvement over 
the usual Cailletet apparatus’? was made by 
connecting permanently the high pressure parts 
of the system to the supply and vacuum sections 
so that it was unnecessary to dismantle the 
capillary for each refilling. The hydrogen chloride 
is generated in flask A from H2SO, and NaCl, 
dried with P.O;, purified by sublimation and 
stored in flask F. The krypton was obtained from 
a commercial source and was not further purified. 
The pure or mixed gases were measured in the 
Ramsay burette K and run into the Toepler 
pump L. The needle valve U was closed and the 
large valve R was opened and the entire system 
to stopcocks d and e evacuated to 10~* mm Hg. 
Now the gas in L was led into the cylinder T 
(1000 cc) and capillary Q. After R was closed U 
was opened and the gas compressed with mercury 
as a confining liquid. The upper part of Y, which 
was connected directly to the gauge tester by 
means of the joint at V, was filled with glycerine; 
the lower end of Y was filled with mercury. This 
arrangement was necessary in order to prevent 
mercury from making contact with the gauge 
tester. The gauges which were calibrated with the 
piston tester are shown at X. 

In making the P—V—T measurements the 
pressure was noted when the mercury meniscus 
just reached one of the eight finely etched 


19K, Jellinek, Lehrbuch d. Phys. Chemie, Vol. I, pp. 
310-936. F. Enke, Stuttgart (1928). 
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calibrating marks on the capillary which was 40 
cm long and 0.2 cm internal diameter. The 
temperature was easily controlled by the vapors 
of ether, acetone, ethyl alcohol, and propy! 
alcohol at their boiling points. The capillary was 
surrounded by a long unsilvered Dewar flask 
(203"), and was calibrated by weighing the 
amounts of mercury between the marks. This 
was done by sealing the end of the capillary and 
filling to the desired mark with mercury from a 
still finer capillary. A correction was made for the 
inverted meniscus. 


DATA AND RESULTS 


We shall only give a summary of a large 
amount of data on krypton and_ hydrogen 
chloride and their mixtures. The PV—P re- 
lations are shown for two temperatures in Fig. 2 
and the needed data are collected in Table II. 


TABLE II. PV—P measurements for krypton-hydrogen 
chloride mixtures. X = mole fraction HCl. 


0 0.263 0.492 0.729 1.000 
95.5°C 0 1.350 1.350 1.350 1.350 1.350 
20 1.328 1.328 1.325 1.301 1.286 
40 1.305 1.310 1.295 1.249 1.207 
60 1.283 L2ii 1.248 1.189 1.112 
80 1.261 1.243 1.203 1.123 1.006 
100 §=1.239 =1.217 1.167 1.066 0.902 
55.5°C 0 1.203 1.203 1.203 1.203 . 1.203 
20 1.170 1.184 1.164 1.144 1.112 
40 1.337 1.144 1.115 1.073 0.997 
60 1.104 1.097 1.052 0.980 0.845 
80 1.072 1.054 0.995 0.881 0.626 
100 1.018 0.935 0.780 — 
t 
1200 
200 
2 
» 


20 P (efmospheres) 


Fic. 2, PV—P isothermals for HCl-krypton mixtures. 
(PV=1 when P=0 and T=0°C.) Data for curves 1-10 
from Table III. 
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TABLE III. Experimental values of B.Kr—HCl. 


COMPLEX FORMATION DUE TO POLARIZATION. I 


— B=0.00114(1 — X)?+-0.00102(1 — X) X +-0.00309.X? 
10 atm. 0 0.263 0.492 0.729 — 1.000 
95.5°C + 0.00111 0.00112 0.00123 0.00192 0.00310 


— B=0,00168(1 — X)?+0,00194(1 — X)X +0,00445.X? 


10atm X 0 0.263 0.492 0.729 1.000 
55.5°C —B 0.00165 0.00169 0.00196 0.00283 0.00447 


The data were calculated from the original 
observations taking PV =1 at P=0 and T=0°C. 
The PV—FP curves are given at two tempera- 
tures: 55.5° and 95.5°C. Each family of these 
curves represents mixtures of different compo- 
sitions for one temperature. It is seen that the 
second virial coefficient is a function of the 
pressure at the temperature considered, with 
exception of pure krypton where curves 1 and 6, 
Fig. 2, are straight lines over the pressure range 
studied. 


DISCUSSION OF RESULTS 


Lennard-Jones®® has shown that for mixtures 
of molecules with spherically symmetrical fields of 
force, the coefficients and bee are given by 


bep= f —1)dp 
0 


where E,g is the interaction energy between 
molecules a and 8 and from experimental data he 
was able to find the constants A, u, m and m in the 
force equation 


OE ap/dr=f(r) 


In the system HCI—Kr however the interaction 
energy is not dependent on r alone but also depends 
on the angles 0 and y because hydrogen chloride is an 
asymmetric molecule of dipole character. The 
simplified treatment given above is therefore not 
possible in the case studied. However, it still is 
possible to show interaction between hydrogen 
chloride and krypton despite the fact that the 
force constants cannot be evaluated. Although 
the second virial coefficient is a function of the 
pressure at any temperature, it is found that at a 
given pressure (10 atmospheres for example) B 
is still a quadratic function of the composition. 


*® Lennard-Jones, Proc. Roy. Soc. A106, 463 (1924). 


broken curve, Broken 
curve, B=b,,(1—x)?+ 
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Since it is of primary interest to show that by. 
does not vanish for the case of the system in 
question, the analytical evaluation of the above 
integral becomes unnecessary. For B, the second 
virial coefficient for the mixture, can be experi- 
mentally determined as a function of compo- 
sition. When the data obtained from the pressure 
volume measurements were used to determine 
the best parabola that could be fitted to the 
experimental points, it was definitely found that 
biz was different from zero in the equation of this 
parabola. It should be noted, that even if b:.=0, 
the resulting form of B is still a parabola, that is, 


B — x)? + 


and this parabola can be obtained directly from 
the pressure volume data for the two pure gases. 
It is immediately found however, that this is in 
marked disagreement with the experimental 
facts. The situation is clearly pictured in Fig. 3, 
where the experimental points are shown by the 
circles, the best parabola that can be passed 
through them by the solid line, and the simplified 
parabola, that is with b;2 set equal to zero, by the 
broken line. This last procedure amounts to the 
omission of interaction between the unlike mole- 
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Fic. 3. B—X curves for HCl-krypton mixtures. Un- 
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cules Kr and HCl, while the interaction between 
like molecules is still considered. It is seen that 
the resulting dotted curves do not fit the experi- 
mental curves nearly as well as do the parabolas 
including interaction between unlike molecules. 
The P— V—T relations studied here give then 
evidence of interaction between a rare gas atom 
and a permanent dipole molecule. The inter- 
action term },2 is seen to be of the same order of 


magnitude and of the same sign as the terms 5, 
and be. thus indicating attractive forces (Table 
III). This result is in keeping with the other lines 
of evidence obtained from band spectroscopy, 
hydrate formation, the general behavior of 
isoelectronic systems (pseudo-atoms) and their 
proton affinity. This behavior can best be 
understood on the basis of modern wave me- 
chanics as studied by London.‘ 
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Complex Formation due to Polarization. II. The System Propane and 
Hydrogen Chloride 


GEorRGE GLOCKLER, DONALD L. FULLER! AND CHARLES P. RoE, University of Minnesota 
(Received July 10, 1933) 


Molecules of saturated hydrocarbons are rare gas-like 
in structure and they may be supposed to behave in many 
respects like rare gases. They may then be expected to show 
complex formation with other molecules or atoms of dipole 
character, depending on their own polarizability (preceding 
article). An experimental study of the P-V-T relations 
of the gaseous mixtures of propane and hydrogen chloride 
shows that interaction does indeed take place between 


these two types of molecules. The expression for the second 
virial coefficient as a function of composition of the gaseous 
mixtures, is shown to require an interaction term between 
unlike molecules in order to agree with experiment. This 
question of complex formation of saturated hydrocarbons 
is further discussed on the basis of Brénsted’s theory of 
acids and bases. 


INTRODUCTION 


N a former paper? the question was discussed 

as to whether or not rare gas atoms had any 
chemical properties in the sense that they could 
react with other atoms or molecules to form 
chemical compounds. It was shown that com- 
pounds of the usual valency type were ener- 
getically unstable if they involved a rare gas 
atom as one of the reaction partners. However, 
evidence was cited to show that the rare gases 
could form loosely-bound structures with other 
molecules and that such complexes could be 
more easily formed if the rare gas in question had 
a large polarizability and if the other reaction 
partner had a high dipole moment. It is, of 
course, well known that methane for example is 
very much like a rare gas.’ It may then be 
expected that rare gas-like molecules such as 
methane should show a similar behavior to the 
rare gases in the sense that methane should show 
no further chemical reactivity by addition. That 
is, primary valency compounds of methane which 
involve an electron pair either by sharing or 
transfer are probably energetically impossible, 
and they are not known. However, secondary 


‘This article is based upon a thesis presented to the 
faculty of the Graduate School of the University of 
Minnesota by Donald L. Fuller in partial fulfillment of the 
requirements for the degree of Doctor of Philosophy. 

* See preceding article, J. Chem. Phys. 1, 703 (1933). 

*Geo. Glockler, J.A.C.S. 48, 2021 (1926). 


complexes of the hydrate type may be expected. 
More generally solvates may be formed or at 
least a further interaction between a hydrocarbon 
of large polarizability and a permanent dipole 
molecule or an excited atom may be possible. 


HyYDRATES OF HYDROCARBON MOLECULES 


Hydrocarbon molecules may then form com- 
plex secondary structures with molecules such as 
water, especially if, they themselves have large 
polarizabilities. Hydrates of the hydrocarbons 
have been described by Frocrand‘ (Table I) and 


TABLE 


CH, C;Hs CoH, 


Polarizability 2.6 4.5 6.3 4.0 3.5 
aX 
Hydrates 6H,O 7H,0 — 6H,0 


these complex secondary structures seem to be 
held together by forces due to the induced 
polarization of the hydrocarbon and the conse- 
quent interaction of the field of the induced 
dipole with that of the permanent dipole of the 
water molecule. 


ISOELECTRONIC SYSTEMS 


Just as the consideration of isoelectronic 
systems in the case of rare gases led to the view 


4M. deFrocrand, Comptes Rendus 135, 959 (1902). 
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that they may have a small residual affinity for 
protons, so again it may be supposed that the 
rare gas-like hydrocarbons should have a similar 
property. For our purposes the important ques- 
tion suggested by the theory of pseudo-atoms is: 
Does the methonium ion exist? That is, is it 
conceivable that a proton may enter the complete 
octet of methane or a similar hydrocarbon in 
accordance with the equation: 


CH,+H*-CH;*. 


At once it is realized that high polarizability of 
the hydrocarbon molecule would be advanta- 
geous, for then the forces of interaction would be 
largest. The idea that a complete structure like a 
saturated hydrocarbon can react further with 
other molecules by addition arises then in two 
ways: (1) because of its polarizability other 
molecules or ions can further interact with it; and 
(2) it may have a proton affinity like an iso- 
electronic system of similar structure. These two 
points of view remind one strongly of the 
situation encountered in discussing so-called pri- 
mary valency compounds on the basis of ordinary 
valency principles, and on the basis of Werner’s 
theory of coordination compounds. 


BRONSTED THEORY OF ACIDS AND BASES 


The question raised as to a possible proton 
affinity of saturated hydrocarbon molecules can 
be put into another language. It can be asked: 
Has the methane molecule any acidic or basic 
properties? On the principle of Brénsted,® that 
any substance able to give up a proton is of acid 
character and any substance able to take on a 
proton is of basic character, the above question is 
quite proper. On these modern views such mole- 
cules as ammonia and water are basic in nature 
when it is considered that they will add protons 
to form ammonium and hydronium ions. But 
they also ionize to a slight extent into hydrogen 
ions and they therefore exhibit also acid prop- 
erties. Ultimately it is found to be the sur- 
rounding medium which determines the pre- 
dominant acidic or basic characteristics of the 
substance in question. In this way it may be 


expected to find that methane acts like a base if 


5 J. N. Broénsted, Chem. Rev. 5, 231 (1928). 


it is brought into contact with a substance which 
gives protons readily, that is, an acid in the 
Bronsted sense. Of course the basic character of 
methane may be very weak and it may be 
difficult to discover evidence of it. However, the 
formation of hydrates noted above seems to 
indicate that even such a symmetrical molecule 
as methane shows the incipient stages of basic 
action. Substitution of other radicals for the 
hydrogen atoms of methane will no doubt have 
great influence upon the basic or acidic character 
of the resulting molecule. Such methano acids 
and bases have been considered by E. C. 
Franklin.® The basic and acidic character of the 
highly symmetrical methane molecule would be 
expected to be very slight but, to deny the 
existence of these properties would be to remove 
a very reasonable method of accounting for these 
hydrates. While the possible basic character of 
saturated hydrocarbons is the main topic for 
consideration in this paper it is of interest to 
recall that the acid character of substituted 
methanes is very marked as pointed out by 
Franklin. For example the substitution of a 
“CH” radical for three of the hydrogen atoms of 
the molecule of methane will make the hydrogen 
atom in the resulting acetylene so “‘acidic”’ that 
it is easily replaceable by metals such as silver 
and copper, and the well-known acetylides result. 
Cyclopentadiene, indene and fluorene are hydro- 
carbons which behave in the same manner. They 
have hydrogen atoms in their structure which are 
replaceable by metals.and Franklin considers 
these structures as substituted methano acids. It 
seems then possible to consider that methane is 
amphoteric in character like ammonia and water. 


NH; 


CH at 
CH,>Ht+CH;- 


While such reactions of methane probably hap- 
pen only to a very slight extent, the possibility of 
such behavior must be admitted. 


6 E. C. Franklin, J.A.C.S. 46, 2137 (1924). 
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BAND SPECTROSCOPY AND QUANTUM MECHANICS 


While evidence of a very interesting nature 
was at hand for ‘“‘rare gas mercury’”’ molecules 
obtained by Oldenberg? who showed the existence 
of band spectra carried by complex molecules, and 
while further progress is being made in their 
interpretation’ no such experimental evidence 
has as yet been procured for rare gas-like mole- 
cules such as the saturated hydrocarbons.* The 
quantum mechanical arguments of London’ on 
the other hand are entirely general, and they are 
therefore applicable to the case at hand, and on 
this basis one might again predict with confidence 
that saturated hydrocarbons can interact with 
permanently polar molecules or with excited 
atoms which have dipole moments due to their 
excitation. 


Mass-SPECTROSCOPY 


Lukanow and Schiitze” mention the possibility 
that ions of mass 17 could be hydroxyl ions 
(OH*) or methonium ions (CH;*) and they find 
ions of mass 18 (H2O* or CH,*) and mass 31 
(C;H;*). In themselves these experiments furnish 
no proof that these ions possess the proposed 
structures, but taking these suggestions together 
with similar ideas regarding complexes between 
rare gases and protons they stand in a better 
light than that of mere speculation. Such 
charged ions should as a matter of fact be really 
more stable than the neutral complexes studied in 
this research, becausé the forces acting between 
an ion and a polarizable molecule will be stronger 
at a given distance of separation than in the 
neutral case. It is also interesting to remember in 
this connection that the existence of the ion H;+ 
is fully established by thé methods of mass- 
spectroscopy" while the evidence for the ex- 
istence of H; molecules is mostly negative. The 
interaction between H.+H+—>H;+ would be 
greater than the interaction between H2+H—Hs3. 
If a neutral H; should be formed in accordance 
with Heitler-London spin-theory or ordinary 


’W. Finkeiburg, Zeits. f. Physik 81, 781 (1933). 

*Such work is being undertaken in this laboratory. 

°F. London, Zeits. f. physik. Chemie 11B, 222 (1931). 
name Lukanow and W. Schiitze, Zeits. f. Physik 82, 610 


"T. R. Hogness and E. C. Lunn, Phys. Rev. 26, 44 
(1925), 


valency consideration, then the Pauli principle 
would have to be considered, and the resulting 
odd molecule would have one of its electrons in 
some excited state. 


INTERMOLECULAR FORCES 


The various lines of evidence presented so far 
make it very plausible that a saturated hydro- 
carbon like ‘‘propane”’ and a dipole like ‘“‘hydro- 
gen chloride” should form complex structures by 
interaction. Propane was chosen for experimental 
investigation because the influence of the two 
methyl-groups was thought to have a tendency to 
make the remaining hydrogens on the central 
carbon atom more “basic” than in the case of 
methane. It was then expected that interaction 
could be demonstrated by studying the pressure- 
volume-temperature relations of mixtures of 
propane and hydrogen chloride. These experi- 
mental data will permit the determination of the 
“second-virial coefficient’”’ for the various gas 
mixtures. The variation of this coefficient with 
composition should then indicate interaction 
between the unlike molecules of the mixture.’ 
The pressure-volume relations should be given by 


PV=RT+B5BP (for one mole of gas) 
where the second virial coefficient is 
(1 —x)?+2d;2(1 —x)x + doox?. (2) 


The quantities by, and are integrals 
showing the interaction of like and unlike 
molecules 


| (3) 


Vas 


The other symbols have the following meanings: 
P=pressure, T=abs. temperature, R=gas con- 
stant, x= mole fraction of gas No. 1. 

E.s= interaction energy between molecules of 
type a and 8 where a and 6 may refer to the same 
or different species. The parabola expressed by 
Eq. (2) is to be determined experimentally for the 
propane-hydrogen chloride system, and it is to be 
decided whether or not the term in })2 is needed in 
order to make experiment fit the curve of Eq. (2). 
It is important to note that the interaction 
integral b;. need not be evaluated in order to find 
out whether or not interaction between unlike 
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molecules exist. It is seen that the integrand is 
“zero” if the interaction energy Ea, is zero. This 
situation applies even if the interaction energy is 
a function of the angle coordinates as well as of 
the distance r from the molecule. In the case of 
interaction involving a dipole like hydrogen 
chloride the interaction energy is evidently a 
function of the angles as well as of distance, but 
still it is possible to use the Eq. (2) to deduce the 
existence of interaction between propane and 
hydrogen chloride. 


APPARATUS, DATA AND RESULTS 


The apparatus used in this study of the 
pressure-volume-temperature relations of pro- 
pane-hydrogen chloride mixtures has already 
been described? and we need only show the 
results of these experiments. The PV—P re- 
lations for propane, hydrogen chloride and their 
mixtures are shown in Fig. 1 and the data are 
given in Table II. The numbers given were 


1000 
t 4250 
° 20 P (atmospheres) 


Fic. 1. PV—P isothermals for HCl-propane mixtures. 
(PV=1, when P=0 and T=0°C). Data of curves 1-9 
from Table III. 


calculated from the original observations taking 
PV=1 at P=0 and T=0°C. Two temperatures 
have been employed: 95.5°C and 116.3°C. From 
these data the second virial coefficient B has been 
obtained at 40 atmospheres pressure at the two 
temperatures for the various mixtures. The best 
parabolas which will fit the B—X experimental 
relations have been calculated by the method of 
least squares. These experimental parabolas are 
then compared with the expression (2). cas tid 


TABLE IT. PV—P measurements for propane hydrogen 
chloride mixtures. X = mol-fraction of HCl. 


\X 

T\, P 0 0.236 0.499 0.850 1.00 
116.3°C O 1.423 1.423 1.423 1.423 1.423 

20. «1.227 «1.290 1.330 1.360 

40 0.956 1.115 1.220 1.290 

60 0.880 1.084 1.202 

80 0.929 1.108 

100 1.007 
95.5°C 0 1.350 1.350 1.350 1.350 — 1.350 
20) «1.133 «1.179 =1.242 1.268 1.270 
40 0.786 0.950 1.102 1.174 1.184 
60 0.913 1.070 1.093 
80 0.688 0.955 994 
100 0.825 886 


TABLE III. Experimental values of B. 


Propane-hydrogen chloride. 


— B=0.017(1 —X)?+0.00398(1 —X)X +0.00388.X? 


40atm. YX 


0 


0.236 


0.499 


0.850 


1.000 


116.3°C —B 0.01710 0.01075 0.00618 0.00398 0.00370 


— B=0.022(1 —X)*+0.00484(1 —X)X +0.00472X? 


40atm. X 
95.5°C 


0 


0.236 
—B 0.0220 0.0138 0.0095 0.00485 0.00445 


0.499 


0.850 


1.000 


0.0IS 


° 


-B(I, 16.3°C) —> 


0.005 


0.010 


-B(I, 95.5°C) 


7 


0.8 


Fic. 2. B—X curves for HCl-propane mixtures. Ut- 
broken curves, B= buy (1—X)?+b12.X (1 Broken 
curves, 
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DISCUSSION OF RESULTS 


The comparison of the experimental B-X 
relation with the terms of Eq. (2) is shown in Fig, 
2. It is at once seen that the interaction term by. 
is needed to make experiment fit the theoretical 
parabola and that it is of the same sign and lies 
between the values of b,, and bes. Therefore these 


experiments indicate that propane molecules 
interact with hydrogen-chloride molecules, and on 
the basis of the various points discussed in the 
introduction it is believed that these results may 
be taken as further proof for the existence of 
complex polarization molecules of propane and 
hydrogen chloride. 
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Binary Systems in Two Phases. I. HCl1—Kr. II. HC1—C;H; 


GeorGE GLOCKLER, DonaLp L. FULLER AND CHARLES P. Rok, School of Chemistry, University of Minnesota 
(Received August 7, 1933) 


The P—T—X relations of the binary mixtures HCI—Kr 
and HCI—C;Hs have been determined. Retrograde con- 
densation was observed with all mixtures studied. Not 


more than one liquid phase was found for any of the 


mixtures. 
$ 


HE behavior of binary systems in two 
phases, gaseous and liquid, has been stud- 
ied by many investigators! and very interesting 
results have been obtained. For example, the 
phenomenon of retrograde condensation has been 
discovered for several cases? and while theoretical 
thermodynamic studies which made use of van 
der Waals equation have not been able to predict 
this phenomenon, still it has been concluded by 
several investigators! that these phenomena are 
not excluded by the theory. 

In order to extend our knowledge of such 
binary systems and in order to build up the 
theory pertaining to them, it seems necessary to 
study many such systems experimentally over 
the complete range of temperature, pressure and 
composition. Very few studies, complete in this 
sense, have been made* and we propose to 
investigate the mixtures of rare gases and 
saturated hydrocarbons (lower paraffins) with a 
dipole gas such as hydrogen chloride. We hap- 
pened to choose these particular cases because we 
were interested in the interaction of these 
respective molecules from the point of view of 
complex formation. Our whole research program 
will include the experimental study of any two 
component system, in the hope that ultimately a 
comprehensive theoretical treatment will be 


possible. 


tvan der Waals, Die Kontinuitdt des gasformigen und 
Jliissigen Zustandes, Vol. II, J. P. Kuenen, Verdampfung 
und Verfliissigung von Gemischen, Verlag, J. A. Barth, 
Leipzig (1881) and (1906). 

2 Cailletet, Comptes Rendus 95, 210 (1880) and Kuenen, 
reference 1. 

3F, Caubet, Zeits. f. physik. Chemie 40, 257 (1902); 
Quint Gzn. Zeits. f. physik. Chemie 39, 14 (1901). 
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APPARATUS 


The apparatus built for this study is of the 
type developed by Cailletet* but it has a few 
novel features which are of general interest. For 
instance, it is possible in our apparatus to fill the 
capillary with the various gas mixtures without 
disconnecting it from the high pressure system. 
This is a distinct advantage and speeds up 
operation considerably. Furthermore the whole 
system could be evacuated to a pressure of 10~‘ 
mm Hg by means of a high vacuum system which 
included a mercury vapor pump, McLeod gauge, 
etc. The manipulation of the apparatus need not 
be described in detail since the various operations 
necessary to obtain the data are obvious, from an 
examination of Fig. 1, to any one carrying out 
such experiments. However, a few of the devices 


Fic. 1, General plan of experimental system. 


4K, Jellinek, Lehrbuch der physikalischen Chemie, Vol. |, 
pp. 310-936, F. Enke, Stuttgart (1928). 


| 
4 
= 
SY | 
| 
th 


the 


ol. 3 


= 


BINARY SYSTEMS IN TWO PHASES 


developed in this research seem to deserve a 
short description. 


Packing for high pressure needle valve 


Considerable difficulty was experienced with 
the packing of the high pressure needle valve. We 
demanded not only that these needle valves be 
tight against high pressure but also that they be 
able to withstand evacuation while open. After 
many trials we finally obtained a Bakelite 
packing? as shown in Fig. 2a which was very 


Gloss to Metal Vecuum dont 
] 


| “Threads 


Gloss to Metal Pressure Jont Q Metal to Metal Pressure Jomt Boiling Vapor System 
ig kd fig ze fig 2t 


Fic. 2. 


satisfactory in every respect. The valve was 
mounted in a vertical position and it was so 
arranged that the mercury used as confining 
liquid would fill any dead spaces. It was possible, 
by means of a Toepler pump, to transfer a given 
quantity of gas to and from the ‘‘observation- 
capillary”’ without loss. 


Magnetic stirrer 

We noticed that the magnetic stirrers used by 
former investigators! would often show a tend- 
ency to stick just at the critical moment when it 
was desirable to stir the gaseous and liquid phases 
contained in the capillary. The shape of the 
magnetic stirrer shown in cross section in Fig. 2c 
was found to be very satisfactory and we had no 
difficulty with it. The stirrer was made of iron 
and it was not covered. By keeping a sample of 


* We wish to thank Mr. F. W. Fullerton for developing 
this packing for us. 
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iron in liquid HCI for a month, we convinced 
ourselves that liquid hydrogen chloride does not 
attack iron to any appreciable extent. The 
amount of gas (presumably hydrogen) formed 
after this treatment was negligible. A small 
solenoid used to actuate the stirrer is not shown 
in the figure. 


Other features 


Various other joints used are shown in Fig. 
2 (b, d, e) which are self-explanatory.® For 
washers we used fiber-board and leather. In 
order to stop the small magnetic stirrer from 
falling out of the capillary, we fastened a small 
platinum wire across the central hole in the lower 
washer, shown in Fig. 2b. 

The measurements were obtained by filling the 
unsilvered Dewar flask surrounding the capillary 
with water at a high temperature and letting the 
system cool. The pressure which the system 
attained was read at frequent intervals of 
temperature. The precise technique of making 
these readings is described in Caubet’s article.’ 


RESULTS 


The results of this investigation are shown 
most conveniently as P—7—X diagrams. Such 
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6 P. W. Bridgman, Physics of High Pressure, The Mac- 
millan Co., New York (1931), 
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1 
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TEMPERATURE ~ CENTIGRADE 


Fic. 4. P—T—X relations, HCIl—C;3Hs. Curve x (Mol frac- 
tion HCl). (Cf. Table II.) 


graphs are given in Figs. 3 and 4 for the systems 
HCI—Kr and HCI—C;3Hs.’ The phenomenon of 
retrograde condensation was observed in every 
case and consisted of the following behavior. As 
the pressure upon the gaseous mixture was 
increased a liquid phase made its appearance. 
Further increase in pressure caused the liquid to 
increase in amount until a maximum was 
‘reached, when still greater pressure caused the 
diminution of the liquid phase until it finally 
disappeared entirely at a higher pressure. Only 


7 The data for the system HCl—Kr are taken from a 
thesis presented to the Faculty of the Graduate School 
of the University of Minnesota by Charles P. Roe in 
partial fulfillment of the requirements of doctor of philos- 
ophy. Similarly the data for the system HCI—C;Hs are 
taken from a doctor’s thesis by Donald L. Fuller, 


one ‘‘gas” phase remained. The temperature and 
pressure ranges in which this phenomenon was 
observed are mentioned in Tables I and II. 


TABLE I. Retrograde condensation, HCl—Kr. 


Composition Temperature Pressure range 
range (°C) (atm.) 

I 1.000 

II 0.962 47.0 to 47.5 84.0 to 85.2 
III 0.879 38.5 to 39.6 87.0 to 90.0 
IV 0.752 23.0 to 28.0 88.0 to 95.4 
Vv 0.634 8.0to 16.0 87.0 to 96.0 
VI 0.489 —10.0 to —1.2 83.5 to 92.7 


TABLE II. Retrograde condensation, HCl—C;Hs. 


Composition Temperature Pressure range 
HCl range (°C) (atm.) 

I 1.000 

II 0.916 50.8 to 51.0 76.6 to 76.8 
Ill 0.846 52.8 to 53.0 74.0 to 74.5 
IV 0.759 56.0 to 56.4 71.6 to 72.6 
Vv 0.617 63.6 to 65.0 67.0 to 69.3 
VI 0.522 67.4 to 69.3 65.0 to 67.2 
VII 0.449 73.6 to 75.0 61.0 to 63.3 
VIII 0.364 78.0 to 79.2 58.0 to 60.4 
IX 0.295 82.0 to 82.7 55.0 to 57.2 
X 0.170 88.0 to 88.3 50.0 to 51.0 
XI 0.090 91.6 to 91.7 46.7 to 47.2 
XII 0.000 


Another point to note is the fact that the 
critical line (marked by Roman numerals in the 
figuies) is quite different for the two systems. 
Furthermore, only one liquid phase was ever 
observed for the two systems studied. This would 
mean that no strong repulsive forces are acting 
between the different species of molecules 
studied.® 


8 J. H. Hildebrand, Solubility, Chem.-Catalog Co., New 
York (1924). 
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The Atomic Weights of Radioactive Substances! 


Forrest WESTERN AND ARTHUR E. Ruark, University of Pittsburgh 
(Received June 16, 1933) 


(1) The isotopic weights of Pb?®, Pb?”, and Pb?, on 
the chemical scale O=16.0000, are obtained by three 
distinct methods, using in each case the assumption that 
they differ by an amount very close to unity. The results 
for Pb? are these: (1) from packing fraction of lead, 
205.96;; (2) from study of the packing fraction curve, 
205.99.; from relative abundance of isotopes in leads of 
known atomic weight, 206.00;. Mean, 205.98;+0.03. By 
taking into account a recent atomic weight determination 
on lead from Bedford cyrtolite, reasonable rounded values 
for the three isotopes are 205.98, 206.98 and 207.98 +0.03. 
(2) From these values the isotopic weights of the radio- 
active substances are obtained by adding the mass (and 
the mass-equivalent of the energy) lost in disintegration. 


The calculated values for radium and thorium, 226.02 
+0.03 and 232.03+0.03, do not agree with the Inter- 
national values 225.97 +0.01 and 232.12+0.01. The calcu- 
lated weights for U™* and U™ are 238.04+0.03 and 
234.04+0.03, respectively, while the heretofore accepted 
atomic weight of U is 238.14+0.01. In view of Aston’s 
conclusion that U does not contain more than 3 percent 
of U* or any higher isotope, it appears probable that the 
accepted value for uranium is somewhat too high. (3) 
W—N is plotted as a function of N, for the three radio- 
active series, W being the isotopic weight and N the nearest 
integer. The graphs obtained are almost linear, and their 
slopes are nearly identical. In an average way they conform 
to the equation W— N= —0.02+0.00214 (N—206). 


I. INTRODUCTORY 


F we know the atomic weight of any member 

of a radioactive series, we can obtain that of 
any other member by taking account of the 
masses of the particles emitted; and the mass 
equivalent of the energy liberated in the disinte- 
grations leading from one of these substances to 
the other. The errors involved are very small, 
now that we have precise values for the masses of 
the electron and the helium atom. Unless other- 
wise stated, we shall use the chemical scale of 
atomic weights throughout this paper. Mass- 
spectrograph results are usually expressed on the 
scale O'=16.0000. To reduce them to the 
chemical scale we subtract a correction of 2.2 
parts in 10‘, based on the relative abundances of 
the oxygen isotopes given by Mecke and Childs.” 
Further, we assume that the atomic weight of 
helium on the chemical scale = 4.00128 +0.0004; 
this is obtained from Aston’s value, 4.00216 on 
the O'* scale.3 


Computations of the atomic weights of radio- 


‘Preliminary communications may be found in Phys. 
Rev. 42, 903 (1932); 43, 205 (1933). 

*Mecke and Childs, Zeits. f. Physik 68, 362 (1931); 
a Nature 130, 22 (1932); Naudé, Phys. Rev. 36, 333 


*See also Bainbridge, Phys. Rev. 43, 103 (1933). 
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active substances were made by Stefan Meyer,* 
but his results differ considerably from those 
reported here, because he did not have available 
all the data now at our disposal and because the 
distinction between the chemical atomic weight 
scale and the O'* scale was not known. Meyer 
treats uranium as a mixture of U™* and 3 percent 
of a hypothetical isotope U?*, thereby computing 
the weight of U8 as 238.011 or 238.010. From 
this he obtains 226.09 or 226.08 for the atomic 
weight of radium, a value almost 0.1 unit higher 
than the accepted one, 225.97. We agree with his 
conclusion that the chemically-determined atomic 
weight of radium is wrong, but his value for U** 
may be questioned, being based on the idea that 
uranium contains a percentage of U*® equal to 
the actinium-uranium branching ratio, whereas 
the latter is a ratio of activities, not of masses. 
Further, Meyer computes the loss of mass 
between U?88 and Pb?®® as 32.07, in close agree- 
ment with our value. He then states that if one 
takes the atomic weight of Pb?°*+3 percent Pb?” 
to be 206.05, and adds 32.07, one arrives at 
238.12 for the atomic weight of the uranium 
complex consisting of U**8 and 3 percent of U*®. 
This lies within the range covered by Hénig- 


4 Meyer, Wien Ber. (IIA) 137, 647 (1928). 
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schmid’s chemical determinations, and Meyer 
concludes that no discrepancy exists in regard to 
the atomic weights of uranium and uranium lead. 
However, the figure 206.05 is the atomic weight 
of a lead extracted from Katanga pitchblende 
(Table I). We know now that it contains 6.7 
percent of Pb?” instead of 3 percent. Further, as 
we shall see, there is some evidence against the 
assumption that uranium contains 3 percent of 
U9, The entire subject, therefore, must be 
reconsidered. 


II. WEIGHTS OF THE LEAD ISOTOPES 


We shall derive the weights of radioactive 
substances by using the weights of lead isotopes 
as a base. The latter will be obtained in three 
ways, namely, from Aston’s packing fraction for 
lead; from the packing-fraction curve; and from 
the chemically-determined atomic weights of 
radioactive leads, together with their isotopic 
compositions. 

(1) Aston® stated that the packing-fractions of 
lead isotopes are indistinguishable from the 
value he gives for mercury, 0.810~*. Later he 
said® that the lead packing fractions would be 
difficult to measure and that they probably lie 
between zero and +1 X10~*. If we place them in 
the middle of this range we obtain 205.965, 
206.965, and 207.965+0.02 for the weights of the 
principal lead isotopes. 

(2) Let W be the weight of an isotope and NV 
the integer nearest to that weight. It is possible 
to determine the value of W—WN for the lead 
isotopes by studying the dependence of this 
quantity on atomic weight. The values of W—N 
are not available for all the individual isotopes of 
heavy elements. The numbers we shall use are 
average values for all the isotopes of an element, 
or values for individual abundant isotopes. 
Accordingly, we shall simply speak of the W—N 
value for an element. Aston’ determined the 
packing fractions of tungsten, rhenium, osmium 
and thallium by comparing their lines, or those of 
their oxides, with the mercury lines; so that any 
change in the value of W—WN for mercury 


5 Aston, Nature 123, 313 (1929). 

6 Aston, Nature 129, 649 (1932). 

7 Aston, Proc. Roy. Soc. A126, 511 (1930); 132, 487 
(1931); 134, 571 (1932). 


requires corresponding changes in the values for 
these other elements. Now, from the relative 
abundances of mercury isotopes and_ their 
packing fractions we find the atomic weight 
200.602+0.05, while the International atomic 
weight is 200.61+0.004.° Thus it seems reason- 
able to increase Aston’s isotopic weights for 
these heavy elements by half the difference 
between these numbers, or 0.004. Then for the 
isotope Hg, W—WN is —0.024. If we plot the 
resulting W— WN values against N, it is found that 
they do not lie even approximately on a straight 
line, and it would be difficult to determine W— N 
for lead by extrapolating the general trend. 
Accordingly, we adopt another procedure, as 
follows. In Seetion III, we shall obtain the 
differences of the values of W—N for radioactive 
elements from disintegration data. If we arbi- 
trarily assign a (W—WN) value to the lead 
isotopes, higher than the correct one by an 
unknown amount, c, a knowledge of these 
differences makes it possible to plot W— N-+< as 
a function of NV. The plotted points can be shifted 
up and down as a body, until a straight line 
representing their average trend passes through 
the plotted point for Hg?*®. This line crosses the 
abscissa 206 at an ordinate of —0.008, so we 
take 205.992+0.03 as a value of the isotopic 
weight of Pb?°®. These relations are shown by the 
line marked XX’ in Fig. 1. However, all the other 
data in this figure are plotted on the basis of the 
value 205.98 for the weight of Pb?°®. This is a 
mean value obtained from all three methods 
outlined in this section. Accordingly, the W—N 
values plotted for tungsten, rhenium, osmium, 
mercury and thallium are 0.004—0.012 or 0.008 
lower than those obtained directly from packing 
fraction data, above. 

(3) Our third method makes use of the relative 
abundances of the isotopes in ordinary and 
radioactive leads,® together with the atomic 
weights of these leads, determined chemically. 
Data are given in Table I. The assumption is 
made that the weights of these three isotopes 
differ by integers on the chemical scale, so that 


8 Hénigschmid, Birckenbach and Steinheil, Ber. 56B, 
1212 and 1219 (1923). 

9 Aston, Nature 120, 224 (1927); 123, 313 (1929); 129, 
649 (1932), 
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Fic. 1. Deviation of atomic weights of heavy elements 
from the nearest integer. Chemical scale, O = 16.0000. 


they may be written 206+, 207+, and 208+., 
respectively.'° We may use Wilberforce uraninite 
as an illustration. From Table I, 206.195 
= 0.859(206-+x) +0.83(207 +x) +0.058(208+-). 
Then x= —0.004; that is, the calculated weight 
of Pb*°* is 205.996. 

On the basis of this assumption the four 
samples of lead in Table I yield the following 
values for the weight of Pb: 


Material Isotopic weight of Pb6 
Ordinary lead 206.030 
Katanga pitchblende 205.981 
ilberforce uraninite 205.996 
Norwegian thorite 206.005 
Average 206.003 +0.035 


Weighting all three methods equally, the final 


Study of data on lighter elements indicates the degree 
to which this assumption may be trusted; the (W—N) 
values determined by the mass-spectrograph for the iso- 
topes of a single element never differ by more than 0.005 
of a unit, and usually the departure is only 0.001 or 0.002. 
In the radioactive domain, however, the differences are 
larger, as we can see from Fig. 1. It appears likely that 
our assumption is correct within 0.01 of a unit of atomic 
weight. It should be considered merely as a convenient 
working hypothesis, 


TABLE I. Isotopic constitution of ordinary and 


radiogenic leads. 

Ordi- Katanga  Wilber- 

nary pitch- force Norwegian 
Source of lead: _ lead blende  uraninite thorite 

Chemical 

Atomic weight: 207.22" 206.048» 206.195° 207.904 
Percentage 206: 27.75 93.3 85.9 4.6 
Percentage 207: 20.20 6.7 8.3 Ls 
Percentage 208: 49.55 (0.02) 5.8 94.1 


Other isotopes: 2.50 


® International value. 

» Hénigschmid and Birckenbach, Ber. 56, 1837 (1923). 

¢ Baxter and Bliss, J. Am. Chem. Soc. 52, 4851 (1930). 

4 Fajans, Sitz. Heidelb. Akad. Wiss. 3 (1918); Hénig- 
schmid, Phys. Zeits. 19, 437 (1918); Zeits. f. Elektrochemie 
25, 91 (1919). 


mean for Pb? is 205.987 +0.03. Simultaneously 
with our announcement of the above results, a 
paper by von Grosse! appeared, in which he 
gives 205.96+0.02 as the weight of Pb?°°. This is 
obtained from the data of Table I; the method of 
computation is not explained. We believe our 
mean is better supported by the data above, than 
the lower value of von Grosse. However, recent 
atomic weight determinations by Baxter and 
Alter” yield a lower value for Pb?°* than any so 
far obtained. They have determined the atomic 
weight of lead extracted from a thorium-free 
cyrtolite occurring at Bedford, N. Y., with the 
result 205.924+0.02. An isotopic analysis of the 
Bedford material is much to be desired, for the 
discrepancy between this value and our mean 
may be due to the presence of the isotopes 203, 
204, or 205, found by Aston." Meanwhile, we 
shall give the cyrtolite value some weight by 
adopting the rounded value 205.98 for Pb?®°. 
This mean depends on the results of three 
independent methods of attack, employing a 
number of atomic weight determinations, in- 
cluding Baxter and Bliss’ work on Wilberforce 
lead, but it is quite possible that a lower value 
may be confirmed by later developments. If so, 
the effect will be to bring our calculated value for 
the atomic weight of radium (see Table II) into 
better agreement with the observed one; but the 
calculated and observed values for uranium will 


still disagree. 


1 yon Grosse, Phys. Rev. 42, 565 (1932). 
2 Baxter and Alter, Science 76, 524 (1932). 
13 Aston, Nature 129, 649 (1932). 
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III. COMPUTATION OF THE ISOTOPIC WEIGHTS OF 
RADIOACTIVE SUBSTANCES 


We shall now consider mass losses in disinte- 
gration. When an alpha-particle is emitted, two 
planetary electrons must leave the atom, so the 
atomic weight of a neutral helium atom plus the 
weight equivalent of the energy loss must be 
subtracted from the atomic weight of the parent 
atom to obtain that of the daughter. However, 
when a beta-particle is emitted, the planetary 
shells of the daughter atom pick up an electron. 
The only mass lost is that associated with the 
kinetic energy of the beta-particle and the 
energy of the gamma-rays. 

In computing the energy loss in beta-ray 
transformations we use the average energy of the 
disintegration electrons rather than the maxi- 
mum. It is well known that the average energy 
obtained with the beta-ray spectrograph is in 
agreement with the energy loss determined 
calorimetrically. The experimental facts require 
that we adopt one of several suppositions: 

Either it is true that (1) energy is conserved 
only statistically in beta-ray transformations; or 
(2) an atom which sends out a slow disintegration 
electron must rid itself of surplus energy by the 
emission of penetrating rays or particles which 
are not detected in the calorimetric measure- 
ments or indeed in any others so far made; or (3) 
nuclei of a given species do not all have the same 
energy. Several other hypotheses, more or less 
complicated, might be mentioned. Our method of 
procedure corresponds to the first assumption, 
but if we were to use the maximum energy of the 
disintegration electrons, none of the atomic 
weights we calculate would be altered by more 
than a few thousandths. This is negligible: in 
comparison with the uncertainty of 0.03 involved 
in the weights of the lead isotopes. 

In the case of isotopes whose gamma-ray 
spectra are known, we have used the energy loss 
corresponding to the shortest gamma-ray re- 
corded, in the hope that its energy is fairly close 
to the average surplus possessed by the excited 
nucleus just after disintegration. A detailed 
understanding of nuclear energy levels is not 
necessary for our purpose because the loss of 
energy by radiation corresponds in most cases to 
a negligible mass-loss. Where the spectrum is not 
known, we neglect the loss of energy by radiation. 


Inspection of known gamma-ray spectra shows 
that this procedure cannot possibly introduce an 
error greater than about 0.005 in any isotopic 
weight. Probably the average error thus intro- 
duced is closer to 0.001 unit. 


The factors used in computation are these: ‘ 


Mass energy of atom having unit chemical atomic 
weight, 1.482 10-* ergs, or 932 10° electron-volts. 
One electron-volt corresponds to 1.074 10-9 
chemical atomic weight units. Energies of alpha- 
particles, recoil atoms and beta-particles are 
taken from Gamow’s" Atomic Nuclei and Radio- 
activity, with few exceptions. The data for U™’, 
U**, and Th are those of Kurie, and of Kurie and 
Knopf,'® while the values for U?*® and U*»® are 
assumed to be the same as those for U*** and 
U**4, respectively. For the hypothetical beta- 
rayers between U*® and U***, we assume the 
disintegration electrons have the same energies as 
those of UX; and UXs, respectively. A similar 
assumption is used in dealing with MsTh, and 
MsThg. No information is available concerning 
the gamma-rays of UY, Ac, An and AcA. We 
assume they are negligible. 


Most of the gamma-ray energies were taken _ 
from Radiations from Radioactive Substances, by ie 


Rutherford, Chadwick and Ellis.'® In a few cases | 


it was necessary to estimate them from ab- 


sorption coefficients given on page 456 of the | 


same work, by use of the Klein-Nishina formula. 
The uncertainties thus introduced are negligible 
for our present purpose. 

The results are shown in Table II.’ The 
isotopic weights are given to a higher accuracy 
than the data warrant, to avoid error in rounding 


off. The outstanding facts which emerge from a ~ 


study of these numbers are: (1) that in a gross 


4 Gamow, Atomic Nuclei and Radioactivity, pages 35 to 
54, Oxford University Press. 

1 Kurie, Phys. Rev. 41, 701 (1932); Kurie and Knop, 
Phys. Rev. 43, 311 (1933). 


16 Rutherford, Chadwick and Ellis, Radiations from | 


Rad: Substances, pp. 360-377, Cambridge University 


Press, 1930. 


7 Phys. Rev. 43, 38-59 (1933) contains a group of | 


papers by F. Allison and his associates, in which large 
numbers of new radioactive isotopes are proposed. These 
are fitted by Miss Bishop into four radioactive series. 
Independently Piggot proposes a somewhat similar ar 
rangement. These schemes are highly hypothetical. We 
base our discussion on the disintegration sequences ge” 
erally accepted. 
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TABLE II. Mass losses and isotopic weights for the radioactive elements. 
(Chemical scale.) 


Mass equivalent of energy loss 


Alpha-particle Beta- Gamma- Total mass Isotopic 

Isotope and recoil atom particle ray loss weight 
a 0.0044 0.0044 +He 238.04 

UX: 0.0001 .0001 234.04 

UX. .0004 0.00015 .00055 234.04 

Us .0050 .0050 + He 234.04 

lo .0050 .0050 + He 230.03 

Ra .0052 .0052 +He 226.02 

Rn 0059 .0059 + He 222.02 

RaA .0065 .0065 + He 218.01 

Ra B .00025 .0005 .00075 214.00 

Ra C .0008 .0024 0032 214.00 

Ra C’ 0084 .0084 +He 214.00 

Ra D 209.99 

Ra E .0004 .0004 209.99 

Ra F .0058 .0058 + He 209.99 

RaG 205.98 

.004(?) .004 + He 239.05 
Beta-Rayer? .0001(?) .0001? 235.04 
Beta-Rayer? .0004(?) .0004? 235.04 

Us .005(?) .005 +He 235.04 
231.03 
Pa 0055 .0003 .0058 + He 231.03 
Ac 227.03 
RdAc .0065 .0007 .0072 +He 227.03 
AcX .0062 .0062 + He 223.02 
An .0075 .0075 +He 219.01 
AcA .0081 .0081 +He 215.00 

Ac B .0002 .0005 .0007 210.99 

Ac C .0005 .0077 + He 210.99 
mec” .0004 .0022 .0026 206.98 

Ac D 206.98 

Th .0040 .0040 + He 232.03 
MsTh, .0001? .0001 228.03 
MsTh» 0004? 0021 0025 228.02 

RdTh .0058 .0001 .0059 +He 228.02 

ThX .0062 .0062 +He 224.02 

Tn .0068 .0068 + He 220.00 

ThA .0073 .0073 -+He 216.00 
ThB .0001 .0003 .0004 211.99 
The .0066 .0066 +He 211.99 
1 Th C” .0007 .0027 .0034 207.98 


207.98 


way, W—N is a linear fraction of atomic weight 
for each radioactive series, and (2) that the mean 
rate of increase of W—N with atomic weight is 
nearly the same for the three series. The equation 
W—N=-—0.02+0.00214(N — 206) is a fair aver- 
age representation of the atomic weights in the 
three series. 

Isotopic weights of UI, Ra, Rn, and Ra F have 
also been computed by von Grosse,'* whose 
results are 0.02 lower than ours. 

The chemically-determined atomic weight of 
radium is 0.05 unit lower than our calculated 
value, and about 0.15 unit lower than the value 


‘von Grosse, Phys. Rev. 42, 565 (1932). 


obtained by calculation, on any reasonable basis, 
from the chemically-determined atomic weight 
of uranium. For thorium the discrepancy is 0.09 
unit. In both cases the discrepancies are larger 
than the sum of the estimated uncertainties of 
the calculated and experimental values. 

The value determined chemically for uranium 
is higher than one would expect on any as- 
sumption which appears probable in view of 
known facts. It can be understood if sufficient 
amounts of isotopes heavier than U** are 
present, but Aston!® found only the isotope U”*, 
and says, “uranium is probably simple to at 


19 Aston, Nature 128, 725 (1931). 
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least two or three percent.’’ Wahlin” also found 
only a single isotope. It is possible, however, that 
several percent of U*® would have escaped 
detection in his experiment. If uranium were 
composed of the isotopes 238, 234 (in negligible 
amount), and 3 percent of 239, the atomic 
weight of the mixture calculated from our 
isotopic weights would be 238.07, which is 0.06 
lower than the chemical value. Of course, small 
percentages of several higher isotopes would 
suffice to bring about agreement, but there is 
no positive evidence for their existence. As to 
U9, we have summarized”! the evidence on the 
origin of the actinium series, concluding that it is 
quite improbable that this isotope exists. 


IV. CoNCLUSIONS 


Our conclusion is that the calculated values for 
the atomic weights of radium and thorium 


20 Wahlin, Phys. Rev. 39, 183 (1932), and private com- 
munication. 
21 Western and Ruark, Phys. Rev. 43, 205 (1933). 


should be used in preference to those determined 
chemically, and that the accepted value for 
uranium appears too high on the basis of present 
information as to the isotopes of this element. In 
making this judgment, we wish to express 
admiration for the work of Hénigschmid and his 
collaborators on these atomic weights. In the case 
of uranium, for example, Hénigschmid and 
Schilz?? made 37 determinations by several 
methods. For all three elements the inner con- 
sistency of the results leaves little to be desired, 
but it appears that systematic errors beyond the 
control of the chemist must have played a part. 
The only other possibilities are to assume that 
Einstein’s energy-mass relation is not strictly 
valid or that energy is not even Statistically 
conserved, in radioactive disintegration. Against 
both of these assumptions, a wealth of experi- 
mental evidence can be adduced. 


2 Hénigschmid and Schilz, Z. anorg. chem. 170, 145 
(1928). 
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An Elementary Deduction of Gibbs’ Adsorption Theorem 


F. O. Koenic R, C. Swain,* Stanford University 
(Received August 3, 1933) 


The Gibbsian treatment of the thermodynamics of 
surfaces is very long and abstract. On the other hand the 
equations designated as ‘‘Gibbs’ adsorption theorem’”’ in 
all texts that have come to the authors’ attention are only 
special cases, incompletely or inexactly derived, of the 
Gibbsian equation for constant temperature, 


dpi. 


An exact elementary deduction of this equation is there- 
fore desirable. Such a deduction is given by means of a 
thermodynamic engine. Various closely related equations, 
all for constant temperature, are also deduced. The relation 
of the equations obtained to the most general Gibbsian 
equations, which include variations of temperature, is indi- 
cated. The assumptions underlying the Gibbsian treatment 
of the thermodynamics of surfaces are discussed briefly. 


A. INTRODUCTION 


HE term “‘Gibbs’ adsorption theorem” has 

been rather indiscriminately applied to 
several more or less exact thermodynamic equa- 
tions relating the surface tension and the adsorp- 
tion at a surface layer to the concentrations, 
vapor pressures or chemical potentials of the 
components in the adjacent phase interiors. 
These equations may all be derived from the 
following general equation first given by Gibbs:! 


K 
do = —s*dT — >> (1) 
i=1 


o is the surface tension, J the temperature, K 
the number of components of the system, yu; the 
chemical potential (per mole) of the ith com- 
ponent, I’; is usually referred to either as the 
amount (in moles) of the ith component adsorbed 
per unit area of surface layer or the surface 
density of the ith component, and similarly s¢ 
is called the entropy per unit area of surface 
layer or the surface density of the entropy; the 
exact definition of the quantities T; and s® is 
given below. For variations at constant tem- 
perature Eq. (1) becomes 


K 
do=— T dpi. (2) 


* National Research Fellow, 1933. 

‘J. W. Gibbs, Collected Works, Vol. I, p. 230, Eq. 508 
(Longmans Green, 1928). The notation in the present 
Paper differs somewhat from that of Gibbs. 


It is this equation which in the present paper will 
be referred to as Gibbs’ adsorption theorem. 
Gibbs’ deduction of Eq. (1) is very general and 
at the same time very exact. Another deduction 
both general and exact has been given by Defay.? 
Because of their completeness these treatments 
are so long and abstract that they have not been 
included in the standard texts on thermody- 
namics and on physical and colloid chemistry. 
The equations given in all texts that have come 
to the authors’ attention are special cases of Eq. 
(2), incompletely or inexactly derived. It is the 
purpose of this paper to give an exact elementary 
deduction of Eq. (2) and some related equations. 
The deduction is carried out with the aid of the 
thermodynamic engine described below. 


B. THE ENGINE 


Let a system composed of two phases a and 8 
separated by a plane surface layer w (henceforth 
designated as the system aw) be contained in 
the thermodynamic engine indicated in Fig. 1. 
The surface layer w, which is in general an in- 
homogeneous region of finite thickness, is 
represented by the shaded portion in Fig. 1. 
C is a geometrical surface whose significance is 
explained below. The material of the engine is 
supposed to be absolutely inert so that the 
properties of its contents remain unchanged 


2R. Defay, Etude Thermodynamique de la Tension Super- 
fictelle, etc., 1°* Partie, reprinted from Bull. Acad. roy. 
Belg. (Classe des Sciences) 1929, 1930, 1931 (Hayez, 
Brussels 1931). 
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Fic. 1. Engine for demonstrating the thermodynamic laws 
of surface layers. 


right up to the engine walls. B, B®, BY are 


frictionless pistons. The cylinders qi%, 


attached to the part of the engine containing the 
phase a contain the separate components of a 
in the pure state. Each pure component in its 
cylinder is separated from a by a membrane 
permeable only to that component. The part of 
the engine containing the phase @ is similarly 
equipped with cylinders qi°, containing 
the components of 8 in the pure state, separated 
from 8 by semipermeable membranes. If all K 
components of the system afw are appreciably 
soluble in both phases, then there are K cylinders 
attached to each phase and the total number of 
the cylinders g,*, g° is 2K. But if some of the 
components are appreciably soluble in only one 
phase, the corresponding cylinders are omitted 
from the other phase, and the total number of 
the cylinders q:*, g is less than 2K by the 
number of components restricted to one phase 
only. By means of the frictionless pistons };%, 
be*+ + -bx* and the pure components 
may be reversibly added to or removed from the 


phases a@ and 8 respectively. The components in 
the cylinders g;* must therefore be in chemical 
equilibrium with the phase a and the components 
in g#® with the phase 8, i.e., the pressures in q;* 
and q,6 must be such that 


and (3.1), (3.2) 


where y,°" and yu, denote the chemical potentials 
of the pure ith component in the cylinders q;* 
and respectively, and and yu, the chemical 
potentials in the phases a and 8 respectively. 
Finally it may be noted that the special case in 
which some of the components are practically 
insoluble in a and 6 but present in the surface 
layer (‘‘insoluble surface film’’) may be treated 
by attaching cylinders g;“ (not shown in Fig. 1) 
for, these components to the engine at such a 
level that their semipermeable membranes are 
in contact with w. The condition for the chemical 
equilibrium between the pure component and the 
surface layer is then analogous to (3.1), (3.2), 
namely 

Mi (4) 


C. THE DEFINITION OF THE EXTENSIVE PRop- 
ERTIES OF THE SURFACE LAYER 


Let C (Fig. 1) be a geometrical surface placed 
parallel to the surface layer at an arbitrary level 
with respect to the latter. C will henceforth be 
designated as the Gibbs surface. Let the two 
volumes into which the system afw is divided by 


C be denoted by V and V? respectively. Let G 


denote the value of any extensive property such 
as energy, entropy or the amount of a component 
for the system afw, and let g* and g® denote the 
value of the property in question per unit 
volume of the homogeneous phase interiors « 
and ~. Then the quantity G* defined by the 
equations. 


(5) 


G8 = yb (6.1), (6.2) 


is to be regarded as the energy, entropy, etc., of 
the surface layer. G* evidently represents the 
difference between the actual energy, entropy, 
etc., of the system afw and that which it would 
have if a and 8 both remained homogeneous 
right up to the Gibbs surface C. 

The extensive properties G* are evidently 
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proportional to the area of the Gibbs surface, 
which will be denoted by @. It is therefore 
expedient to define surface densities G*/Q which 
correspond formally to the volume densities 
g*, g? in the phase interiors. Thus the surface 
densities of the entropy and of the ith com- 
ponent, appearing in Eq. (1) and designated by 
s* and I; respectively, are defined by 


s* = S“/Q, 


where S* is the surface entropy and n,° is the 
number of moles of the ith component adsorbed 
at the surface. 

Since the position of the Gibbs surface C is 
arbitrary, the quantities G* and their surface 
densities G*/Q are not completely determined for 
a given state of the system a$w and a given value 
of 2, but are functions of the position of C 
relative to the surface layer w. The only G* and 
I’; values of physical interest, however, are those 
for which the Gibbs surface is placed within or 
very near the inhomogeneous region w. Mathe- 
matically the position of the Gibbs surface 
relative to the surface layer may be fixed for 
each state of the system by adopting certain 
(arbitrary) conventions with respect to the 
values of ';. A number of useful conventions of 
this kind have recently been proposed and the 
relations between them discussed in detail by 
Guggenheim and Adam.* Throughout the fol- 
lowing discussion the level of the Gibbs surface and 
the level of the surface layer are taken to mean the 
respective levels relative to the engine wall. 

Finally it is to be noted that since according 
to the above definitions the volume V of the 
system afw is given by 


V=Ve4 (8) 


the extensive quantity V* is necessarily zero and 
thus occupies a special position among the G*. 


D. THe AssuMPTIONS 


The deductions given in this paper depend 
upon elementary principles and upon the five 
assumptions stated below. Of these assumptions 
only the last three are necessary for the validity 
of the equations deduced. The first two are made 


*E. A. Guggenheim and N. K. Adam, Proc. Roy. Soc. 


(7.1), (7.2) 
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merely in order to simplify the working of the 
engine sufficiently; they thus indicate the extent 
to which the present deduction of Eq. (2) is less 
general than that of Gibbs. 


1. The composition of the two phase interiors a 
and 6 and of the surface layer w can be com- 
pletely described in terms of the amounts 
(in moles) n;*, n#, n;* of K electrically 
neutral components 

By the general method of Gibbs it can be 
shown that Eq. (2) holds also when the com- 
position of the system must be described in 
terms of electrically charged components, if for 
each such component the electrochemical poten- 
tial is introduced in the same way as the 
chemical potential yu; of an electrically neutral 
component. It may be noted that because of the 
electroneutrality of all phase interiors (at ordi- 

nary temperatures and pressures), assumption 1 

is always fulfilled whenever the state of the 

surface layer is completely determined by the 
state of the two adjacent phase interiors. The 
simplest case in which it is not fulfilled is that 
in which one of the two phases is a conductor, 
the other a dielectric (e.g., NaCl solution | air) 

and the surface layer between them carries a 

net electric charge. Another case in which it is 

not fulfilled is that of the systems giving rise to 
electrocapillary curves elsewhere discussed by 
one of the authors.5 


2. The surface layer is plane 


By the general method of Gibbs it can be 
shown that Eq. (2) applies to curved as well as 
to plane surface layers under two conditions: 
(i) the radii of curvature of the surface layer are 
very large compared with its thickness; (ii) the 
Gibbs surface (whose position relative to a plane 
surface layer is entirely arbitrary) is placed within 
or very near the (curved) surface layer. There 
are no particular disadvantages connected with 
the latter restriction because, as already men- 


4On the concept of the electrochemical potential see 
E. A. Guggenheim, J. Phys. Chem. 33, 842 (1929); 34, 
1540 (1930); Modern Thermodynamics by the Methods of 
Willard Gibbs, in print (Methuen, London, 1933). 

5F, O. Koenig, J. Phys. Chem., in print; Wien-Harms, 
Handbuch der Experimentalphysik., Vol. XII, Part 2, in 
print (Akademische Verlagsgesellschaft m.b.H. Leipzig). 
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tioned, the only G* and IT; values of physical 
interest are those for which the Gibbs sur- 
face is placed within or very near the surface 
layer anyway. If the radii of the surface layer 
become of the same order as its thickness, then 
the whole problem becomes more complicated 
and Eq. (2) is no longer applicable. 


3. Macroscopic electric, magnetic and gravita- 
tional fields are absent or of negligible 
influence 

In the absence of gravity the pressure and 
composition of any phase interior are constant 
throughout all its parts. In a gravitational field 
the pressure and composition are functions of 
the gravitation potential, i.e., they are variable 
in the direction of the field. In the earth’s field 
however the effect of this variation on the values 
of surface tension chemical potential and 
adsorption [; in the systems usually studied 
is entirely negligible. Consequently any theo- 

retical result connecting the quantities o, w:, T; 

(Eq. (2)) and obtained on the assumption that 

gravity is absent is also valid in gravitational 

fields of terrestrial magnitude. 


4. The first and second laws of thermodynamics 
These laws are used in the form of the theorem 
that in any isothermal process the gain in the 
Helmholtz free energy of a system is equal to the 
total reversible work done on the system in the 
process. They are also implied in the equilibrium 
conditions stated below (assumption 5). 


5. The system afw is in thermal, mechanical and 
chemical equilibrium 

The condition for thermal equilibrium is that 
the temperature TJ be the same throughout the 
system. 

The condition for the mechanical equilibrium 
of two phases separated by a plane boundary 
layer (assumption 2) in the absence of gravity 
(assumption 3) is that the pressures P* and P® 
in the two phase interiors be equal, i.e., 


Pi =P. (9) 


The condition for chemical equilibrium is that 
for each electrically neutral component (assump- 
tion 1) appreciably soluble in both phases, the 
chemical potentials and in the two phase 


- amounts, dV*, dV®, dQ, simultaneously moving 


interiors be equal, i.e., 
wit = = Mi. (10) 


From this equation and (3.1), (3.2) it follows 
that for any component soluble in both a@ and 8, 
the chemical potentials in the corresponding 
cylinders q:%, g:° are also equal, i.e., 


08 


and hence that the pressures ;*, p,° in the 
cylinders q;*, gi° are equal, i.e., 


= pF = pi. (12) 


For any component, say the Ath, soluble only in 
a and not in 6 or vice versa, the upper suffix a 
or 6 may without loss of generality be omitted 
in all following formulas and the notation of 
Eqs. (11), (12) preserved by writing: 


case in which some of the components are 
soluble only in the surface layer (insoluble surface 
film, Eq. (4)), so that for any such component, 
say the jth, 

(15) 


= pi. (16) 


E. Tue DeEpuctTIon oF EQ. (2) 


Eq. (2) may be deduced by carrying out the 
following two processes with the engine, equating 
the gain in the Helmholtz free energy (abbrevi- 
ated Hfe.) of the entire engine (i.e., of the system 
oBw+the cylinders q:*, qi (and g,*)) to the total 
work done on the entire engine in each process 
and combining the resulting equations. 


Process 1 

By moving the pistons B*, B*, Be and the 
level of the Gibbs surface C infinitesimally, 
change V2, V8 and Q by arbitrary infinitesima 


the pistons 6,“, b,8 (and 6, in the case of insoluble 
surface films so as to change the number 0 


moles 1; of the components in afw by arbitrary 4 


infinitesimal amounts dn,. 
The gain in the Hfe. of the entire engine ! 


(13.1), (13.2) | 
or (14.1), (14.2) | 


The same convention is applicable also to the | 
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this process may be calculated as follows. From 
Eq. (5) it follows that the gain dA in the Hfe. of 
the system a®w in any infinitesimal process 
whatever is given by 


dA =dA*+dA’+dA*. (17) 


The gain in the Hfe. of the cylinders q,*, gf 
(and g;*) is evidently equal to 


K 

(18) 

i=1 

where A;° is the molal Hfe. of the 7th pure com- 
ponent at 7, p; and dn; represents the amount 
of the component transferred from both cylinders 
git and qg (or from q;*) to the system aBw. The 
gain in the Hfe. of the entire engine is therefore 


K 
dA*+dA8+dA*—> (19) 
i=1 


The work done by the entire engine consists of 
the volume work done by B+, B®, Be: 


—PdV*—Pdvé (20) 
the surface work done by B*: 
odQ (21) 
and the volume work done by 6;*, 6° (0,°): 
K 
(22) 


where v,° denotes the molal volume of the pure 
ith component at 7, p;. Hence for process 1 


K 
dA*+dA8+dA° Avdn; 


i=1 


K 
= —PdV*—PdV8+od2+> pwidn; (23) 


i=1 


dA*+dA8+dA* = — 
K 
(24) 
i=1 


A? is related to the Gibbs free energy F;° per 
mole of pure component by 


(25) 


and F,’ is in turn equal to the chemical potential 
ue. Hence by Eqs. (11), (13), (15), 


Ai+ pve =p; (26) 


so that Eg. (24) becomes 
dA*+dA'+dA* = 


+E uidn;. (27) 
Since by Eq. (5) 
(28) 
Eq. (27) becomes 
dA*+dA*+dA* = —PdV*—PdV*+odQ 


K K K 
+2 midniet+ widn P+ (29) 
The state of the system afw is completely 
determined by the temperature 7, the volume 
V(= V*+ V*), the area 2 and the number of 
moles n;(=n;*+n+n,") of the components. In 
process 1 all of these quantities except T have 
been varied in an arbitrary way. Hence Eq. 
(29) applies to any isothermal change whatever 
to which the system afw in equilibrium can be 
subjected. 


Process 2 

Imagine the piston B* to have been pushed 
as far to the left (Fig. 1) as possible, so that 2 
is zero and the phases a and £8 are isolated 
mechanically and chemically. Then by moving 
the pistons B* and b,* change V* and n,“ by 
arbitrary infinitesimal amounts dV@ and dn,*. 

The gain in the Hfe. of the entire engine in 
this process is 


Ai*dn;* (30) 
and the work done on the entire engine is 
—Pd pwviidn;*. (31) 
Hence 


K K 
dA*—¥> afdn;*= —PdV*+> (32) 
int 


i=1 


K 
dA*=—Pd ve+>. (33) 


i=1 


and, by Eq. (26), 
K 
dA*= —PdV*+> uidn;*. (34.1) 
i=1 
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The state of any phase a is completely deter- 
mined by the temperature 7, the volume V2 
and the number of moles ;* of the components. 
In process 2 all of these quantities except T 
have been varied in an arbitrary way. Hence 
Eq. (34.1) applies to any isothermal change 
whatever to which the phase a in equilibrium 
can be subjected. The corresponding equation 
for the phase 8 is evidently 


K 
dA® = —PAVP+ pidné. (34.2) 
i=1 


Combination of the equations resulting from 
processes 1 and 2 
Since Eqs. (29), (34.1), (34.2) are perfectly 
general for any infinitesimal isothermal change 
in which equilibrium is preserved, the subtraction 
of (34.1) and (34.2) (derived by process 2) from 
(29) (derived by process 1) leaves 


K 
dA*=od2+> (35) 


i=1 


and this equation must also hold for any in- 
finitesimal isothermal change to which the 
system afw in equilibrium can be subjected. 

For a given state of the system afw the in- 
tensive quantities in Eq. (35) are fixed, but Q 
is still variable, and the extensive quantities 
A*, n“ are proportional to © for a given position 
of the Gibbs surface relative to the surface layer. 
Consequently Eq. (35) may be integrated, keep- 
ing the intensive quantities and the position 
of the Gibbs surface relative to the surface 
layer constant, and allowing @ to vary, the 
result is 


K 
(36) 
inl 


where the integration constant ® is a function 
of T, o, uw; and the position of the Gibbs surface 
relative to the surface layer. When is zero, A* 
and n,* are also zero and therefore ® is zero, 
leaving as the complete expression for the Hfe. 
of the surface layer 


K 
une. (37) 


i=1 


The differentiation of (37) gives 


K K 
dA* =odQ+ n@dpj. (38) 

i=1 i=1 
This equation applies to any infinitesimal 
change whatever, and therefore in particular to 
the general isothermal change to which Eq. (35) 
refers. Consequently the subtraction of (35) 
from (38) gives 


K 
0=Qde+)Dd (39) 
i=1 
which, like (35), holds for any infinitesimal 
isothermal change in which equilibrium is 
preserved. Solving (39) for do and _ taking 
account of Eq. (7.2) gives Gibbs’ adsorption 
theorem, Eq. (2): 


K 
do=— Tidus. (2) 


i=l 


F. SoME FURTHER EQUATIONS DIRECTLY 
DERIVABLE FROM THE ENGINE 


The important relation (35) was obtained by 
a mathematical operation carried out upon 
Eqs. (29), (34.1), (34.2) which were directly 
derived from the engine. It is perhaps worth 
while indicating how some of the consequences 
of (35) can be obtained directly from the engine. 

Thus it is possible to examine the effect of a 
change in the relative positions of the Gibbs 
surface and the surface layer, the state of the 
system and Q remaining constant. For this 
change (35) gives 


K 
dA* =>) (40) 
i=1 
which, since the uw; are constant, may be in- 
tegrated, giving for a finite shift of the Gibbs 
surface relative to the surface layer 


K 
piAn,’. (41) 

i=1 
Each of the Eqs. (40), (41) is directly derivable 
from the engine. To derive (40) it is only neces- 
sary to shift the level of the surface layer 
infinitesimally with respect to the engine wall, 
the level of the Gibbs surface with respect to the 
engine wall, V*, V8, 2 and all intensive variables 
being kept constant. This can be done by 
keeping B*, B®, Be stationary and operating 
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b,*, b# so as to add to one of the phases, say a, 
its components in the proportions in which they 
are present in the interior of a, simultaneously 
removing from 8 its components in the propor- 
tions in which they are present in the interior 
of 8, the amounts added to a and removed from 
8 being such that P remains constant. From 
Eqs. (5), (6.1), (6.2) it follows that for any 
infinitesimal change whatever 


dG =d(g*V2) +d(g6 V8) +dG*. (42) 


Since in the change in question V*, V® and the 
intensive variables g*, g° remain constant, 


dG =dG* (43) 

so that in particular 
dA =dA* (44) 
(45) 


The gain in the Hfe. of the entire engine is 
therefore 


dA®— (46) 


K 
=1 


and the work done on the entire engine is 


K 
> (47) 


i=1 
Hence 


K K 
dA*—> Aiidn; => pvidn? (48) 
i=1 i=1 
which by Eq. (26) leads directly to (40). To 
obtain (41) the process is made finite instead of 
infinitesimal and the operator d in the formulas 
(42)-(48) replaced by A. 

Another important relationdirectly derivable 
from the engine is Eq. (37). For this purpose 
cause B* to generate the finite surface area Q, 
simultaneously moving B*, B®, Be, bf (and 
b in the case of insoluble surface films, in 
such a way as to keep all intensive variables, 
the volumes V2, V8, the level of the Gibbs 
surface and the level of the surface layer con- 
stant. This process is best carried out in four 
Steps: (i) generation of 2 by B*, whereby V*, 
Ve, P and the concentrations c,*, cf are changed; 
(ii) restoration of V« and V® to their original 
values by B« and Bé respectively ; (iii) restoration 
of c*, c (and the surface densities of any com- 


ponents present only in the surface layer) and 
therefore P to their original values by means of 
(and (iv) restoration of the surface 
layer to its original level by means of );*, 5, as 
in the foregoing derivation of Eq. (40). From 
Eq. (42) it follows that for any finite change 


AG = V2) +A(g8 V8) +AG*. (49) 


Since in the change in question V*, V® and all 
intensive variables remain constant and the 
surface 2 is generated, 


AG = AG* =G* (50) 

so that in particular 
AA=A* (51) 
An;=n (52) 


where An; is the total amount of the ith com- 
ponent added to or removed from the system 
aBw in steps (iii) and (iv) above. The gain in 
the Hfe. of the entire engine is therefore 


K 
Aine (53) 


i=1 


and the work done on the entire engine is 


K 
oQ+ (54) 


i=1 


Hence 


K 
i=1 
which by Eq. (26) leads directly to (37). 

It may furthermore be noted that just as 
Eq. (35) for a surface layer was integrated to 
give Eq. (37), so may Eq. (34.1) for a single 
phase a@ be integrated keeping all intensive vari- 
ables constant and allowing only V* to vary; 
the result is 


K 
At=—P Ve+>. (56) 


i=1 
This equation, like (37), may also be derived 
directly from the engine by making process 2 
above finite and keeping all intensive vari- 
ables constant. Finally, Eq. (29) may be inte- 
grated in a manner similar to Eqs. (35) and 
(34.1), (or (34.2)); the result is the complete 
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expression for the Hfe. of the system afw: G. COMPARISON OF THE EQUATIONS DERIVED 
FROM THE ENGINE WITH THE GENERAL 
= a B B 
A GIBBSIAN EQUATIONS 


K K K 
+¥ wn’ (57) The differential equations for the Hfe. derived 
= i ates in this paper apply only to isothermal changes, 
but are otherwise quite general for the cases 


K 
= —PV+o2+) mini (58) defined by the above assumptions 1 and 2. 
It therefore seems worth while comparing these 

and these equations too are immediately obtain- equations with the completely general equations cz 

able from the engine by allowing it to generate which take account of changes of temperature. wv 


the whole system a$w keeping the intensive The latter equations, obtainable by the method 
variables constant. Eqs. (57), (58) are of course of Gibbs, differ from the isothermal equations 


also obtainable by combining (37) with (56) only by an additive term of the form SdT. Thus " 
and with the integrated form of (34.2). the general equation corresponding to (29) is , 
st 
K K K al tk 
dA*+dA8+dA* = —S¢dT — — — Pd V*— Pd V8+od24+ (59) | 
i=1 i=1 i=1 
a 
or, more compactly, Finally, just as the combination of (37) with . “ 
K the isothermal equation (35) leads to the _ ‘ 
dA = —SdT—Pd V+odQ+> widn;. (60) isothermal Eq. (2), so does the combination | 
7 of (37) with the general Eq. (61) lead to the | 
The general equation corresponding to (35) is: general Eq. (1). . 

dA* = — S*:dT+od2+ (61) H. SUMMARY 

i=1 

and to (34.1): An exact elementary deduction of Gibbs’ mi 
» adsorption theorem, be 
dA*= —S*dT—PdV*+ > pidn;*. (62) x an 
.: hg It is worth noting that the integration of the es a 
es general Eqs. (59), (60), (61), (62) keeping all and some related equations is given by means of ” 
an ¢ intensive variables constant leads to Eqs. (57), a suitable thermodynamic engine. The assump- ” 
ane (58), (37), (56) already obtained in this paper; tions underlying the Gibbsian treatment of the | ” 
ian the integrated equations derived from the engine thermodynamics of surfaces are discussed briefly. | oi 
ee: therefore do not suffer from the fact that the The relation of the equations derived to the cor- a 
corresponding differential equations apply only responding general Gibbsian equations is ex | ~ 


to isothermal changes. plained. 
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A Calculation of the Equilibria in Keto-Enol Tautomerism 


G. W. WHELAND, Gates Chemical Laboratory, California Institute of Technology* 
(Received July 20, 1933) 


A calculation of the equilibria obtaining in some simple 
cases of keto-enol tautomerism has been carried through 
with the aid of bond and resonance energies derived from 


empirical data. The results obtained are in qualitative and 
semi-quantitative agreement with experiment. 


HE energy of a molecule which can be 

represented by a single Lewis electronic 
structure is given with considerable accuracy by 
the sum of the energies of the various bends 
which it contains. On the other hand, the energy 
of a molecule which can be represented ade- 
quately only in terms of two or more structures 
is, as a result of resonance, greater than would be 
predicted for even the stablest of them, on the 
basis of bond energies alone. The amount of this 
extra energy has been calculated from empirical 
data by Pauling and Sherman! for a number of 
cases, and has been shown to be fairly constant 
for a given type of resonance. We can, conse- 
quently, calculate the energy of an arbitrary 
molecule if we know the energy of each of its 
bonds and in addition the contribution made by 
any resonating systems that may be present.? 
Unfortunately, however, we do not have suffi- 
cient accurate data of this kind to proceed very 
far except in a few simple cases, and so we must 
usually resort to rough estimates and content 
ourselves with general results of a more or less 
qualitative nature. The purpose of the present 
paper is to apply this type of treatment to the 
problem of keto-enol tautomerism. The results 
obtained are quite satisfactory qualitatively and 
semi-quantitatively, and are of such nature as to 
lead one to expect that a more rigorous treatment 


*Communication No. 373. 

"L. Pauling and J. Sherman, J. Chem. Phys. 1, 606 
(1933), All bond and resonance energies given without 
reference will be taken from this paper. 

* This is rather an over-simplification since it leaves out 
of account the energy resulting from interaction between 
the dipoles in the molecule. This effect is usually small, 
and so we shall neglect it in the following discussion, except 
in a few cases where it becomes of special importance. 


would give complete agreement with experiment. 
Since, however, this will have to wait until more 
powerful methods than are now available have 
been developed, it seems of interest in the mean- 
time to carry the rougher and more qualitative 
methods as far as possible. 

The calculations are simplified to a certain 
extent by the fact that the entropy change in the 
reaction keto->enol is apparently very small.’ 
We can consequently neglect it entirely and put 
AH=AF. This is convenient since it makes 
possible a direct comparison between the values 
of AH given by the calculation and the values of 
AF derived from experimental data. In only a 
few cases are experimental A// values known. 


SIMPLE KETO-ENOL SYSTEM 
H-—C—C=0O C=C-—-O-H 
I 


In the keto grouping I there is no possibility 
of resonance, and so the energy is simply the sum 
of the bond energies. This amounts to 357.2 Cal. 
if the C=O group is aldehydic, and to 360.6 Cal. 
if it is ketonic. In the enol grouping II, resonance 
between the two structures 

70:8 
is possible, but presumably the second structure 
does not make a very large contribution and the 
resonance energy is small. That some resonance 


does occur, however, is made probable by the 


3J. B. Conant and A. F. Thompson, Jr., J. Am. Chem. 
Soc. 54, 4039 (1932). Most of the data regarding the equi- 
libria with 8-diketones and f-ketonic esters are taken from 


this paper. 
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fact that the enols are considerably stronger acids 
than the alcohols. The bond energy for the enol 
grouping is 339.9 Cal. if the energy of the C—O 
bond is the same as for primary alcohols. 
Actually it is probably somewhat higher since an 
enol bears a rather closer resemblance to a 
secondary or tertiary than to a primary alcohol. 
The total energy of the enol grouping is then 
slightly greater than 339.9 Cal—probably about 
346.8 Cal. since in phenol similar effects raise the 
energy by 6.9 Cal. Thus ketones should be 
approximately 14 Cal. and aldehydes 10 Cal. 
more stable than the corresponding enols. This 
is in qualitative agreement with experiment; 
since the enols of this type are too unstable ever 
to have been isolated, a quantitative comparison 
cannot be made. 


PHENOL 
One enolic and two ketonic forms are possible. 


H H H 

if i if 

11 IV V 


(In the structure III the benzene ring is written, 
for the sake of convenience, as if it possessed the 
simple Kekulé form. Actually, of course, it 
resonates among several structures, of which this 
is only one.*) As far as the keto and enol group- 
ings are concerned, the two keto forms IV and V 
should be more stable than the enol III by 14 
Cal. The resonance energy remains to be con- 
sidered, however. For the enol form this amounts 
to 39 Cal. (resonance within the benzene ring), 
and for the two keto forms to considerably less, 
probably between 13 and 19 Cal. The reason for 
this is as follows. The resonance energy of the 
system C=C—C=C can be calculated by the 
method described by Pauling and Wheland, and 
is found to be approximately 8 Cal. The reso- 
nance energy of the system C=C—C=O is 
probably a little less than 10 Cal.5—that is to 
say, practically the same as for C=C—C=C. 


*L. Pauling and G. W. Wheland, J. Chem. Phys. 1, 362 
(1933). 

5 This value is derived from data quoted by W. Hiickel, 
Theoretische Grundlagen der Organischen Chemie, Vol. I, 


Consequently it seems reasonable to assume that 
the resonance energies of the systems C=C~C 


C=Q 
=C—C=O and poe lie close to the 


values 17 Cal. calculated for C=C—C=C-—-C 


=C and 15 Cal. calculated for Pew, re- 
C= 

spectively. It follows then that the enol form of 

phenol is more stable than either of the keto 

forms by 6-12 Cal. Again the agreement with 

experiment is qualitatively satisfactory, but a 

quantitative comparison is impossible on account 


of lack of experimental data. 


RESORCINOL 


There is one dienolic and one diketonic form. 


O-H 
H 
H He 
VI VII 


(The intermediate keto-enolic forms turn out to 
be less stable than either of the above.) If the 
resonance energy is neglected the diketonic form 
VII would be expected to be more stable than 
the dienolic VI by about 28 Cal. When the 39 
Cal. resonance energy for the benzene ring and 
the 10 Cal. for the system C=C—C=O are 
considered, however, the dienol becomes more 
stable by 1 Cal. This is in qualitative agreement 
with the fact that resorcinol shows more of 4 
ketonic behavior than phenol, but that its 
phenolic character is still predominant. 


PHLOROGLUCINOL 
There is one trienolic and one triketonic form. 
O-H O 
| 
H- ) -H He SH: 
O= =0 
H He 
VIII IX 


p. 304, Leipzig, 1931, regarding the heats of combusti | 
of isomeric a8 and By unsaturated ketones. 
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EQUILIBRIA IN KETO-ENOL TAUTOMERISM 


A calculation similar to the above gives the 
result that the triketonic form IX should be 3 
Cal. more stable than the trienolic form VIII. 
This is again in good qualitative agreement with 
the fact that phloroglucinol possesses both 
phenolic and ketonic character, but the quanti- 
tative aspect is less satisfactory, as the phenolic 
character is apparently still predominant. The 
agreement would have been somewhat better if 
the dipole interactions had been considered. No 
great weight should be attached to the actual 
numerical values derived in this and the pre- 
ceding calculations as they could hardly be 
expected to be accurate within a few kilocalories. 
The relative values, however, are probably 
significant. 


ACETYLACETONE 


This is chosen as a typical example of a £- 
diketone. The keto and the enol forms (X and 
XI, respectively) are both known, and the free 
energy change accompanying the conversion of 


CH;—C—CH=C—CHs 
Xx XI 


the keto into the enol form has been found to be 
—1 Cal. With compounds of this type an 
additional complication enters in the fact that 
the enolic modification does not have the simple 
structure given above, but contains a chelate 
ring with a hydrogen bond. The resulting com- 
pound is then free to resonate between the two 
equivalent structures: 

O— H O-— H 


| | | 
CH;C-CH=C-—CH; CH;C=CH—C-—CH; 


The amount by which this bond formation and 
resonance stabilizes the enol form can be calcu- 
lated from the fact mentioned above that AF for 
the reaction keto>enol is — 1 Cal. The keto form 
is readily seen to be 4 Cal. more stable than the 
open-chain enol form, and hence the formation 
of the chelate ring must stabilize the latter by an 
additional 5 Cal. This value is reasonable enough, 
but unfortunately cannot be checked by an 
independent calculation. 
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There is no provision in this treatment for the 
effect of substituent groups in the molecule 
(provided that at least one enolizable hydrogen 
is left), and in fact the calculation leads to the 
conclusion that the enolization should be inde- 
pendent of the presence of substituents. Actually 
this is not the case, but the effect of substitution 
is very small—less than 2 Cal. in all cases—and 
lies well within the accuracy of the method. 


OXYMETHYLENE KETONES OR 
B-KETONIC ALDEHYDES 


In compounds of this type enolization can 
take place in either of two ways and consequently 
there are (before formation of the chelate ring) 
two enol forms (XIII and XIV) in addition to 
the one keto form XII. 


R-C-—C-C-H R-C-C=C-H 
|/ 

R R 
XII XIII 
R-C=C-C-H 
R 
XIV 


Since the energy of the C=O bond in aldehydes 
is 3.4 Cal. less than in ketones, it follows that the 
stability of the enol XIII is greater than that of 
XIV and nearly as great as that of the keto form 
XII. When the hydrogen bond is formed the 
resonance in the chelate ring will be less than in 
the case of the diketones as a result of the non- 
equivalence of the two structures 


R-C=C-C-H 
| 7 
R R 
XV XVI 


The effect of this will be twofold. In the first 
place, the energy of the ring will be less than 5 
Cal., but probably still enough to make the enol 
form more stable than the keto; and in the 
second place, the structure XV will predominate 
over XVI since it represents the more stable 
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configuration. Experimentally, compounds of 
this type (for example, oxymethylene camphor) 
are found to possess the oxymethylene ketone 
structure, corresponding to the enol XIII or XV, 
in qualitative confirmation of the calculation. 
There are no accurate experimental data of 
quantitative nature for purpose of comparison. 


ACETOACETIC ESTER 
A £-ketonic ester can give rise to two (open- 
chain) enolic modifications (XVIII and XIX) in 
addition to the ketonic (XVII). 


O O-H O 


|| | 
CH;—C—CH,—C-—OEt CH;—C=CH—C-—OEt 
XVII 
O H-O 


| | 
CH;—C—CH=C-OEt. 


XIX 


This problem is complicated by the fact that the ° 


carbethoxy group has a resonance of its own 
which makes the relatively large contribution of 
28 Cal. to the energy of the molecule, and that 
this resonance is interfered with to an uncertain 
extent by any other types of resonance that may 
be present. In the keto form this effect does not 
enter since there is no other resonance than that 
of the carbethoxy group. In the enol form XVIII 
it is apparently of very minor importance since 
in a few cases for which data are available® the 
resonance energy of af unsaturated acids is 
essentially the same as if the C=C—C=O and 
the carboxyl groups were entirely independent of 
each other. In the enol form XIX, however, the 
carbethoxy group and its resonance energy have 
been destroyed, and hence this structure repre- 
sents a relatively unstable molecule. The con- 
clusions to be drawn from these considerations 
are: first, that the 8-ketonic esters should be 
somewhat less enolic than the 8-diketones, as a 
result mainly of damping of resonance in the 
chelate ring; and second, that the first of the 
resonating structures of the chelate compound 


CH;C =CH—C—OEt CH;—C—CH=C-—OEt 


* For example, the heats of combustion of af and By- 
pentenoic acids differ by 8.5 Cal. See M. S. Kharasch, 
Bur. Standards J. Research 2, 359 (1929). 


G. W. WHELAND 


should predominate over the second. Both con- 
clusions are in agreement with experiment. 
Acetoacetic ester is less enolic than acetyl 
acetone by about 1 Cal. (AF=O Cal.), and the 
position of the C=C double bond has been 
determined by ozonization.’? The 8-ketonic esters 
are perhaps somewhat more enolic than might be 
expected from the discussion, but the agreement 
with the calculation is really quite satisfactory. 


TRIACETYLMETHANE 
In this case there is only one keto and one enol 
form (XX and XXII, respectively). 
O O O H-O 
| | | | 


CH;C —CH —C—CHs CH;C—C=C—CH; 


| 
C=O C=0 


| | 
CH; CH; 
XX XXI 


We see that this compound will be more enolic 
than acetylacetone by the amount by which the 
C=O 

exceeds 
C=O 
that of the system C=C—C=O. An estimate of 
about 4-8 Cal. seems to be reasonable for this 
figure and to agree satisfactorily with what 
little experimental data there are. The compound 
is usually considered to be completely enolic, but 
the equilibrium has never been measured accu- 
rately. 


resonance of the system C =c¢ 


ACETYLMALONIC ESTER 


C=O 
| 
O O 
Et Et 
O 
Et 


7 J. Scheiber and P. Herold, Ann. d. Chemie 405, 295 
(1914). 
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EQUILIBRIA IN KETO-ENOL TAUTOMERISM 


This compound would be expected to be less 
enolic than triacetylmethane (for the same 
reason that acetoacetic ester is less enolic than 
acetylacetone) and to be more enolic than 
acetoacetic ester. Experimentally, AF of enoliza- 
tion is found to be about —0.5 Cal. A value 
between —4 and —8 Cal. would seem more 
reasonable, but the qualitative agreement is 


good. 


TRICARBOMETHOXY METHANE 


O O 
Me Me 
XXII XXIII 


The stability of the enol form XXIII is 
decreased very materially by the fact that in it 
one of the carbomethoxy groups is no longer 
intact. Just how great this effect will be, however, 
cannot be predicted with any assurance since, as 
a result of the resonance within the chelate ring, 
the relationships involved are very complicated. 
All that can be said definitely is that this 
compound should be less enolic than acetyl- 
malonic ester. This is indeed the case, but the 
difference is very much less than might be 
expected. The enolic form is less stable than the 
ketonic by only about 1 Cal., while a value of 
10-20 Cal. would seem more reasonable. The 
quantitative aspects of the present treatment are 
consequently very unsatisfactory at this point, 
although the qualitative agreement is still good. 
Two possible explanations for this discrepancy 


‘Suggest themselves at once. In the first place, the 


value of 28 Cal. for the resonance energy of the 
carbomethoxy group may be considerably too 
high. This figure was arrived at! on the doubtful 
assumption that the energies of the C=O and 
the C—O bonds were the same as in ketones and 
primary alcohols, respectively, and consequently 
may very well be in error by several kilocalories. 
In the second place, the neglect of the large 
interaction between the dipoles in the carbo- 
methoxy group is probably not legitimate. 
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MALOonIc ESTER 


Et—O—C—CH;—C—OEt 
O H-O 


| | 
Et—O—C—CH =C-—OEt 


This case differs from the preceding only in 
the substitution of the resonating system C=C 


C= 
—C=O for c=c¢ in the enol form. Con- 
C=O 


sequently the stability of the enol relative to that 
of the keto should be decreased by the difference 
between the resonance energies of the two 
systems, i.e., by 4-8 Cal. This is in good agree- 
ment with the fact that AF for the reaction 
keto—enol is known certainly to be greater than 
6 Cal. and probably to be less than 14 Cal.® 


CYCLOHEXADIONE 1,2 


Three structures are possible. 


Hi H, Hx SH SH 
H: He 2 Hz He 
XXIV XXV XXVI 


The calculation is simpler than in the case of 
B-diketones since the enol forms of a-diketones do 
not form chelate rings, partly for steric reasons 
and partly on account of the impossibility of 
resonance within the ring. The resonance energy 
of the system O=C—C=O is not known, but 
presumably is not very different from that of 
C=C—C=O or of C=C—C=C. Consequently 
it would appear that the diketo form XXIV 
should be considerably more stable than the 
keto-enol XXV, which in turn should be more 
stable than the dienol X XVI. Actually, however, 
the compound appears to exist only in the keto- 
enol form. The reason for this discrepancy is 
probably to be found in the potential energy 
which results from the dipole repulsions of the 
carbonyl and the hydroxyl groups present, and 
which cannot be neglected in this case. The two 


8 J. B. Conant and G. H. Carlson, J. Am. Chem. Soc. 
54, 4054 (1932). 
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dipoles are close together and are held in the cis 
position by the ring formation. The magnitude of 
this effect cannot be determined with any 
accuracy on account of the uncertainty in regard 
to the position and orientation of the dipoles and 
the dielectric constant of the medium. A rough 
calculation, however, shows that the effect is in 
the right direction and is of the correct order of 
magnitude.® It is probable that similar effects 
are responsible for other discrepancies in quanti- 
tative comparisons noted above. 


DIACETYL 


This case differs from the preceding in that the 
molecule is not kept rigid by the presence of a 
ring. Hence the dipoles are free to orient them- 
selves in the more stable trans position where 
their effect upon the stability of the molecule is 
relatively small. This is in accord with the fact 
that the molecule exists only in the diketo form. 


CONCLUSION 


The results which have been obtained for the 
simple tautomeric systems treated above are in 
complete qualitative agreement with experiment 


® For a somewhat similar calculation see H. M. Small- 
wood, J. Am. Chem. Soc. 54, 3048 (1932). 


G. W. WHELAND 


and furnish considerable justification for the 
semi-empirical treatment used. In only two cases 
(tricarbomethoxy methane and malonic ester) 
are the quantitative aspects of the calculation 
unsatisfactory. These two compounds differ from 
the rest in that, with them, enolization necessarily 
destroys the characteristic resonance of a car- 
boxyl group. As has been pointed out, this fact 
makes the problem very complicated, and so it is 
not surprising that quantitative agreement is 
lacking at this point. It is doubtlessly significant, 
however, that the calculations for the two com- 
pounds are in good quantitative agreement with 
each other. 

This same treatment could be extended to a 
number of other cases, but since its success 
depends upon the possibility of estimating the 
resonance energies of various groupings in the 
molecules concerned, it rapidly becomes of less 
value as the complexity of the resonating systems 
increases. For this reason it does not seem 
profitable to attempt the treatment of such cases 
as dibenzoylmethane, etc., where the presence of 
the phenyl groups complicates the problem 
tremendously. The success obtained with the 
simpler systems, however, makes it seem quite 
probable that these could also be treated satis- 
factorily by the same method if sufficient 
empirical data were at hand. 
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Further Observations Concerning the Propagation of Sound in Nitrogen Tetroxide 


T. RicHarDs* AND JAMEs A. Rein, Princeton University 
(Received June 12, 1933) 


The absorption of sound in nitrogen tetroxide is dis- 
cussed, and it is pointed out that absorption measurements 
should be capable of distinguishing whether a given 
absorptive region is due to failure of the dissociation reac- 
tion or of the heat capacity to follow the adiabatic cycle 
of the sound wave. Appropriate expressions for the absorp- 
tion are derived from Einstein’s theory, and are briefly 
discussed. An experimental method of determining the 
absorption maximum is described in which the tempera- 
ture is varied at constant pressure and frequency. Under 
these conditions the acoustical variables are maintained 


very nearly constant. An attempt to detect the absorption 
maximum at 450.1 kc and 260 mm in nitrogen tetroxide 
was unsuccessful because of the inaccuracy of the absorp- 
tion measurements. It is concluded that dispersion meas- 
urements provide at present an easier and more accurate 
method of investigating the acoustical properties of rapidly 
dissociating gases than measurements of the absorption 
coefficient. Occasion is taken to discuss the tube corrections 
applied to the previous dispersion measurements at greater 
length than was before considered necessary. Several 
corrections to the previous work are also appended. 


INTRODUCTION 


RECENT communication from this labo- 

ratory' described measurements on the 
dispersion of sound in nitrogen tetroxide, and 
discussed in considerable detail an argument by 
which it was believed legitimaté to interpret 
them as due to failure of the dissociation reaction 
to attain equilibrium in the sound wave. The 
chief purpose of this note is to extend the 
argument to acoustical absorption. 

Our report has been followed by an article by 
Teeter? in which he devotes much space to 
unfavorable criticism of our work. We shall take 
this opportunity to answer such parts of Teeter’s 
attack as appear to us to require discussion. The 
bulk of his criticism is, however, based on 
arguments which we consider that we have 
already adequately refuted, and it would serve no 
useful purpose to labor them further. 


Part I. A DiscussIon OF THE ABSORPTION OF 
SouND IN DtssociATING GASES 


The interpretation of absorption measurements 
in terms of dissociation rates 


It has already been demonstrated! that suita- 
ble measurements of the dispersion of sound in a 

*At present Fellow of the John Simon Guggenheim 
Memorial Foundation. 


' Richards and Reid, J. Chem. Phys. 1, 114 (1933). 
Teeter, J. Chem. Phys. 1, 251 (1933). 
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dissociating gas are capable of showing whether 
(a) the dissociation reaction, (b) the heat 
capacity, or (c) both the dissociation reaction and 
the heat capacity concomitantly are failing to 
attain equilibrium in the adiabatic cycle. This 
demonstration is not wholly free from objection, 
since it is based on a somewhat oversimplified 
picture of the conditions in the path of the 
sound-wave.* However, since it appears to repre- 
sent the most satisfactory description of the 
phenomenon which is at present available, it is 
perhaps worth while to extend the treatment to 
absorption. In order to do this with as little 
complication as possible it is proposed to treat 
quantitatively only cases (a) and (b), leaving the 
somewhat more doubtful case (c) for later 
qualitative discussion. In developing the argu- 
ment it will be convenient to review that 
previously given for cases (a) and (b) in con- 
nection with the dispersion. 

Calling « and A the complex pressure and 


3 The assignment of a single characteristic time to the 
rate of adjustment of the heat capacity of a gas in which 
several quantum states are excited is at present without 
theoretical justification. See Bourgin, Phys. Rev. 42, 721 
(1933). It seems, nevertheless, to be a reasonable experi- 
mental expedient. See Richards and Reid, Nature 130, 739 
(1932). A more completely deduced description of the 
concomittant failure of dissociation reaction and heat 
capacity may, for the present, be expected to yield results 
of substantially the same form as that given here and in 
the previous communication. 
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density variations of the sound wave, the theory 
of Einstein‘ gives for a dispersive gas 


« Pr AotifA, 
A p Bot+iB,, 


(1) 


where P and p are the equilibrium pressure and 
density, respectively, and 


f=o/¢ (2) 


if w is the frequency of the sound wave in radians, 
and ¢ is the characteristic frequency of the 
process which produces the dispersion. In the 
case of a unimolecular reaction ¢=k,, the rate 
constant. For a lag in heat capacity ¢ is pro- 
portional to the number and effectiveness of 
the collisions which excite an active state of 
the molecule. The quantities 1+A o/Bo and 1 
+A,,/B., represent the ratio of the isothermal to 
the adiabatic compressibilities when f sensibly 
vanishes and when f becomes sensibly infinite, 
respectively. The subscripts 0 and ~ will be used 
in this sense throughout. 

Since, as will later become apparent, the 
absorption due to the process characterized by ¢ 
becomes negligible at very low and at very high 
frequencies, the real part (r.p.) of +/A becomes 
equal to the square of the phase velocity V. Thus 
from (1) the result is obtained that 


Vo? = (3) 


and the useful abbreviation Q, which is inti- 
mately connected both with dispersion and ab- 
sorption, is defined. The relationship (3) is alone 
sufficient to determine whether (a) or (b) has 
produced a given dispersive region, since the 
maximum change in velocity will be quite 
different for the two unless the dissociation heat 
is small. In the case of nitrogen tetroxide at 
25°C and 260 mm, for example, the degree of 
dissociation a is about 0.3, and 


4 Einstein, Sitz. Ber. Akad., 380 (1920). 

5 Kneser, Ann. d. Physik 11, 761 (1931), introduced the 
rate constant which we employed in the previous work, 
and shall continue to use here. Because the notation of 
Einstein has been retained wherever possible Kneser’s 8 
has been replaced by 6, both meaning identically the 
average lifetime of the excited molecule. Thus ¢@=1 in 
our notation. In Kneser’s, =1. 
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Case (a) 

2D _\1-a 

Ay= ) 
1+a 

A,.=Ap=R =2 cal. 

4 

D? 1-a 5(— 


2-—a 
+R(—*) =56.4 cal. 
a 


By=—— 


RT* 1+a 
B..=Bp= =10 cal. 


) =651 cal. 


a 


according to the result of Einstein. This is on the 
assumption that the heat capacity continues to 
follow the wave without lag, as in case (a). For 
case (a) then, Q?=Qp*= 1.105. In case (b) A» and 
By remain unchanged, but the heat capacity 
alone drops out at high frequencies. To describe 
this it is sufficiently definite to write Co=C+C, 
where Cy represents the average heat capacity of 
the N2O,—NO, mixture at equilibrium and at 
constant volume, and C and C, are its frequency- 
independent and -dependent parts, respectively. 
Cy is probably about 10 calories: mole, and the 
utmost change of the heat capacity in regions of 
frequency which are at present accessible may be 
expected to be a drop to half this value. Thus C 
will be taken as 5 calories: mole, and A,, and B, 
become 


Case (b) 
[Ay 256.4 cal. 
2D _\1-a 2-—a 
A,=A-.= 259.1 cal. 
a 
By 2651 cal. 
D? 1-a _/2-a 
B,.=B,.=— ) =623 cal. 
RT? 1+ea a 


and for case (b) @?=Q.? =1.008. 

Thus it is seen that approximately ten times as 
much change in velocity results from case (a) 
as from case (b). A change of velocity greater 
than 0.4 percent, ((Q.—1) 100), must, there- 
fore, involve a partial failure of the dissociation 
reaction. It was an extension of this argument to 
case (c) by means of a simplified model for the 
heat capacity which made possible the calcu- 
lation of rate constants from the measurements 
reported previously.! 
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PROPAGATION OF SOUND IN N:2O, 


In order to evaluate the corresponding maxima 
of the absorption coefficient it will be necessary 
to examine the structure of the complex pressure 
and condensation variations, 7 and A. For these 
Einstein‘ has written 


T= (4) 
A = (5) 


where the amplitudes are moe’* (complex) and 
Ay (real), ¢ is the time, x the distance in the 
direction of propagation, and a@ a complex 
quantity involving both velocity and absorption. 
It is apparent that the imaginary part of a 
becomes the only real part of the two exponents, 
and hence is a damping coefficient. If the symbol 
6B denotes the absorption coefficient for the 
amplitude expressed as a reciprocal length, then a 
becomes 


a=(1/V—1B/w) (6) 


where 2xV=wd and i=(—1)}. The familiar 
relationships, in which u denotes displacement in 
the direction of the wave normal, 


dn/dx= —pd*u/dt? and A=.—pdu/dx (7) 


yield a wave equation of the usual form in which 
7/A is the square of a complex plane wave 
velocity, and has the value 


(8) 


Eq. (1) is that from which the velocity and 
absorption of sound are to be determined in 
terms of the physical properties of the gas. It 
may readily be seen that the real part of 7/A is 
more easy to obtain, from (1), than the real part of 
a (the real part of (A/x)!). Thus the convenient 
solution for the velocity is 


1—6°V?/w? 
| 
A (1+6°V?/w*)? 


(9) 


Similarly the imaginary part (i.p.) of A/z gives 
conveniently the absorption coefficient, as is 
apparent from (1), (6) and (8). Thus 


ip. A/r= —i28/wV. (10) 


Up to this point nothing has been added to the 
argument of Einstein. 


If, now, the absorption coefficient for the 
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intensity expressed as a reciprocal number of 
wave-lengths is introduced according to the usual 


definition 
1, = = 


(10) becomes 
h=—(22rV?/i)(i.p. A/7). (11) 


The imaginary part of the reciprocal of (1) gives 


Bi(A,+B..)f —B..(AotBo 
p (12) 


h=2rV? — 
P\ 


For the sake of convenience the réle of absorption 
in the real part of A/z will be neglected, and it 
will be taken that 


rp. A/r21/V? 


and, consequently, from the reciprocal of (1), 
that 


p 


(13) 
By(Aot+ Bo) +B,,(A.+B.)f2) 


The approximation so introduced amounts to 
neglecting terms in h’, or, practically, to setting 


(r.p. r/A)(r.p. A/r) 21, 20? =0*+1 


which, in the present condition of interferometric 
absorption measurements, will be without influ- 
ence.® Substituting (13) into (12) gives finally 


(14) 


where y, which has been suggestively called the 
“relative frequency” by Luck,® is defined as 


y (14a) 
Herzfeld and Rice,’ Bourgin,®’ Henry® and 


IIS 


® Luck, Phys. Rev. 40, 440 (1932) has carried through 
the solution for a dissociating gas without a heat capacity 
lag (case (a)) without making this approximation. Our 
(14) and (20) for this case are a slightly less accurate 
version of his 


The approximation does not enter (15), and (16) has been 
given in both the approximate and the exact forms. 
7 Herzfeld and Rice, Phys. Rev. 31, 691 (1928). 

8 Bourgin, Phil. Mag. (vii) 7, 821 (1929). 
® Henry, Proc. Camb. Phil. Soc. 28, 249 (1932). 
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Kneser!® have given similar equations for non- 
dissociating gases with variable heat capacity. 
The exposition used above was clearly implied by 
Einstein, and has been here repeated only in 
order to give it somewhat more general signifi- 
cance. 

The maximum of the absorption coefficient h 
may of course be found by setting dh/dw=0. It 
occurs at a frequency 


Wmax = £Bo/QB,, 
and has the value 
Iinax (16) 


as Luck,® Henry® and Kneser'® have pointed out 
in appropriate cases. 
If (1) is written without regard for terms in h? 


w 2(6Bo/B..) {(Va?— (17) 


and @max from (15) is substituted therein, there is 
obtained for the velocity at the point of maxi- 
mum absorption 


V + V0?Q?) /(Q? + 1) 


which Luck,’ Henry? and Kneser® have also 
reported. 

Remembering that at 260 mm and 25°C in 
nitrogen tetroxide for case (a) Q=Qp = 1.051 and 
for case (b) Q=Q.=1.004, the maximum ab- 


(15) 


By—w0C+i[ 
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sorption coefficient for the dissociation reaction 
hpanax) is seen to be about 0.312 when the heat 
capacity remains fully active! and the maximum 
absorption coefficient for the heat capacity 
he«max) is approximately 0.030 if the dissociation 
reaction continues to follow the wave without 
lag. Here again case (a) gives about ten times the 
effect of case (b), and an experimental distinction 
between the two should easily be possible. 

Turning now to case (c) where both heat 
capacity and dissociation contribute simultane- 
ously to the absorption we shall make use of the 
same approximate model as was previously em- 
ployed for the dispersion. It will be supposed 
that a single time constant @ is sufficient to 
describe the dispersive behavior of the heat 
capacity, and that the frequency dependence of 
the heat capacity conforms to the expression 
first suggested by Kneser® where C, 


C+C.(1/(1 +iw)) 


is the contribution of the heat capacity to the 
complex velocity of sound at a frequency »w. 
Substituting this in the equation of Einstein in 
the usual way gives the complex velocity of 
sound as in Eq. (6) of the previous communi- 
cation. For the present purpose the reciprocal of 
this may be written 


+o(R+Co)] 


| (18) 


where A», Bo, A-, and B, have the values given above. The absorption coefficient in case (c) is 


then given by the expression 


hoc 


This is not only unsightly, but has the disad- 
vantage that it does not readily yield expressions 


10 Kneser, Ann. d. Physik 16, 337 (1933). Kneser has 
given (14) for the case Byp=C;+C, and B,=C, in his 
notation. He has obtained (14) from the tangent of the 
phase angle ¢ by the approximation 2” tan g=h. This 
may be seen to be identical with the approximation here 


introduced if tan g is obtained from Eqs. (1) and (2) of 


our previous communication by neglecting terms in £*. 


+B.) +o(R+Co)]— Co) [hi (A 0+Bo) —w0(R+ (19) 
0+ Bo) —w0(R+C)]+ +B.) +o(R+Co) 


of the type of (15) and (16). Nevertheless it must 
be briefly discussed in view of the experimental 
work which follows. 

In case (a), when the heat capacity is without 
lag, 


[6=0; yo =(w/ki)-(Co/Bo)], (19) reduces to 
hp 22ryp{ } (20) 


1 Luck (reference 6) was the first to point this out. 
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(18) 


is 


and in case (b) when the dissociation reaction is 
without lag, 


(ki>>w, 0-1; ye-=w0B./Bo), (19) reduces to 
he (21) 


If the dissociation is so slow that it plays no part 
in the compressibility ratio (k;=0, f=w6), (19) 
gives Kneser’s expression'® applied to a gaseous 
mixture of very particular properties: 


C Q?-1 
Col1+Q*%(Cf/Co)? 
This case is without interest at present, since 
frequencies so high that the dissociation reaction 
has been wholly eliminated have not yet been 
employed to study nitrogen tetroxide. The 
magnitudes of the absorption maxima corre- 
sponding to (20) and (21) have already been 
discussed. 

It appears best at present to limit the treat- 
ment of intermediate cases to a qualitative 
paragraph. In general, as might be expected, two 
absorption maxima, corresponding to two sepa- 
rate dispersive regions, appear when k,; and @-! 
are different by more than two or three orders of 
magnitude. The two maxima are then roughly 
independent of each other. If 2; is small compared 
to 6! the two separate peaks are of about the 
same height, corresponding to the two appro- 
priate forms of (16). If, on the other hand, the 
heat capacity lag occurs at a lower frequency 
than that of the dissociation reaction, i.e., 6~' is 
small compared to k;, the absorption from the 
former is small, but that due to the latter 
considerably greater than (20) indicates, since 
the situation will then approximately be 

wo C 


= 
ky B. 


where Q has a value greater than Qp by an 
amount nearly proportional to the ratio (Co/C)!. 
As ki and 6 approach each other the two 
absorption peaks merge into one which is 
conspicuously higher than either alone. It ap- 
pears reasonable to suppose that in the limiting 
case where both dissociation and heat capacity 
fail at an identical frequency of sound the 
Maximum of the absorption coefficient will be 


wa 
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given by 


B(R+C) 
Apc max) 1/Q), where 


If this is true, Apcanax) Will rise to 0.8 in nitrogen 
tetroxide at 25°C and 260 mm provided that C is 
reduced to 5 calories. In case the heat capacity 
fails at lower frequencies than the dissociation 
reaction or at the same frequency, wpqnax) will 
shift to higher frequencies in rough conformity to 
(15). An indirect consequence of this has already 
been illustrated by the dip in the curves of Fig. 1 
of the previous communication. 


An attempt to determine the absorption maxi- 
mum by experiment 


Since it appears from the previous discussion 
that the magnitude and position of the ab- 
sorption maximum is influenced by heat capacity 
effects, it becomes reasonable to inquire whether 
absorption measurements may not conveniently 
supplement the information already gained from 
dispersion. Two publications have already re- 
ported observations concerning the absorption of 
sound in nitrogen tetroxide, but they do not, 
unfortunately, permit even qualitative con- 
clusions. 

Kistiakowsky and Richards’ quote measure- 
ments of the logarithmic decrement coefficient of 
the d.c. plate current of the amplifier coupled to 
the magnetostriction oscillator as a function of 
pressure and frequency. This has been shown to 
be unreliable as a measure of absorption unless 
the conditions are very carefully controlled." 
Their results will none the less be reproduced in | 
Table I as the best at present available. It ' 


TABLE I. Decrement coefficient of the electrical reaction ex- 
pressed as a reciprocal number of wave-lengths. 


Reprinted from Kistiakowsky and Richards." 


Frequency Pinmm 
760 400 150 
10 -- 0.12 0.48 
42 0.04 0.12 _— 
80 0.06 0.10 0.40 


® Kistiakowsky and Richards, J. Am. Chem. Soc. 52, 
4661 (1930). 

8 Klein and Hershberger, Phys. Rev. 36, 1262 (1930). 
Hubbard, ibid. 41, 523 (1932). 
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is apparent that these coefficients can have little 
relation to h as defined above, since both at 10 kc 
and 80 kc and 150 mm they exceed the theoretical 
value of Apimax). It appears rather that they 
represent some result of the mass reaction of 
the gas on the oscillator, and that, therefore, the 
amplitude of the oscillator is not independent of 
the position of the reflector. It must be pointed 
out that the 10 kc and 42 kc measurements are 
further invalidated by a damping coefficieit 
which corresponds in magnitude to the tube- 
dispersion.’ It seems, therefore, unwise to 
attempt to form any conclusions from them. 

The observations of Teeter? are even less 
worthy of confidence. His electrical decrement 
coefficient at 25°C, 550 mm, and 54 kc is 0.82, 
which is some eight times as large as that 
reported by Kistiakowsky and Richards” at 
higher frequency and lower pressure (Table I 
above, 80 kc, 400 mm) and can, therefore, hardly 
bear any useful relation to the absorption 
coefficient. Teeter’s failure to obtain evidence of 
transmission of sound above 54 kc is strikingly in 
contradiction both with the prediction of theory 
and with the observations of previous ex 
perimenters.”: ! 

It is apparent that in both these investigations 
the comparison of measurements is made pre- 
carious by discontinuous changes in acoustical 
conditions as the pressure and frequency are 
changed. Each different oscillator introduces a 
new set of acoustical conditions, and each pres- 
sure change involves a significant alteration of 
the mass-reaction of the gas on the oscillator. A 
method of locating the absorption maximum by 
means free from these objections was, therefore, 
devised. By using the same oscillator and reflector 
for all measurements, and by holding the pressure 
constant, large changes of the type indicated 
above are eliminated. The gas may then be taken 
through its absorption maximum by varying the 
temperature, since the dissociation rate constant 
of nitrogen tetroxide varies with temperature in 
proportion to a factor exp —14,000/RT7. The 
density of the gas will, of course, vary inversely 
with the degree of dissociation, but over a small 
temperature range it seems reasonable to suppose 
this influence negligible. 

See Lamb, The Dynamical Theory of Sound, 2nd Ed., 
p. 196 (published by Arnold). 


If the rate constant of 4.8X10* sec.—' for 
nitrogen tetroxide at 25°C and 260 mm is 
accepted from the dispersion measurements, and 
the constants Ao, Bo, Co are taken as being those 
given in the previous section, wpqmax) will be 
about 470 kc at this temperature and pressure. 
This frequency must not be considered a pre- 
diction, since it is of course extremely sensitive to 
Co, and Cy has been arbitrarily given the value of 
10 cal. It suffices, however, to show that the 
450.1 kc oscillator already described (Fig. 4 and 
corresponding text of the work on dispersion, 
reference 1) is suitable to explore the absorption 
maximum at 260 mm in the neighborhood of 
25°C. At 30°C the rate constant is about 6.7 X 104 
and at 15°C about 2.6X10* at this pressure so 
that variation of the temperature between these 
limits should be sufficient to carry the gas through 
the absorption maximum provided that no heat 
capacity lag is present. The change in the 
absorption to be expected in this range is not, 
however, great, and it is difficult to predict 
whether measurements of the logarithmic decre- 
ment coefficient of the electrical reaction may be 
made with sufficient accuracy to detect it. 

Table II summarizes a series of measurements 
at 451 kc and 260 mm, using this method. The 
first column, which is headed x/\, is the distance 
from the reflector to the oscillator measured in 
wave-lengths. Thus the face of the oscillator 
corresponds to x/A=0, and the first, third, fifth, 
etc., points of maximum electrical reaction, 
which do not correspond to an integral number of 
wave-lengths, are not reported. The relative 
electrical reaction is expressed as the ratio 
between the microammeter deflection at a point 
x/k=n divided by that at the point x/A=1. 
Only 15 wave-lengths are recorded, although by 
suitable amplification nearly twice that number 
could be detected. | 

In argon, at a density approximating that of 
N2O, at 25°C and 260 mm, the electrical reaction 
of the 451 kc oscillator decreased to 62 percent of 
its initial value after 15 wave-lengths. 

That this attempt to locate the absorption 
maximum has been unsuccessful may be seen by 
looking horizontally across the table (for example 
at x/A=5, 10 and 15). No significant diminution 
of the relative electrical reaction is noticeable at 
any temperature, and it is apparent that an 
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TABLE II. Relative electrical reaction at various distances from the oscillator in nitrogen tetroxide at 260 mm, 451 kc, and 


various temperatures. 

x/r 17.0° 18.3° 19,2° 20.0° 21.0° 22.0° 23.0° 24.1° 25.1° 26.1° 27.0° 29.0° 
1 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 
: 0.64 0.61 0.66 0.72 0.76 0.64 0.60 0.59 0.68 0.73 0.64 0.63 
3 46 46 49 48 De A5 41 46 46 49 42 48 
4 32 .40 at .30 33 .28 33 .36 33 
5 Py .28 ag .27 .21 .24 .24 .20 
6 18 18 .20 21 19 .16 17 Ad .16 18 15 18 
8 «AZ .08 .08 10 .08 | | .09 10 
9 09 .06 09 .10 08 07 09 07 07 09 07 07 

10 .08 05 .08 .08 .07 05 .08 .06 05 .06 .07 .06 

11 .06 .03 .07 .06 05 .03 07 .04 05 05 05 
12 04 02 .06 05 04 03 .06 04 03 04 04 04 
13 .03 01 .05 .04 .04 .03 .04 .03 .03 .03 04 04 
14 .03 01 .04 .03 .03 .02 .03 03 01 .03 03 .02 
15 02 01 .02 .03 02 01 .03 02 01 .02 02 02 


accuracy which is greater by at least one order of 
magnitude is necessary to secure the desired 
information. The measurements have neverthe- 
less been reported, since such scant conclusions as 
they permit are considered distinctly favorable to 
the interpretation of the dispersion which has 
previously been given. 

It will be noticed that the intensity of the 
electrical reaction drops with fair uniformity to 
2 percent of its initial value after 14 wave-lengths. 
This is in quantitative agreement with the 
prediction of Luck based on the supposition that 
the dissociation reaction alone is failing to follow 
the sound wave. Although the electrical decre- 
ment coefficient'® and the acoustical absorption 
coefficient may be considered identical only 
under experimental conditions far more circum- 
spectly controlled than those here described, this 
agreement is considered encouraging. A dimi- 
nution to 2 percent after 14 wave-lengths is 
conspicuously less than that to be expected in 
case (c), if the dissociation reaction and all of 
the vibrational heat capacity are simultaneously 
failing in the sound wave, and Q=Qpc=1.14. 
The possibility of case (b) has been conclusively 
eliminated by the dispersion measurements, and 
need not be discussed. Table II appears to 
indicate, therefore, that case (a), where the 
dissociation reaction alone is producing the 
absorption, is that which approximately obtains 
in nitrogen tetroxide at 25°C and 450 kc. 
Furthermore, the fact that decrement of the 


°A single decrement coefficient is not quantitatively 
satisfactory for describing the results in Table II. 


electrical reaction is not greater may be taken as 
support for the already well grounded conclusion 
that the tube correction in this apparatus is 
small. Such tentative conclusions as are per- 
missible from Table II are, therefore, in full 
accord with those previously based on the far 
more quantitative dispersion measurements. 

At present the heat capacity lag constants 
Oxo, and @yx,o, must, therefore, be considered 
unknown.!® 

In closing the section which deals with ab- 
sorption the situation may perhaps be sum- 
marized. It appears that at present measure- 
ments of the dispersion of sound are both easier 
to perform and more trustworthy in their 
implications than measurements of the ab- 
sorption. It is, in fact, useless to attempt to 


16 Teeter (reference 2) has stated that the dispersion 
which we have observed in nitrogen tetroxide may in- 
telligibly be interpreted as wholly due to a heat capacity 
effect, and justifies this by a calculation in which a mean 
lifetime @ of 5X10- sec. is assigned to the NO.— NO, 
mixture in rough analogy to Kneser’s results on carbon 
dioxide. The calculation is made on the basis of an equation 
which is the dispersive analogue of (22) above, and is, 
therefore, applicable only when the dissociation reaction 
has been wholly eliminated from the adiabatic cycle. 
Since this cannot be the case under the conditions of our 
measurements, we disagree with his conclusion. Further- 
more, if 0=5 X10 as Teeter supposes possible the corre- 
sponding absorption maximum would be at about 3 kc 
at 30°C and 260 mm from Eq. (15) above. This is not con- 
sistent with his own conclusion that the absorption maxi- 
mum “‘comes somewhere near 700 kc.”” The basis for his 
conclusion is not given, and since neither temperature nor 
pressure are specified it cannot be compared with any 
deductions which are here supported. 
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obtain quantitative information concerning the 
absorption with the apparatus at our disposal in 
this laboratory. Nevertheless, it is to be hoped 
that some investigator experienced in absorption 
measurements, who has at his disposal a large 
number of calibrated quartz crystals, will ulti- 
mately bring this research to a positive con- 
clusion. The position and magnitude of the 


absorption maximum can provide information of 


great theoretical interest concerning heat ca- 
pacity effects in dissociating gases. This informa- 
tion is more difficult to obtain from dispersion 
measurements, for the reason that the position 
and magnitude of the absorption maximum 
indicate the value of V,,, whereas dispersion 
measurements can obtain V,, only at much 
higher frequencies by direct measurement. If the 
dispersive region is compounded of two or more 
concomittant types of time lag the terminal 
velocity at high frequencies corresponding to 
each separate process may be clearly deducible 
from the absorption maximum although wholly 
inaccessible to direct measurement. 


Part II. SUPPLEMENTARY DATA CONCERNING 


THE DISPERSION MEASUREMENTS 
PREVIOUSLY REPORTED 


Further discussion of the tube-correction 


The application of a tube-correction deter- 
mined by experiment considerably increased the 
accuracy of the measurements of dispersion 
which were previously reported. We feel that the 
method by which this correction was obtained 
has already been sufficiently discussed, but that 
the account of its application is too brief to be 
immediately intelligible. Further, parts of Table 
III and all of Table IV were inadvertently 
printed with the measurements uncorrected, in 
direct contradiction to what was stated in the 
text. This has caused an obscurity concerning the 
tube correction which is wholly unnecessary” 
and we welcome an opportunity to dispel it. 
It must, however, be emphasized that nothing 
will be written in this section which was not 
previously considered in calculating the rate 
constants. 


" Teeter (reference 2) has pointed out the inconsistencies 
in Tables III and IV, and has correctly estimated their 
origin. 


TABLE III. 


Linear Linear 
~— Velocity frequency Velocity AV 
Cc kc) ( 


m-sec.~!) m-sec.~') (m-sec.~') 
1.0 9.265 170.8 92.07 172.8 1.8 
10 9.341 171.4 92.14 173.0 
7.0 9.343 176.4 92.06 177.2 0.8 
8.0 92.06 178.3 0.8 
8.1 9.340 177.5 : 

10.0 9.331 179.2 92.02 179.9 0.7 

13.9 9.336 183.05 92.00 183.6 0.55 


25.0 9.298 195.3 91.89 195.7 03 
25.0 9.298 194.8 91.89 195.0 . 


30.0 9.213 202.4 91.73 202.6 0.2 
30.0 94.00 202.1 450.7 205.6 in 


94.00 202.0 450.7 205.5 
94.00 202.2 450.7 205.4 


TABLE IV. 


P Frequency Velocity Frequency Velocity AV 
(mm) (kc) (m-sec.-!) (kc) m-sec.~!) (m-sec.~') 


30°C 
342 94.00 198.4 450.7 200.7 2.3 
260 94.00 202.1 450.7 205.4* 3.4 
235 94.00 203.4 450.7 207.9 4.5 


25°C 
668 9.301 185.2 91.89 185.3 0.1 
260 9.301 195.1 91.89 195.4 0.3 
162 9.301 202.0 91.89 202.5 0.5 
1.0°C 
260 9.340 171.1 92.10 172.9 1.8 
156 9.340 175.4 92.10 179.1 3.7 
132 = 9.340 176.7 92.10 181.0 4.3 


* This number is the average of the three corresponding 
measurements in Table III and should evidently be printed 
205.5. It has not been corrected in order to facilitate 
comparison with the previous printing. 


In order to be consistent with the text of the 
previous work, Tables III and IV should have 
been printed as they are here given. A careful 
comparison with the previous printing will show 
that the tube-corrections have been applied in 
the manner stated to all measurements in which 
the difference between the two velocities quoted 
did not accord with the sixth column and that, 
with the exception of a single misprint, the 
agreement is now complete and the sixth column 
unaltered. One further change has been made in 
reprinting the tables in that the sixth column 
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has been headed AV to show that it represents 
the difference in velocities at the two appropriate 
frequencies instead of V.—Vo as previously. 
These quantities may, of course, be quite differ- 
ent, but in the tables happen to be the same 
within experimental error. 

Before describing in further detail the calcu- 
lation of the rate constants from these measure- 
ments it should be emphasized that at 1°C and 
132 mm (Table IV) AV is 4.3 m-sec.—', and the 
velocity at 9 kc, 176.7 m-sec.~!. At 7°C and 260 
mm (Table III), on the other hand, the velocity 
is almost identical at 9 kc and AV is only 0.8 
m-sec.-!, These two measurements alone are 
considered sufficient refutation of any funda- 
mental objections to our conclusions based on the 
tube-correction, or a possible major effect due to 
heat-capacity dispersion. It was not previously 
considered necessary to emphasize this point. 

The calculation of rate-constants from these 
measurements will now be considered in greater 
detail. It has already been stated that the chief 
cause for numerical disagreement between the 
measured velocity of sound and that calculated 
from Einstein’s theory lies in uncertainty con- 
cerning the density of the gaseous mixture. A 
method of calculation was, therefore, employed 
in which the pressure and density were wholly 
eliminated by the use of the velocity of sound at 
very low frequency. This was effected by writing, 
in the notation of this paper, 


ky = (wCo/Bo) { (Vn? / (Vu? } # 


where the symbol v now refers to an experi- 
mentally determined phase velocity. V,.2 was 
calculated from the experimental velocity by 
means of the relation 


V..2 Bo(Cot+R)/Co(Ao+Bo)}. 


This procedure is applicable quite generally to 
measurements on the dispersion of sound. Not 
only does it very largely remove systematic 
errors because of imperfections of the gas, in- 
Correct determination of its thermodynamic con- 
stants, etc., but it makes possible the calculation 
of reasonably accurate time lag and dissociation 
rate constants when the measured values of the 
velocity accord relatively, but are greatly differ- 


ent from their values calculated from theoretical 
considerations alone. 

For this reason no great concern was felt when 
the measured and calculated velocities differed 
by as much as 2 m-sec.~'. The tube corrections 
were applied merely to make AV obtained from 
Vo2— Vo a self-consistent quantity, and similarly 
to make comparable the quantity V4;1— Vos. The 
method of calculation made it unnecessary that 
Vo4 and Vo should be comparable, or that either 
should agree with the calculated velocity to 
anything like the precision claimed for AV. It 
will be remarked, for instance, that at 30°C and 
260 mm Voo is 202.6 m-sec.~!, and Vou, 202.1 
m-sec.—'. A difference of this kind is to be ex- 
pected, for no effort has been made to apply a 
tube-correction which will make the two com- 
parable. It is further a matter of indifference 
which of the two is accepted provided that the 
calculation of the rate constant is consistently 
made. There are two reasonable methods of mak- 
ing the calculation, (A) to add 0.5 m-sec.~' to Vo4 
and to V4; to make them comparable to V», and 
Vo2, or (B) to subtract 0.2 m-sec.~ from Vo4 to 
account for the dispersion and to use this for Vo. 
The result obtained is the same within experi- 
mental error 


wCo we 
By | (206.0)?— (202.4)? 
wCo( (211.0)2—(205.5)2) wo 
(B) | =1.250— 
By |(205.5)?— (201.9)? 0 


Nevertheless an attempt was made to ascertain 
whether Vo2 or Vo, was the more worthy of 
confidence by measuring the velocity of sound in 
argon, a quantity which can be calculated with 
some precision since Cp/C, = 1.667 on theoretical 
grounds. On this basis V9: appeared the more 
accurate, as already reported. 

One further aspect of this argument which 
appears to have caused confusion is the determi- 
nation of the difference Vs— Vo, because of the 
entry of the gas into the dispersive region. This 
was shown negligible in all cases by a calculation 
at 1°C and 260 mm. Here V,, =177.9, Vo2= 172.9, 
Vo=171.1 m-sec.—'. is, therefore, approxi- 
mately 
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2n(9.2X104)Co 


Bo (172.9)?— (171.1)? 
=7.0X 10%. 


In order to find if Vs is appreciably different from 
V, the calculation is then made 


7.0X 108 = 


Bo Vo? 


and it is found that V;— Vo=0.03 m-sec.~! and 
therefore inconsiderable. It was on this basis that 
we have written, ‘The accuracy of the rate 
constants in Table VI must now be considered. 
It may be shown from these by (9)”’ [the equa- 
tion used above ] “that the velocity of sound at 
very low frequency at 1°C is less than 0.1 m-sec.— 
below that quoted at 9 kc. Error from this cause 
is, therefore, eliminated.” 

One further point concerning the tube-correc- 
tion requires more extended discussion. We have 
suggested that the tube-correction may undergo a 
change in the dispersive region of frequency, 
since a literal application of the Helmholtz- 
Kirchhoff expression for the damping of plane 
waves in small tubes appears to indicate a 
certain increase in damping as the compressibility 
ratio of the gas is increased. We were content to 
show that this effect, if it occurred, would be 
_ negligible. Teeter has, however, based a criticism 
of our work on this subject, and the argument 
must, therefore, be more fully stated. 

For this purpose it will be sufficient to use a 
very approximate statement of the Helmholtz- 
Kirchhoff relation. The tube correction for a 
particular diameter of tube, a particular gas, and 
a particular frequency will be called 7 when 
V..—Vo=+0. If the conditions in the gas are 
changed so that V,—V» becomes appreciable 
(for example by lowering the temperature or 
pressure) the Helmholtz-Kirchhoff expression 
leads to the expectation that a larger tube- 
correction 7,’ would be found.'® A much sim- 
plified approximation 


Te /Te = Vo? 


18 Since it is desired only to find the maximum error due 
to this change the legitimacy of such a literal interpretation 
of the expression for the tube-dispersion need not be 
discussed, 
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reproduces the numerical implications of the 
Helmholtz-Kirchhoff expression with quite suffi- 
cient accuracy for the present purpose. Since 
V..2/ Vo? might change from 1.00 to 1.25 from the 
beginning to the end of the dispersive region an 
increase of 0.5 m-sec.~! might occur in a tube- 
correction of 2.0 m-+sec.~!. Since V,,— Vo is about 
10 m-sec.~! this represents a 5 percent change in 
the differences in velocity due to the tube- 
correction. With a tube-correction of 2 m-sec.~', 
then, the tube-correction will be approximately 
Tot+1.05(V.— Vo), and will equal 7, for small 
values of V..— Vo. 

This estimate of the error may easily be shown, 
however, to be unnecessarily generous when 
applied to our measurements. In the tubes which 
we used the tube-correction was always less than 
2 m:sec.~! at high frequencies, and only at high 
frequencies did V.,.—Vo become of important 
magnitude. At lower frequencies the tube- 
correction was great, but V.,— Vo small enough to 
be negligible. An estimate of the error possibly 
introduced by change in the tube-correction is 
appended in Table V. V.2/V? from appropriate 


TABLE V. Estimate of error introduced into previous measure- 
ments on dispersion by a possible change in the 
tube-correction. 


Maximum 


Conditions 


9ke 1°C 260 mm 1.001 0.8 y 
92 kc 1°C 260 mm 1.002 2.0 0.004 
94 kc 30°C 260 mm 1.010 0.4 0.004 
451 kc 30°C 260 mm 1.050 0.6 0.030 


measurements gives the column headed 7,'/7.. 
The column headed 27, is the most probable 
value of the tube-correction when V,,—Vo=0 
multiplied by two for good measure. The last 
column shows twice the error probably intro- 
duced by neglecting the possible change in the 
tube-correction in the dispersive region, and has 
been called the maximum error. 

It is perhaps sufficient comment on this table 
to repeat that our claimed error was never less 
than 0.2 m-sec.—'. We can therefore only reiterate 
our previous statement that the error introduced 


by neglecting the tube-correction change “is 


4 
| 
error 
& 2 — Tw) 
Te! * Te m-sec.-! m-sec.! 
4 «x 
a 
| 
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believed to be considerably inferior to that due to 
other causes.” We do not understand on what 
grounds Teeter has disregarded this statement 
without attempting to refute it. 

In concluding this section it should be noted 
that nothing which has been said above must be 
considered as a fortiori justification of our 
previous work, for it is an amplification of an 
argument already indicated. Neither must it be 
considered a defense, for our previous position 
has not been attacked, but wholly disregarded. 
We must express our regret that it has been 
necessary to discuss in such intimate detail 
matters of very small scientific interest. 


The effect of deviations from ideal gas laws and 
of absorption on rate constants computed 
from dispersion measurements 


We had previously supposed that the approxi- 
mate agreement of V» and of dV)/dP at 260 mm 
with the values obtained from Einstein’s equa- 
tion’? combined with the method we have used 
in calculating the rate constants, removed any 
serious error due to non-ideality of the mixed 
gas. Teeter’s demonstration, from the argument 
given by Luck, that this is not- necessarily the 
case, is therefore welcome. His estimate of the 
error introduced at the absorption maximum is, 
however, conclusive only if Q (=V,,/Vo) is a 
measured quantity. When, as in the case we have 
discussed, V,, is calculated, k; at the absorption 
maximum becomes from (3) and (15) 


ky =w{ Co(Co+R)/Bo(Aot+Bo)}! 


and A» is also involved. Luck has already 
pointed out that the coefficients involved in the 
calculation of Ao from his development of the 
theory cannot at present be estimated. It might, 
of course, be argued that the experimental 
agreement of V» with the result calculated when 
using ideal gases indicates that A» and Bo are 
affected equally by departures from ideal be- 
havior. This would introduce, at the absorption 


* This is well illustrated at 25°C by Teeter’s Fig. 5. It 
should perhaps be pointed out that our measurement at 
9 ke and 668 mm has been incorrectly plotted in this 
figure, although it is correctly quoted in the corresponding 
table. The error is minor, but it gives an unmerited appear- 
ance of incoherence to our work. 


maximum, about the error which Teeter has 
estimated, but would involve as a corollary our 
original assumption that rate constants calcu- 
lated from small values of V..— Vo are without 
error due to this cause. It appears impossible, at 
present, to give this question a more satisfactory 
answer. This is especially the case because 
experimental determinations of the density are 
lacking below 25°C and, even above this temper- 
ature, do not define the value of P/p with 
accuracy sufficient for acoustical purposes. 

Teeter has pointed out that errors in rate 
constants due to departure from ideal behavior 
are of opposite sign to those due to the neglect of 
the effect of absorption on the velocity. He has 
estimated the latter effect variously as 20 percent 
and as 4 percent.”° We agree with the latter of 
these two estimates. In the only case of quanti- 
tative importance for our argument, that at 30°C, 
451 kc and 260 mn, this effect may readily be 
calculated. Here Q, when using the average value 
from the “maximum” and “minimum” heat 
capacity, is about 1.048, Cy is 9.83 and Bo (using 
a=0.359, D=14.0 kg cal-mole) is 661. This 
gives h=0.267 from (14). Then 


Vo -{ 1—(h/47)? 


Vo 

where v, is the measured velocity of sound 
(27v.=wd) and V, is the square root of the real 
part of +/A. This equation is of course derived 
from (9), and was first suggested by Luck.’ The 
measured velocity of sound is, therefore, about 
0.14 m-sec.~! below that given by r.p. +/A under 
these conditions. The relation of the true rate 
constant k, to that calculated neglecting the ab- 
sorption ke is then 


= 0.99933 


ky Vo? 
V..2— Vo? 


—= = 1.04 
kr 
which agrees with Teeter’s estimate. This is well 
within the experimental error which we have 
allowed, and is probably overbalanced by the 
error due to the departure of the gas from ideal 


behavior. 


20 Teeter, reference 2 and J. Am. Chem. Soc. 54, 4111 
(1932). 
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There appears to be no reason at present for 
revision of the rate constants which we have 
given previously from either of the causes dis- 
cussed in this section. 


Corrections of the previous communication 


Unfortunately, a number of misprints and one 
actual error appear in the previous work. 
The following corrections should be noted. The 


argument was developed throughout by using 
the correct and not the incorrect statements. 

Page 115. Lines 13 and 14, left-hand column. 
The real part of the two exponents should be 
— Bx and not —8. 

Page 116. Lines under Eq. (8) should read ‘‘y is 
the ratio of the isothermal to the adiabatic 
compressibilities.’’*! 

Page 117. Eq. (9) should read” 


The calculations of the rate constants were 
made using the convenient grouping of constants 


printed by Gruneisen and Goens,” after a 


comparison with a properly reduced version of 
our (8). It may easily be verified that the rate 
constants we have given are free from the error 
involved in (9) as originally printed. 

Page 117. Fig. 2 should have as abscissae 
co 4567 0 making the point of inflection of the 
w= 10° curve nearly at 

Page 119, 22nd line from top, right-hand 
column, should read “pure plane ‘sinusoidal 
wave.” 

Page 123. Tables III and IV should read as 
printed above. 

Page 123, lines 17-15 from bottom, right-hand 
column. The statement ‘The second equation of 
(3) neglects the influence of temperature and 
pressure on the degree of dissociation” is non- 
sense, and should be deleted. 

Page 125, line 10, right-hand column, should 


read 


and it should be emphasized that V» refers to the 
experimental velocity. 


21 Teeter has pointed out this slip. By stating that we 
have the ratio inverted he is apparently, however, referring 
to (8) which is correctly printed. His statements that 
y=Cp/C, and that it is a complex have no relation to 
anything we have said. It is necessary to call attention to 
this since, as he has not reprinted y, the impression is 
distinctly misleading. The relation between Einstein’s 6 
and Luck’s / is of course 6] =1. 


Page 127. The formula for 1/A should omit J/, 
from the denominator of the last term on the 
right. 

The only systematic error which Teeter has 
pointed out must now be discussed. We have 
incorrectly used 14.6 kg-cal. as the heat of 
dissociation at constant volume derived from the 
measurements of Verhoek and Daniels* instead 
of the correct value of 14.0 kg-cal. The effect of 
this error is to make the rate constants previously 
reported uniformly about 9 percent too low. It is 
not surprising that this figure falls within the 
estimated experimental error of 10 percent be- 
cause uncertainty concerning the dissociation 
heat was of course included in the estimate. 
Evidently it is necessary to consider that the 
dispersion measurements indicate a reaction rate 
constant of 5.3104 sec.-! at 260 mm and 25°C, 
and similarly to raise by 9 percent all the rate 
constants quoted. The exact absolute values of 
these constants have, however, little contempo- 
rary importance. It is more significant that the 
activation energy which was previously given 
(13.9+0.9 kg-cal.) is unchanged, and is now so 
near to the most probable value of the dis- 
sociation heat (14.0 kg-cal.) that the agreement 
between the two must be considered in part 
fortuitous. 


* We are indebted to Teeter for calling attention to this 
error. 

23 Gruneisen and Goens, Ann. d. Physik 72, 193 (1923). 

* Verhoek and Daniels, J. Am. Chem. Soc. 53, 1250 
(1931). 
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LETTERS TO THE EDITOR 


This section will accept reports of new work, provided these 
are terse and contain few figures, and especially few halftone 
cuts. The Editorial Board will not hold itself responsible for 
opinions expressed by the correspondents. Contributions to 
this section must reach the office of the Managing Editor 


not later than the 15th of the month preceding that of the 
issue in which the letter is to appear. No proof will be sent 
to the authors. The usual publication charge ($3.00 per 
page) will not be made and no reprints will be furnished free. 


Spectroscopic Evidence for the Molecule PN 


In a discharge through air in a tube which previously had 
been used with phosphorus we found a rather extensive 
system of simple diatomic bands in the region 2400—2900A 
which had not been recorded before. Further experiments 
showed that these new bands do not occur in a similar 
discharge through air or nitrogen without the presence of 
phosphorus. On the other hand they have never been ob- 
served in discharges through phosphorus vapor alone. 
The bands do however occur strongly in a discharge 
through a mixture of pure nitrogen and pure phosphorus 
vapor (at the same time strong clean-up occurs). We 
therefore conclude that the band system is due to the 
diatomic molecule PN. This is confirmed by a closer 
analysis of the band system (see below). ° 

A search through the literature showed that the only 
well-known compound of phosphorus and nitrogen is 
P;N;, discovered by Stock. In addition Moldenhauer! 
found several years ago that in a discharge through nitrogen 
and phosphorus a yellow powder is formed which contains 
equal amounts of phosphorus and nitrogen and which is 
remarkably stable. The free molecule PN in the vapor 
phase has not been observed up until now. The above 
experiments show that this molecule at least temporarily 
exists in gas discharges through a mixture of Nz and Py. 

We have so far measured the heads of 24 bands. These 
are represented rather accurately by the formula 


»=39,699.3 + (1094.75v’ —7.25v””) — (1329.380"” — 6.980"). 


The fine structure of the bands is fairly well resolved on 
our plates (second order of 3 m grating). Each band con- 


sists of one P, one Q and one R branch, thus the system is a 
1]]-—'Z transition (most probably the lower state is '2). 
Therefore the remote possibility that the system is due to 
PO or CP, which might have been present as impurities, 
is excluded, because these can only have doublet (or 
quartet) terms. A preliminary rotational analysis gave the 
constants: 


=0.78 cm, Bo’ =0.73 cm=; =1.49A, =1.54A. 


These, as well as the vibrational constants (see the above 
formula) are similar to those for the corresponding systems 
of CS and SiO, as is to be expected. This is a further 
argument for PN as carrier of the band system. 

It seems probable that the lower electronic level 'D of the 
bands is the ground state of the molecule. The energy of dis- 
sociation of this state comes out by linear extrapolation to be 
7.8 volts. This of course is to be regarded as an upper limit. 

A full report of this work will appear shortly elsewhere. 
A further investigation of this molecule, both in emission 
and absorption, seems promising especially with respect 
to a comparison with Nz and P». 

J. Curry 
L. HERZBERG 
G. HERZBERG 
Physikalisches Institut, 
Technische Hochschule, 
Darmstadt, Germany, 
August 19, 1933. 


1W. Moldenhauer, Chem. Ber. 59, 926 (1926). 


1-4 Addition and Subtraction Reactions in Hydrocarbon Chemistry' 


Eyring, Sherman, and Kimball? have treated 1-4 addi- 
tion to conjugate double bonds by means of the model 


cc c:¢ 


c 
H:H “HOH” 


where each dot represents an electron considered in the 
calculations. Using 74.5, 92, and 102.4 kcal. for the 
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strengths of the C—C, C—H, and H—H bonds, they 
find an activation energy of 64 kcal. for the addition of 
hydrogen to conjugate double bonds. These bond strengths, 
however, give a heat of reaction of only 7.1 kcal., whereas 
heats of combustion require 32 kcal. This difference must 


1 Published by permission of the Director, U. S. Bureau 
of Mines. (Not subject to copyright.) 

2 Eyring, Sherman and Kimball, J. Chem. Phys. 1, 
586 (1933). 
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be made up mainly by decreasing the strength of C—C, 
which denotes actually the difference between C=C and 
C—C. The calculated heat of activation then will be some- 
what lower, probably about 50 kcal. Sherman and Eyring,* 
using a more suitable value for C=C, previously calculated 
an activation energy of 51.5 kcal. for the hydrogenation of 
ethylene. Here again the reaction is approximately 30 
kcal. exothermic. The theoretical activation energy for 
either 1-2 or 1-4 dehydrogenation is thus 80 kcal.; the 
experimental value in the former case is 70 kcal., and a 
similar figure seems reasonable for the 1-4 reaction. 

There is also a second possible type of 1-4 dehydro- 
genation:! 

CH: CH:+CH: CHa. 


This type actually corresponds somewhat more closely to 
the six-electron model than does the production of 
butadiene from butylene. It is approximately 20 kcal. 
more endothermic, however, since it takes more energy to 
break C—C than it does to reduce C=C to C—C. The 
difference in activation energy will be considerably less than 
20 kcal.; accurate predictions are impossible, but it appears 
worth while to look for this type of decomposition of 
hydrocarbons with a straight chain of at least four carbon 
atoms. By analogy with the 1-2 case, we shall expect to 
find 1-4 loss of methane, ethane, etc., and also, when the 
parent hydrocarbon is large enough. 

The recent work of Frey and Hepp*® may be interpreted 
as showing that this type of reaction occurs. For n-butane 
itself the evidence is somewhat conflicting; 1-4 dehydro- 
genation may constitute up to 5 percent of the primary 
reaction. For n-pentane the data may be interpreted as 
indicating the following primary products. 


Percent 


18.5 
12.2 


Percent 


H2+CsHwo...... 5.2 
CHs+C.Hs..... 30.9 
CsHe+CsHe.... 31.5 
CsHs+CeHy.... 1.7 


69.3 30.7 


For n-hexane, 60 percent of the reaction seems to go by 
way of 1-4 unsaturation, and only 40 percent by 1-2. 
For 2-methylbutane and 2,3-dimethylbutane, 20 percent 
and 12 percent, respectively, of the total reaction involves 
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1-4 unsaturation. These results have been interpreted by 
Frey and Hepp, and by Rice,’ in terms of Rice’s free 
radical mechanism. A detailed analysis shows that the 
free radical theory and the present unsaturation theory 
predict the same products in approximately the same 
yields. For very long chains it would be necessary to 
postulate 1-6, 1-8, etc., unsaturation, with resulting 
scission into four, five, or more primary products, to match 
exactly the free radical predictions; there is no experi- 
mental evidence on this point. The free radical mechanism 
does predict a decrease in the mean molecular weight of 
the primary products as the pressure is lowered. On the 
unsaturation theory such a decrease should not occur in 
the region where the reaction rate is strictly first order. 
No existing data even approach being good enough to 
make this test. 

There is little doubt that free radicals are formed from 
metal alkyls, and in a few other cases where there are un- 
usually weak bonds, but it appears possible to give a 
nearly complete account of the decomposition reactions of 
organic chemistry in terms of 1-1, 1-2, and 1-4 unsatura- 
tion. It is hoped to discuss this subject in more detail later. 
For the present, it is sufficient to mention a few further 
examples of 1-4 unsaturation: the formation of ethylene 
and butadiene from cyclohexene, of carbon dioxide and 
allyl alcohol from glyceryl oxalate, and of methane, 
ethylene, and formaldehyde or hydrogen, ethylene, and 
acetaldehyde from ethyl ether. In this last case, the chief 
reaction is 1-2 splitting to ethane and acetaldehyde. 

Louis S. Kasset’ 

Physical Chemistry Section, 

Pittsburgh Experiment Station, 
U. S. Bureau of Mines, 
Pittsburgh, Pennsylvania, 
August 24, 1933. 


3Sherman and Eyring, J. Am. Chem. Soc. 54, 2561 
(1932). 

4Cf. Gault and Hessel, Ann. chim. [10] 2, 376 (1924), 
who postulated a similar process to account for the pro- 
duction of hydrogen in the pyrolysis of hexadecane. 

5 Frey and Hepp, Ind. Eng. Chem. 25, 441 (1933). 

6 F. O. Rice, J. Am. Chem. Soc. 55, 3035 (1933). 

7 Associate physical chemist, U. S. Bureau of Mines, 
Pittsburgh Experiment Station, Pittsburgh, Pa. 


On the Concentration of Deuterium by Electrolysis 


Bleakney and Gould! have raised the question as to 
the reason for the low concentration of deuterium in 
commercial electrolytic cells as is shown by the analysis of 
hydrogen obtainable commercially. They point out that 
this low concentration would imply that this tank hydrogen 
has all come from apparently new cells. 

The specific gravity of water which was distilled from 
electrolyte which was kindly furnished by the Pittsburgh 
plant of the Burdett Hydrogen and Oxygen Company has 
been obtained. This water was carefully purified by 
distillation through a still which is used for the production 
of conductivity water, and then deaerated. The specific 
gravity was obtained by means of a buoyancy balance 


similar to that described by Lewis and Macdonald.’ Six 
cells were sampled and the specific gravity obtained was 
1.000013 compared to ordinary conductivity water a 
18.2°C. The results on all samples agreed to about one 
part in a million. 

The cells from which the electrolyte was taken had been 
in operation for about five years without the electrolyte 
having been renewed. A simple calculation showed that 


1 Walker Bleakney and Austin J. Gould, Phys. Rev. 4: 
268 (1933). 

2G. N. Lewis and R. T. Macdonald, J. Chem. Phys. !. 
341 (1933). 
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the amount of water electrolyzed during this period of 
operation was about 30 times the amount the cells con- 
tained. It would appear from this calculation that these 
cells had been in operation for a long enough period of 
time to have nearly reached their equilibrium concentra- 
tion, this equilibrium concentration being given by the 
relation deduced by Lewis and Macdonald? as 1/a times 
the concentration of the water added to the cells, where a 
is the separation factor. This separation factor a must be 
a function of the temperature at which the electrolysis is 
taking place, as is shown by the following calculation. 

If we assume Dalton’s law of partial pressures to apply, 
the fractional part of a mole of water vapor removed per 
mole of water electrolyzed will be 3Pw/2(Pa—Pw), 
which quantity we shall call «. Where Pw is the partial 
pressure of water vapor at the temperature 7, and Pa is 
the atmospheric pressure. Thus the amount of water re- 
moved from the cell will be 1+ moles per mole electro- 
lyzed. If 6; and 62 represent the amounts of H' and H? 
evolved per mole of hydrogen produced, and C, and C; 
are the respective concentrations of the two kinds of 
hydrogen in the electrolyte expressed as mole fractions, 
then 52C;/52.C2=c is the true electrolytic separation factor. 
The hydrogen lost by vaporization of the water will be « 


_ moles per mole electrolyzed. Thus the a of the Lewis and 


Macdonald? equation obviously becomes, 
a= C2. 


Assuming Bleakney and Gould's! value (1 part in 5000) 
for the abundance of deuterium in ordinary water, and 
using Lewis and Macdonald's’ value for the density of 
pure deuterium oxide, the value of o can be calculated if 
we further assumed that the electrolytic hydrogen 
analyzed by Bleakney and Gould! came from a typical 
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old cell of the type which was referred to above. This 
calculation gives a value of 0.12 for « which would agree 
fairly well with the provisional value of 0.10 given by 
Urey.’ By taking this value (0.12) for o and calculating a 
for this cell, (T=160°F) the value of @ obtained is 0.37. 
The experimental value of @ for such a cell will be con- 
siderably higher than this calculated value, due to the 
fact that evaporation from the surface of the electrolyte, 
which will be considerable at this temperature, has not 
been taken into account in the calculation. Calculating for 
Lewis and Macdonald’s? cells (T=35°C) the value of a@ 
obtained is 0.16 which is in fair agreement with the experi- 
mental value of 0.20 when it is considered that the 
mechanical loss of water during the electrolysis is neglected, 
which loss would tend to increase the calculated value 
of a. 

From the above considerations it would appear that 
commercial cells which are run hot in order to attain high 
power efficiency must necessarily have a low equilibrium 
concentration of deuterium, and also that the hydrogen 
produced by such cells will have a correspondingly low 
concentration of the heavier isotope. 

The author is indebted to Dr. J. G. Aston of the Chem- 
istry Department for the use of the thermostat and con- 
ductivity still which made the determinations of specific 
gravity possible. 

D. H. RANK 

Physics Laboratory, 

Pennsylvania State College, 
August 29, 1933. 


3G. N. Lewis and R. T. Macdonald, J. Am. Chem. Soc. 
55, 3057 (1933). 
4H. C. Urey, Rev. Sci. Inst. 4, 425 (1933). 


The Effect of Water Containing the Isotope of Hydrogen upon Fresh Water Organisms 


Since the announcement by Washburn and Urey! that 
in electrolysis the residual water is enriched in heavy 
hydrogen, there has been a greatly increased interest in 
the separation and properties of isotopes. That the chemical 
reactivity of two isotopes with such different atomic 
weights would differ very considerably was of course clear, 
and the opinion that heavy water would have striking 
biological effects was generally held. 

Lewis? has recently demonstrated that the seeds of 
tobacco will not germinate in pure H*H2O whereas seeds 
in ordinary distilled water sprouted within two days. We 
have recently performed a series of experiments to test 
the action of heavy water upon the vital activities of fresh 
water organisms. Our results are briefly summarized as 
follows: 

(1) Water containing the isotope of hydrogen in high 
Concentration (92 percent) was found to be toxic for the 
animals tested. 

(a) Tadpoles of the green frog Rana clamitans died 
within an hour after being placed in the water. Control 
animals in 30 percent heavy water, and in ordinary 
distilled water, remained unaffected at the end of twenty- 
four hours, 


(b) The common aquarium fish Lebistes reticulatus is 


killed by the water within an interval of two hours. 
Neither 30 percent heavy water, or distilled water exerted 
an effect over a twenty-four hour interval. 

(c) Flatworms, planaria maculata are destroyed within 
three hours in the 92 percent heavy water, but were 
uninjured by the 30 percent heavy water over a period of 
three days. 

(d) The protozoan, paramoecium caudatum is killed by 
the 92 percent heavy water within 48 hours. Water 
containing 15-20 percent of the heavy isotope of hydrogen 
did not affect the organisms during a three day interval. 
Paramoecium showed more resistance to the toxic action 
of the water than the highly organized animals studied. 

H. S. TAYLor 
W. W. SWINGLE 
H. EyrinG 
A. H. Frost 
Chemical and Biological Laboratories, 
Princeton University, 
September 15, 1933. 


1 Washburn and Urey, Proc. Nat. Acad. Sci. 18, 496 


(1932). 
2G. N. Lewis, J. Am. Chem. Soc. 55, 3503 (1933). 
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